Nuclear fusion device building earthquake spectrum automatic fitting and structure checking system

Abstract

The application discloses a kind of nuclear fusion device construction earthquake spectrum automatic fitting and structure checking system, comprising: GB parameter library and index module, for using index key retrieval and return specified data domain composition parameter set;Spectrum generation and fitting module, according to internal data of parameter set, call GB database to extract the parameter required for seismic design, automatically truncate negative value and control spectrum tail change consistent, generate the horizontal seismic influence coefficient curve after fitting, vertical seismic influence coefficient curve and corresponding frequency acceleration spectrum;Static equivalent and FRS determination module, the first natural frequency is obtained in fitting data by modal analysis and is judged according to preset threshold, respectively corresponding to static equivalent acceleration and response spectrum analysis operation;Data export and interface module, for output readable load spectrum file.The application can realize the automatic generation export operation of seismic influence curve, frequency acceleration spectrum, and support the seismic spectrum analysis and static strength checking of nuclear fusion device.

Description

TECHNICAL FIELD

[0001] The application relates to the field of nuclear fusion device building anti-seismic technology, and particularly relates to a nuclear fusion device building earthquake spectrum automatic fitting and structure checking system based on Matlab. BACKGROUND

[0002] With the growing global demand for clean energy, nuclear fusion can be a potential form of efficient and clean energy in the future. As a new type of clean energy device, the core structure of a nuclear fusion device usually includes a vacuum chamber, a superconducting magnet, a central column support system and peripheral building facilities. These structures need to withstand complex static and dynamic loads during the service period, and the earthquake action is one of the important working conditions that must be considered. In order to ensure the safety and reliability of the device and its key components under earthquake conditions, the earthquake response spectrum fitting and structure strength checking need to be carried out in accordance with the relevant national standards.

[0003] In the current engineering practice, the earthquake action analysis and structure checking for nuclear fusion devices mainly rely on the method of manual query and artificial fitting based on GB standards: first, the engineers manually consult the appendix tables and influence coefficient curves in the "Code for Seismic Design of Buildings", relevant nuclear power plant standards and other specifications according to the geographical location of the device, the structure type and the design earthquake grouping, and then modify them combined with the damping ratio. Afterwards, these data need to be manually or with the help of general software to be curve fitted and the earthquake influence curve is drawn. These standards provide earthquake influence coefficient curves and calculation methods under different seismic fortification intensities, site categories and structure damping ratios.

[0004] However, the following problems still exist in the actual use process:

[0005] (1) The data interface is not unified, which increases the secondary conversion work and causes the problem of earthquake parameter acquisition automation

[0006] The existing method needs to manually consult GB standards and manually input parameters, which is low in efficiency and prone to errors. The application aims to realize the automatic extraction of parameters such as device seismic fortification intensity, design earthquake grouping, site category, structure damping ratio, maximum horizontal earthquake influence coefficient and characteristic period by establishing a GB database interface.

[0007] (2) The manual table lookup and curve fitting are low in efficiency and redundant in process, which leads to the extension of project cycle and the decline of efficiency

[0008] The engineers often need to manually consult the tables and curves in the appendix of GB standards according to the location of the device and the design grouping, and then draw the earthquake influence curve combined with the correction formula of damping ratio. This process is tedious and prone to human error, which is difficult to meet the requirements of high efficiency and accuracy for the structure checking of complex nuclear fusion devices.

[0009] (3) The generation of the seismic spectrum is disconnected from the structural analysis, and it is difficult to realize full-process automation and error control

[0010] In the existing process, there is no direct interface between the generated seismic influence curve and the structural analysis software, and usually manual secondary data processing is required to obtain the seismic acceleration spectrum that can be used for spectral analysis. This method is not only inefficient, but also may introduce data bias.

[0011] (4) Static equivalent analysis lacks automation support

[0012] For different components in the modal analysis results, the first natural frequency is higher or lower than 30Hz, and the specification requires different seismic action calculation methods. However, the existing process relies on manual judgment and calculation, and lacks unified and automated tools to quickly select static equivalent acceleration or call floor response spectrum (FRS) for spectral analysis.

[0013] (5) The existing technology is not suitable for the particularity of nuclear fusion devices

[0014] Nuclear fusion device structures involve vacuum chambers, superconducting magnets, support systems and other components, and the working environment is extremely harsh, requiring higher accuracy and automation of seismic effect analysis. However, the existing technology mainly targets general buildings or equipment, and lacks effective support for such multi-type and parameter-rich structures.

[0015] Therefore, the present application proposes a nuclear fusion device building seismic spectrum automatic fitting and structure checking system to solve the above technical problems. SUMMARY

[0016] The main purpose of the present application is to provide a nuclear fusion device building seismic spectrum automatic fitting and structure checking system, which can realize automatic retrieval of seismic parameters, automatic generation of seismic influence curve and automatic derivation of frequency acceleration spectrum based on specified data standards and in combination with device geographical location and structure type. The system has built-in static equivalent acceleration automatic determination and FRS (floor response spectrum) spectral analysis calling logic for modal frequency threshold (30Hz) to support seismic spectrum analysis and static strength checking of nuclear fusion devices, thereby solving the technical problems proposed in the background technology.

[0017] The present application solves the above technical problems by using the following technical solutions:

[0018] A nuclear fusion device building seismic spectrum automatic fitting and structure checking system is realized based on the Matlab platform, comprising:

[0019] (1) GB parameter library and index module, which stores specified data domain (including seismic fortification intensity, design seismic grouping, site soil category, structure damping ratio , maximum horizontal seismic influence coefficient , characteristic period , curve tail end , and index keys (including province / city / county (or latitude and longitude mapping), site category, device structure type, etc.), for automatically retrieving and returning specified data domains using index keys by location, site, and structure type conditions Parameter set {intensity, grouping, site category, structure damping ratio , maximum horizontal seismic influence coefficient , characteristic period };

[0020] (2) Spectrum generation and fitting module, built-in damping correction function and piecewise function, for extracting parameters required for seismic design according to device location, site category, and structure type data in the parameter set, automatically truncating negative values and controlling the trend of the end region of the seismic spectrum to remain smooth and consistent, and finally generating the fitted horizontal seismic influence coefficient curve , vertical seismic influence coefficient curve , and corresponding frequency acceleration spectrum , ;

[0021] (3) Static equivalent and FRS determination module, from the fitting data of the spectrum generation and fitting module, the first natural frequency is obtained through modal analysis , if ≥ 30Hz, directly use the static equivalent acceleration of the corresponding frequency; if < 30Hz, call FRS and perform response spectrum analysis, and provide a set of frequency-banded static equivalent parameters for checking;

[0022] (4) Data export and interface module, outputting readable load spectrum files (Excel / CSV) including seismic influence curves ( ), frequency acceleration spectrum ( ), and data transmission through frequency spectrum analysis / static analysis software interface (text / table format).

[0023] Preferably, the specific calculation process of the spectrum generation and fitting module to generate the fitted horizontal seismic influence coefficient curve, vertical seismic influence coefficient curve, and frequency acceleration spectrum includes:

[0024] S1. Correct the maximum horizontal seismic influence coefficient in the parameter set according to the actual structure damping ratio to obtain the actual maximum horizontal seismic influence coefficient , and calculate the control coefficient of the corresponding spectrum shape, with the calculation formula being:

[0025]

[0026]

[0027]

[0028]

[0029] wherein, is the structural damping ratio, is the decay index of the descending segment, is the slope coefficient of the straight descending segment, is the damping correction coefficient, is the maximum horizontal seismic influence coefficient at the table value, is the maximum horizontal seismic influence coefficient after damping correction;

[0030] S2. The horizontal seismic influence coefficient curve is calculated by a set of preset piecewise functions, and the piecewise functions include:

[0031]

[0032] At this time, the non-negative truncation operation is performed simultaneously, and the calculation formula is:

[0033]

[0034] wherein, is the vibration period, is the characteristic period, is the spectral tail length, and the default value is 6;

[0035] S3. The vertical seismic influence coefficient curve is calculated by the engineering default value , and has:

[0036]

[0037] wherein, the engineering default value is consistent with the common rules of nuclear engineering;

[0038] S4. Let the frequency vector be (unit: Hz), is the preset value, and the default value is Hz, for , take , combine the horizontal seismic influence coefficient curve and the vertical seismic influence coefficient curve to calculate and generate the frequency acceleration spectrum, and the calculation formula is:

[0039] ​

[0040]

[0041] wherein, is the frequency, is the weight acceleration, taken as 9.81 m / s 2 , is the horizontal seismic frequency acceleration spectrum, is the vertical seismic frequency acceleration spectrum.

[0042] Preferably, in the static equivalent and FRS determination module, for the first natural frequency , if ≥ 30 Hz, directly use the corresponding frequency static equivalent acceleration; if < 30 Hz, call the fixed reference system FRS and perform response spectrum analysis, while giving a set of frequency banded static equivalent parameters for checking.

[0043] Preferably, the formula for calculating the static equivalent acceleration is:

[0044]

[0045] wherein, is the static equivalent acceleration;

[0046] After the system generates a static working condition file (such as the *CLOAD card of Abaqus, the F file of ANSYS, and the SPCD load card of NASTRAN), directly write the equivalent acceleration acting on the component mass node or reference point.

[0047] Preferably, the specific operation process of the response spectrum analysis includes:

[0048] Through the built-in floor response spectrum database of the system, according to different parts of the nuclear fusion device (such as the central column support, the vacuum chamber, and the superconducting magnet support structure), preset appropriate FRS parameters in the database;

[0049] When analyzing low-frequency components, automatically match the corresponding FRS parameters and call the response spectrum analysis method in the background to generate a set of static equivalent acceleration parameters suitable for the natural frequency of the component.

[0050] Preferably, after generating the horizontal seismic influence coefficient curve, the vertical seismic influence coefficient curve, and the frequency acceleration spectrum, the spectrum generation and fitting module also calls the time history synthesis algorithm (such as the weighted inverse Fourier transform or the iterative adjustment method) for synthesizing the output seismic acceleration time history record consistent with the spectrum shape. The specific calculation process of the called time history synthesis algorithm includes:

[0051] L1. Select the target seismic acceleration spectrum , Define target spectrum, determine the length of time history to be synthesized (typically 20-40s) and sampling step (typically 0.005-0.01s), calculate total sampling points is:

[0052]

[0053] L2. The system uses the inverse Fourier transform method to generate the initial time history, and the calculation formula is:

[0054]

[0055] where,

[0056]

[0057] is the th frequency component, is the initial amplitude set according to the target spectrum envelope, is a random phase angle, and the value range is ;

[0058] L3. Acceleration spectrum iterative fitting of the current time history is used until the preset requirements are met, and the iterative algorithm used in the acceleration spectrum iterative fitting includes:

[0059] L31. Calculate the response spectrum of the current time history and compare it with the target spectrum to obtain the error function , and the calculation formula is: ;

[0060] L32. According to the error distribution, correct the frequency domain amplitude and update the time history again;

[0061] L33. Repeat the L31-L32 operation until the error meets the tolerance requirement: is a set of minimum values, and typically =0.05, i.e. within 5% error.

[0062] L4. Generate horizontal seismic components in the plane direction (X, Y direction) using mutually independent random phases;

[0063] L5. Use the vertical target spectrum to generate the vertical seismic component;

[0064] L6. Finally output three sets of mutually orthogonal seismic components as the ground motion time history (X, Y, Z) for time history analysis.

[0065] Preferably, the spectrum generation and fitting module generates horizontal seismic influence coefficient curves, vertical seismic influence coefficient curves and frequency acceleration spectra, and simultaneously generates corresponding seismic influence curves and acceleration spectra for = 0.02, 0.03, 0.04, 0.05, 0.07, etc. The envelope spectrum is obtained by taking the maximum value of all damping ratio results, and finally forms a conservative spectrum, which is used to support automatic docking with structural analysis software. Through the built-in fixed reference system (FRS) database of the system, the seismic time history record consistent with the target response spectrum is synthesized.

[0066] Preferably, the data export and interface module automatically outputs a complete report containing input parameters, calculation process, result curves, equivalent load and log after the generation of seismic spectrum and the completion of structure checking, so as to ensure the reviewability and traceability of the process.

[0067] Preferably, the data export and interface module is used for converting the spectrum curve in the cycle domain to the acceleration spectrum in the frequency domain, and outputting the standardized data interface file in the Excel / CSV / text format as a readable load spectrum file, which is directly called by structural analysis software such as ABAQUS, ANSYS and NASTRAN.

[0068] In another aspect, the application also discloses a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to make the processor execute the steps of the above method.

[0069] In still another aspect, the application also discloses a computer device, which comprises a memory and a processor, and the memory stores a computer program, and the computer program is executed by the processor to make the processor execute the steps of the above method.

[0070] From the above technical solution, the application provides a nuclear fusion device building seismic spectrum automatic fitting and structure checking system. Compared with the prior art, the application has the following advantages:

[0071] 1. The application can quickly and accurately obtain key parameters such as seismic fortification intensity, design seismic grouping, site category, damping ratio and characteristic period by establishing a GB specification database and an automatic retrieval mechanism, replacing the traditional manual table lookup and manual fitting method, significantly improving the efficiency of the seismic design process, reducing the error risk caused by manual operation, improving the efficiency and reducing human error.

[0072] 2. The present application can automatically generate horizontal and vertical earthquake influence coefficient curves and frequency acceleration spectrum meeting the requirements of the specified building seismic design code standard by setting the damping correction formula and the specification segmentation function in the spectrum generation and fitting module, realize automatic fitting of the earthquake spectrum, and avoid the uncertainty that may be caused by traditional manual fitting.

[0073] 3. The present application can automatically convert the generated earthquake frequency acceleration spectrum into a load input file that can be directly recognized by mainstream structural analysis software such as Abaqus, ANSYS, NASTRAN, etc., saving the steps of manual format conversion and data entry, reducing the deviation introduced by secondary processing, realizing seamless connection with the structural analysis software, and thus improving the automation level of seismic checking analysis.

[0074] 4. The present application can automatically select different seismic loading methods according to the modal analysis results through the built-in 30Hz judgment logic in the system, ensure that components in different frequency ranges can meet the specification requirements, reduce manual judgment, improve the reliability of calculation, realize the intelligentization of modal frequency threshold judgment and loading method.

[0075] 5. The present application can automatically generate reports after completing the analysis, ensure that the entire seismic design process is auditable and traceable, and realize the traceability of the whole process and the automatic generation of reports.

[0076] 6. The present application can realize automatic retrieval of seismic parameters, automatic generation of earthquake influence curves, and automatic export of frequency acceleration spectrum, and support seismic spectrum analysis and static strength checking of nuclear fusion devices.

[0077] It should be understood that the content described in this part is not intended to identify key or important features of the embodiments of the present application, nor is it used to limit the scope of the present application. Other features of the present application will become apparent through the following description. Of course, any product implementing the present application does not necessarily need to achieve all the advantages described above. BRIEF DESCRIPTION OF DRAWINGS

[0078] The drawings accompanying the specification of this application form a part thereof, serve to provide further understanding of the present application, and together with the description of the exemplary embodiments of the present application and their description serve to explain the present application, and do not constitute improper limitations on the present application. In the drawings:

[0079] Fig. 1 The figure is a schematic diagram of the nuclear fusion device building earthquake spectrum automatic fitting and structure checking process of the system of the present application;

[0080] Fig. 2 The figure is a schematic diagram of the earthquake influence coefficient curve provided by the embodiment of the present application;

[0081] Fig. 3 The figure is a schematic diagram of the earthquake frequency acceleration spectrum provided by the embodiment of the present application. DETAILED DESCRIPTION

[0082] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. The embodiments in the present application and the features in the embodiments can be combined with each other without conflict. Based on the embodiments in the present application, all other embodiments obtained by a person of ordinary skill in the art without creative labor fall within the protection scope of the present application.

[0083] In the embodiments, refer to the detailed description Figs. 1 to 3 .

[0084] The embodiments of the present application aim at the problems of low efficiency, lack of automation, non-uniform data interface and lack of targeted support in the process of anti-seismic design and structure checking of existing nuclear fusion devices, and propose a nuclear fusion device building seismic spectrum automatic fitting and structure checking system, which is realized based on a Matlab platform. The specific system modules include:

[0085] (A) GB parameter library and index module, which stores specified data fields (including seismic fortification intensity, design seismic grouping, site soil category, structure damping ratio , maximum horizontal seismic influence coefficient , characteristic period , curve tail , etc.) and index keys (including province / city / county (or longitude and latitude mapping), site category, device structure type, etc.), which are used to automatically retrieve and return the parameter set {intensity, grouping, site category, structure damping ratio , maximum horizontal seismic influence coefficient , characteristic period } composed of specified data fields according to the device location, site category and structure type data in the parameter set using the index keys;

[0086] (B) spectrum generation and fitting module, which internally has a damping correction function and a segmented function, which is used to extract the parameters required for seismic design according to the device location, site category and structure type data in the parameter set, automatically truncate negative values and control the tail of the spectrum to change uniformly, and finally generate the fitted horizontal seismic influence coefficient curve , vertical seismic influence coefficient curve and corresponding frequency acceleration spectrum , ;

[0087] (C) static equivalent and FRS determination module, which obtains the first natural frequency from the fitting data of the spectrum generation and fitting module through modal analysis, and determines whether the first natural frequency ≥ 30Hz, directly using the static force equivalent acceleration of the corresponding frequency; if < 30Hz, calling FRS and performing response spectrum analysis, and giving a set of frequency-banded static force equivalent parameters for checking;

[0088] (D) Data export and interface module, outputting readable load spectrum files (Excel / CSV) including seismic influence curves (F(t)) ), frequency acceleration spectrum (S(f)) ), and transmitting data through spectrum analysis / static force analysis software interface (text / table format), and also being able to realize automatic docking with various analysis software.

[0089] Based on the above system, the present application also proposes embodiment 1, a data processing method performed based on the above system, as shown in Fig. 1 , the processing flow includes:

[0090] STEP1. Selecting position and structure conditions

[0091] The user selects or inputs the construction site (administrative area or latitude and longitude), design earthquake grouping (Group ID), site category (Site Class), and device structure type, etc.

[0092] STEP2. Automatically retrieving GB parameter set

[0093] The system establishes a database of national standards (GB50011, GB50267, etc.), and can automatically extract key parameters such as seismic fortification intensity, design grouping, site category, damping ratio, , , etc. according to the device location, site category, and structure type, realize full-process automation, and avoid manual table lookup and errors.

[0094] The system retrieves the corresponding parameter set from the built-in GB database, including: seismic fortification intensity (unitless), design earthquake grouping (unitless), site soil category (I / II / III / IV), structure damping ratio (unitless, typically 0.02-0.08), maximum horizontal seismic influence coefficient (unitless, based on g ), characteristic period (unit: s), and spectrum tail length (unit: s, default 6).

[0095] At this time, by establishing GB specification database and automatic retrieval mechanism, the traditional manual table lookup and manual fitting method is replaced, which can quickly and accurately obtain key parameters such as seismic fortification intensity, design earthquake grouping, site category, damping ratio and characteristic period, significantly improve the efficiency of seismic design process, and reduce the error risk caused by manual operation.

[0096] STEP3. Damping correction and spectrum amplitude calculation

[0097] The maximum value of the horizontal seismic influence coefficient in the parameter set According to the actual structure damping ratio The actual maximum horizontal seismic influence coefficient is obtained by correction And calculate the spectral shape control coefficient 、 、 :

[0098]

[0099]

[0100]

[0101]

[0102] where, is the structure damping ratio, is the attenuation index of the descending segment, is the slope coefficient of the straight descending segment, is the damping correction coefficient, is the tabulated maximum horizontal seismic influence coefficient, is the maximum horizontal seismic influence coefficient after damping correction.

[0103] STEP4. Horizontal seismic influence coefficient curve

[0104] Let the period (s), define the piecewise function:

[0105]

[0106] And make non-negative truncation:

[0107]

[0108] where, is the vibration period, is the characteristic period, is the spectrum tail length, which is 6 by default.

[0109] The system is based on Matlab programming to realize the automatic generation of standard spectrum segmentation function (rising section, platform section, curved descending section, straight descending section), support damping correction, non-negative truncation and vertical spectrum correction, to ensure the strict consistency of the spectrum curve with GB standard.

[0110] At this time, the system has built-in damping correction formula and standard segmentation function in Matlab platform, which can automatically generate horizontal and vertical earthquake influence coefficient curve and frequency acceleration spectrum that meet the requirements of "Code for Seismic Design of Buildings" (GB50011) and "Standard for Seismic Design of Nuclear Power Plants" (GB50267), to ensure the strict consistency of the results with national standard and avoid the uncertainty that may be caused by traditional manual fitting.

[0111] STEP5. Vertical earthquake influence coefficient curve

[0112] The default value of the project is 2 / 3 (this value is consistent with the commonly used rules for nuclear projects):

[0113]

[0114]

[0115] STEP6. Frequency acceleration spectrum generation

[0116] Let the frequency vector be (unit: Hz, default Hz). For , take:

[0117]

[0118] Combined with , , we get:

[0119]

[0120]

[0121] where is the frequency, is the weight acceleration, which is 9.81 m / s 2 , is the horizontal earthquake frequency acceleration spectrum, is the vertical earthquake frequency acceleration spectrum.

[0122] At this time, the system automatically generates the frequency acceleration spectrum in the range of 0-100 Hz based on the period domain spectrum , And output Excel / CSV standard data interface that can be directly called by structure analysis software, to realize seamless connection between spectrum generation and structure analysis.

[0123] STEP7. Data export

[0124] Export table file (Excel / CSV) containing Frequency_Hz, A_horizontal_mps2, A_vertical_mps2, which can be directly read by spectrum analysis or static analysis software.

[0125] STEP8. Static equivalent acceleration automatic determination

[0126] (1) If ≥ 30Hz:

[0127] Static equivalent acceleration :

[0128]

[0129] The system automatically inputs the as static equivalent load into static structure analysis.

[0130] (2) If < 30Hz:

[0131] The system will automatically perform the following operations:

[0132] (2a) Select the appropriate FRS (map by location / level / device area)

[0133] (2b) Perform response spectrum analysis (component / connector oriented);

[0134] (2c) Output a set of static equivalent acceleration parameters (can be weighted by frequency band or modal participation coefficient to form an envelope) according to natural frequency and FRS characteristics for structure checking. This automatic switching logic reduces manual judgment and ensures compliance with nuclear engineering specifications.

[0135] STEP9. Report and traceability

[0136] Generate automatic report of spectrum graph ( , ), key parameters, determination process and load list to ensure the auditability of the checking process.

[0137] At this time, the system can also automatically generate a report, including input parameters, damping correction calculation process, spectrum curve, acceleration spectrum, modal determination results, load file path and calculation log, to ensure that the seismic design process has review and traceability functions.

[0138] Based on the method mentioned in system embodiment 1, in specific embodiments, the following data processing scheme is further adopted:

[0139] Embodiment 2: Seismic spectrum automatic fitting and acceleration spectrum generation method based on Matlab

[0140] Based on the above method, input the specified conditions and sequentially execute database retrieval, damping correction calculation and seismic influence coefficient curve generation, and as shown in Fig. 2 and Fig. 3 , the system outputs the following results:

[0141] (1) Graph (horizontal and vertical same figure);

[0142] (2) Acceleration spectrum graph (horizontal and vertical same figure);

[0143] (3) Excel / CSV file (column name Frequency_Hz, A_horizontal_mps2, A_vertical_mps2).

[0144] At this time, the system automatically generates a report, including:

[0145] (a) Input parameter table (intensity, grouping, site category, , , , etc.);

[0146] (b) Damping correction calculation process;

[0147] (c) Graph and frequency spectrum;

[0148] (d) Export file path;

[0149] (e) System version number and specification version number.

[0150] Embodiment 3: Automatic docking with structural analysis software and equivalent load generation method

[0151] Based on the generation of seismic acceleration spectrum, this embodiment can also automatically dock mainstream structural analysis software (such as Abaqus, ANSYS, NASTRAN), and automatically convert it into the input file of the structural analysis software (such as Abaqus, ANSYS, NASTRAN). After the system generates the frequency acceleration spectrum, it can be automatically converted into the load input format required by the target software, and combined with the first natural frequency to automatically determine the loading mode, reduce the manual processing link, and improve the efficiency and accuracy of structural checking.

[0152] The specific implementation steps include:

[0153] (1) Modal analysis results import

[0154] (1a) The system reads the modal results file exported by the structure analysis software, such as OP2 (NASTRAN), ODB (Abaqus), and RST (ANSYS).

[0155] (1b) Automatically extract the first natural frequency (Hz) and store it in the database.

[0156] (2) Using the same static equivalent acceleration automatic judgment rule mentioned in example 2 above, if ≥ 30Hz, the system determines that the component is a high-frequency component, and uses the static equivalent acceleration method, if < 30Hz, the system determines that the component is a low-frequency connection component, and uses the response spectrum analysis method, and calls the FRS response spectrum for frequency spectrum input, and then imports the results into the structure software according to the component classification, reduces the manual data conversion steps and error probability, wherein:

[0157] (2a) The system directly takes the frequency corresponding to the horizontal acceleration value according to the acceleration spectrum generated in example 2:

[0158]

[0159] The system generates a static working condition file (such as Abaqus *CLOAD card, ANSYS F file, NASTRAN SPCD load card), directly writes the equivalent acceleration, and acts on the component mass node or reference point;

[0160] (2b) The system generates a complete response spectrum input file according to the generated acceleration spectrum , Combined with the characteristics of FRS, wherein the file format supports mainstream solvers, such as:

[0161] Abaqus *RESPONSESPECTRUM statement;

[0162] ANSYS SPECTRUM command;

[0163] NASTRAN RLOAD and DLOAD card.

[0164] At this time, the system can segment the spectrum values in the frequency range or generate an envelope spectrum as input to ensure calculation safety.

[0165] Here, through the system built-in 30Hz judgment logic, different seismic loading methods can be automatically selected according to the modal analysis results, ensuring that components in different frequency ranges can meet the specification requirements, reducing manual judgment, and improving the reliability of the calculation.

[0166] (3) Automatically generate interface files: the system generates standardized input files according to the user's selected software environment:

[0167] Abaqus:.inp file, containing CLOAD or RESPONSESPECTRUM paragraph;

[0168] ANSYS: Command stream file (.cdb / .txt), containing F, SPECTRUM instructions;

[0169] NASTRAN:.bdf file, containing SPCD, RLOAD, DLOAD cards.

[0170] (4) The system automatically generates a log file when exporting the file, including: modal frequency determination results, corresponding calculation methods (static equivalent or response spectrum), output file path and name, key data of load value or spectrum value.

[0171] At this time, the generated earthquake frequency acceleration spectrum is automatically converted into a load input file that can be directly recognized by mainstream structural analysis software such as Abaqus, ANSYS, and NASTRAN, which can save the steps of manual format conversion and data entry, reduce the deviation introduced by secondary processing, and improve the automation level of seismic checking analysis.

[0172] Example 4: Low-frequency component response analysis method based on floor response spectrum (FRS)

[0173] For some connecting components in nuclear fusion devices with natural frequencies below 30Hz, the existing specification requires that their seismic action should be analyzed by frequency spectrum analysis using floor response spectrum (FRS).

[0174] The purpose is to achieve rapid response analysis and equivalent load generation for low-frequency components through the system's built-in FRS database and automatic calling mechanism, avoiding manual selection and secondary calculation.

[0175] Therefore, the embodiment can also have a built-in floor response spectrum (FRS) database. The system presets appropriate FRS parameters in the database according to different parts of the nuclear fusion device (such as the central column support, the vacuum chamber, the superconducting magnet support structure, etc.). When the user performs low-frequency component analysis, the system can automatically match the corresponding FRS and call the response spectrum analysis method in the background to generate a set of static force equivalent acceleration parameters suitable for the component natural frequency. This function not only improves the analysis efficiency, but also ensures that the processing method conforms to the specifications and the actual situation of the device. The specific implementation steps include:

[0176] (1) Component frequency identification

[0177] (1a) The system automatically reads the modal analysis results and filters out the first natural frequency <30Hz components.

[0178] (1b) Mark these components as "low-frequency components" and enter the FRS analysis process.

[0179] (2) FRS database calling

[0180] The system has a built-in floor response spectrum (FRS) database that can automatically match the low-frequency response requirements of different device parts, which includes:

[0181] (2a) Typical device areas (such as central column support, vacuum chamber, magnet support structure, hoisting beam, etc.);

[0182] (2b) Reference seismic spectrum or acceleration parameters corresponding to each area (defined according to GB50011, GB50267, etc.);

[0183] (2c) User-defined FRS that can be extended by engineering projects.

[0184] When the component position is identified, the system automatically matches the corresponding FRS entry.

[0185] (3) Response spectrum loading mode selection

[0186] The system automatically selects the appropriate loading mode according to the natural frequency of the component and its relationship with the FRS spectrum band:

[0187] If fall within the main frequency band of the FRS spectrum: directly use the FRS spectrum value as input;

[0188] If located outside the spectrum: the system performs interpolation or extrapolation to obtain the modified spectrum value.

[0189] (4) Structure analysis file generation

[0190] The system will automatically generate the calculated response spectrum, and automatically import the load file into the structure analysis software:

[0191] Abaqus: load in *CLOAD or *DLOAD format;

[0192] ANSYS: input in F command or SPECTRUM keyword;

[0193] NASTRAN: realized by RLOAD / DLOAD card.

[0194] (5) Result output and report

[0195] The system automatically generates a "FRS analysis report", which includes: component name and its natural frequency , the matched FRS entry and corresponding parameters, the seismic response spectrum file for structure analysis, the automatically generated structure software input file path, the whole process of determination and calculation log.

[0196] Example 5: Method for synthesizing earthquake time history consistent with response spectrum

[0197] In some complex conditions, only relying on the seismic influence coefficient curve or frequency acceleration spectrum for response analysis is not enough to meet the engineering needs. For example, nuclear fusion devices need to perform nonlinear time history analysis to evaluate the dynamic response of the structure under earthquake action during hoisting, transportation or extreme conditions. Therefore, this embodiment provides a method for synthesizing earthquake time history consistent with the response spectrum based on generating the target acceleration spectrum, which can synthesize earthquake time history consistent with the response spectrum. After generating the seismic acceleration spectrum, the system can call the time history synthesis algorithm (such as weighted inverse Fourier transform or iterative adjustment method) to synthesize the earthquake acceleration time history record consistent with the spectrum shape, so that the time history analysis and spectrum analysis remain consistent and traceable. This function is suitable for complex conditions that require time history analysis, and users can directly import the generated time history curve into the structure dynamic analysis, avoiding relying on external uncertain sources of earthquake records, thereby improving the reliability of nuclear fusion device dynamic response analysis.

[0198] The specific implementation steps include:

[0199] (1) Target spectrum definition

[0200] (1a) Select the target seismic acceleration spectrum generated by Example 2 、 as the target spectrum;

[0201] (1b) Determine the time history length to be synthesized (generally 20-40s), and the sampling step (generally 0.005-0.01s).

[0202] (1c) Calculate the total number of sampling points :

[0203]

[0204] (2) Initial time history generation

[0205] The system uses the inverse Fourier transform method to generate the initial time history:

[0206]

[0207]

[0208] where,

[0209] is the first frequency component, is the initial amplitude set according to the target spectrum envelope, is a random phase angle, with a value range .

[0210] (3) Iterative adjustment process

[0211] Use the acceleration spectrum fitting iterative algorithm:

[0212] (3a) Calculate the response spectrum of the current time history ;

[0213] (3b) Compare with the target spectrum , get the error function:

[0214]

[0215] (3c) According to the error distribution, correct the frequency domain amplitude , update the time history;

[0216] (3d) Repeat the iteration until the error meets the tolerance (generally take =0.05, that is, within 5% error).

[0217] (4) Synthesis of horizontal and vertical components

[0218] For horizontal seismic components (X, Y direction), independent random phase generation method can be used to ensure the independence between directions;

[0219] For vertical seismic components, generate according to the vertical target spectrum ;

[0220] ​Finally, three orthogonal seismic time histories (X, Y, Z) are output, which meet the input requirements of time history analysis.

[0221] (5) Result export

[0222] The system output results include: time history curve (X, Y, Z three directions), fitting before and after the spectrum comparison chart (target spectrum vs synthetic spectrum), time history data file (support TXT / CSV format, including time column and acceleration column), report file, explain the number of iterations, convergence error, final spectrum consistency test results.

[0223] In summary, this method can also synthesize seismic time history records consistent with the target response spectrum, supporting nonlinear time history analysis.

[0224] Example 6: Seismic spectrum envelope generation method under multiple damping ratio conditions

[0225] In the seismic design and structural checking of nuclear fusion devices, the structural damping ratio (D) ) is often not a unique value, and different components and different working conditions may take different damping ratios. The common damping ratio range is 0.02-0.08. Existing methods generally only generate spectrum curves for a single damping ratio, which cannot fully cover the possibility of multiple working conditions. This embodiment provides a method for automatically generating seismic spectra under multiple damping ratio conditions and forming a conservative spectrum through envelope processing for structural strength and stability checking.

[0226] Therefore, this embodiment can also batch generate seismic spectra under different damping ratio conditions. The system can simultaneously generate corresponding seismic influence curves and acceleration spectra for =0.02, 0.03, 0.04, 0.05, 0.07, etc. in one calculation, and give "envelope spectrum" in the output results. The envelope spectrum takes the maximum value of all damping ratio results as a reference input for conservative design, providing multiple working condition coverage for engineering checking, and improving the safety margin and integrity of the design

[0227] The specific implementation steps include:

[0228] (1) Input damping ratio set

[0229] The user can manually input or the system can automatically call the common damping ratio set:

[0230]

[0231] The system will generate seismic influence coefficient curves and acceleration spectra for each respectively.

[0232] (2) Damping correction calculation

[0233] For each , the maximum influence coefficient is corrected according to the following formula:

[0234]

[0235]

[0236]

[0237]

[0238] wherein, is the structural damping ratio, is the decay index of the descending section, is the slope coefficient of the straight descending section, is the damping correction coefficient, is the maximum influence coefficient of the table value horizontal earthquake, is the maximum influence coefficient of the horizontal earthquake after damping correction.

[0239] (3) Spectrum curve generation

[0240] For each damping ratio , the horizontal / vertical earthquake influence coefficient curve , is generated.

[0241] Calculate the frequency acceleration spectrum:

[0242]

[0243]

[0244] (4) Envelope spectrum calculation

[0245] The system compares all spectra point by point, and takes the maximum value as the envelope spectrum. The envelope spectrum ensures that the calculation results are not less than the envelope value under any damping ratio condition, meeting the conservative design requirements.

[0246] (5) Result export

[0247] System output: spectrum curve under each damping ratio, comparison chart of envelope spectrum and single damping ratio spectrum, and envelope spectrum data file (Frequency_Hz, A_horizontal_env_mps2, A_vertical_env_mps2).

[0248] As mentioned above, this method can also generate spectra under different damping ratios and output envelope spectra in batches, meeting the conservative design requirements.

[0249] In yet another aspect, the present application also discloses a computer readable storage medium storing a computer program, which, when executed by a processor, causes the processor to perform the steps of the above method.

[0250] In yet another aspect, the present application also discloses a computer device comprising a memory and a processor, wherein the memory stores a computer program, and the computer program, when executed by the processor, causes the processor to perform the steps of the above method.

[0251] In yet another embodiment provided in the present application, a computer program product containing instructions, which, when run on a computer, causes the computer to execute the system for automatic fitting of seismic spectrum of fusion device construction and structure checking in any of the above embodiments.

[0252] It can be understood that the system provided in the embodiments of the present application corresponds to the method provided in the embodiments of the present application, and the related content explanation, examples and beneficial effects can refer to the corresponding part in the above method.

[0253] The present application also provides an electronic device comprising a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory complete mutual communication through the communication bus,

[0254] The memory is used for storing a computer program.

[0255] The processor is used for executing the program stored on the memory, and realizes the above.

[0256] The communication bus mentioned in the above electronic device can be a peripheral component interconnect bus or an extended industry standard architecture bus, etc. The communication bus can be divided into an address bus, a data bus, a control bus, etc.

[0257] The communication interface is used for communication between the above electronic device and other devices.

[0258] The memory can include a random access memory, and can also include a non-volatile memory, for example, at least one disk memory. Optionally, the memory can also be at least one storage device located away from the aforementioned processor.

[0259] The processor mentioned above can be a general-purpose processor, including a central processing unit, a network processing unit, etc.; can also be a digital signal processor, an application specific integrated circuit, a field programmable gate array or other programmable logic device, a discrete gate or transistor logic device, a discrete hardware component.

[0260] It should be further explained that the electronic device also includes a terminal device, which can also be referred to as a terminal, a user equipment, a mobile station, a mobile terminal, etc. The terminal device can be a mobile phone, a smart television, a wearable device, a tablet computer, a computer with wireless transceiver function, a virtual reality terminal device, an augmented reality terminal device, a wireless terminal in industrial control, a wireless terminal in unmanned driving, a wireless terminal in remote surgery, a wireless terminal in smart grid, a wireless terminal in transportation safety, a wireless terminal in smart city, a wireless terminal in smart home, etc. The embodiments of the present application do not limit the specific technology and specific device form of the terminal device.

[0261] In the above embodiments, all or part can be realized by software, hardware, firmware or any combination thereof. When realized by software, all or part can be realized in the form of a computer program product. The computer program product includes one or more computer instructions. When loaded and executed on a computer, the computer program instructions produce the processes or functions described in the embodiments of the present application.

[0262] The above only describes the preferred embodiments of the present application and is not intended to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

[0263] In addition, it should be noted that if the embodiments of the present application involve directional indications (such as up, down, left, right, front, back, etc.), the directional indications are only used to explain the relative position relationship, motion condition, etc. between components in a certain posture, and if the certain posture changes, the directional indications also change accordingly.

[0264] In addition, if the embodiments of the present application involve descriptions of “first”, “second”, etc., the descriptions of “first”, “second”, etc. are only for description purposes and cannot be understood as indicating or implying the relative importance of the indicated technical features or implicitly indicating the number of the indicated technical features. Therefore, the features limited by “first” and “second” can explicitly or implicitly include at least one of the features. In addition, the meaning of “and / or” appearing throughout the text includes three parallel solutions. Taking “A and / or B” as an example, it includes A solution, or B solution, or A and B solutions. In addition, in the embodiments of the present application, “multiple” means more than two. In addition, the technical solutions of each embodiment can be combined with each other, but it must be based on the realization of a person of ordinary skill in the art, and when the combination of technical solutions contradicts each other or cannot be realized, it should be considered that the combination of technical solutions does not exist and is not within the protection scope of the present application.

Claims

1. A nuclear fusion device building seismic spectrum automatic fitting and structure checking system based on a Matlab platform, characterized in that, Comprise: The GB parameter library and index module stores the specified data field and index key, which is used to automatically retrieve and return the parameter set composed of the specified data field using the index key according to the specified conditions; The spectrum generation and fitting module is used to call the GB database to extract the parameters required for seismic design according to the device location, site category and structure type data inside the parameter set, automatically truncate negative values and control the trend of the end region of the seismic spectrum to remain smooth and consistent, and finally generate the fitted horizontal seismic influence coefficient curve, vertical seismic influence coefficient curve and corresponding frequency acceleration spectrum; The static equivalent and FRS determination module obtains the first natural frequency from the fitting data of the spectrum generation and fitting module through modal analysis and determines according to the preset threshold, respectively corresponding to the static equivalent acceleration and response spectrum analysis operation; The data export and interface module outputs the readable load spectrum file including the seismic influence curve and frequency acceleration spectrum, and transmits data through the spectrum analysis / static analysis software interface; The specific calculation process of the spectrum generation and fitting module to generate the fitted horizontal seismic influence coefficient curve, vertical seismic influence coefficient curve and corresponding frequency acceleration spectrum includes: S1. The maximum horizontal seismic influence coefficient in the parameter set is modified according to the actual structural damping ratio According to the actual structural damping ratio The maximum horizontal seismic influence coefficient is modified to obtain the actual maximum horizontal seismic influence coefficient And the control coefficient corresponding to the spectrum shape is calculated, and the calculation formula is: Wherein, is the structural damping ratio, is the decay index of the descending section, is the slope coefficient of the straight descending section, is the damping correction coefficient, is the maximum influence coefficient of the table value horizontal earthquake, is the maximum influence coefficient of the horizontal earthquake after damping correction; S2. Calculate the horizontal seismic influence coefficient curve by a set of pre-set piecewise functions ; S3. Calculate the vertical seismic influence coefficient curve by engineering default values There are: The engineering default value Consistent with the rules commonly used in nuclear engineering; S4. Set the frequency vector , , take , take , combine the horizontal seismic influence coefficient curve and the vertical seismic influence coefficient curve , calculate the generated frequency acceleration spectrum, and the calculation formula is: wherein is the frequency, is the weight acceleration, is the horizontal seismic frequency acceleration spectrum, is the vertical seismic frequency acceleration spectrum.

2. The nuclear fusion device construction seismic spectrum automatic fitting and structure checking system of claim 1, wherein, The piecewise function includes: At this time, the non-negative truncation operation is performed simultaneously, and the calculation formula is: wherein, is a vibration period, is a characteristic period, is a spectral tail length.

3. The nuclear fusion device construction seismic spectrum automatic fitting and structure checking system of claim 1, wherein, In the static equivalent and FRS judging module, for the first natural frequency , if ≥ 30 Hz, directly use the static equivalent acceleration corresponding to the frequency; if < 30 Hz, call the fixed reference system FRS and perform response spectrum analysis, and give a set of frequency banded static equivalent parameters for checking.

4. The nuclear fusion device construction seismic spectrum automatic fitting and structure checking system of claim 3, wherein, The calculation formula of the static equivalent acceleration is: wherein is the static equivalent acceleration; After the system generates the static working condition file, the equivalent acceleration is directly written into the file and acts on the component mass node or reference point.

5. The nuclear fusion device construction seismic spectrum automatic fitting and structure checking system of claim 4, wherein, The specific operation process of the response spectrum analysis includes: Through the built-in floor response spectrum database in the system, according to the different parts of the nuclear fusion device, the appropriate FRS parameters are preset in the database; When analyzing low-frequency components, automatically match the corresponding FRS parameters and call the response spectrum analysis method in the background to generate a set of static equivalent acceleration parameters suitable for the natural frequency of the component.

6. The nuclear fusion device construction seismic spectrum automatic fitting and structure checking system of claim 1, wherein, After the spectrum generation and fitting module generates the horizontal seismic influence coefficient curve, vertical seismic influence coefficient curve and frequency acceleration spectrum, it also calls the time history synthesis algorithm to output the seismic acceleration time history record consistent with the spectrum shape, and the specific calculation process of the calling time history synthesis algorithm includes: L1. Select the target seismic acceleration spectrum , Define the target spectrum, determine the length of time history needed to be synthesized and the sampling step , calculate the total number of sampling points as: L2. The system uses the inverse Fourier transform method to generate the initial time history, and the calculation formula is: wherein is the th frequency component, is the initial amplitude set according to the target spectral envelope, is a random phase angle in the range ; L3. Acceleration spectrum iterative fitting is used to fit the current time history until the preset requirements are met; L4. Independent random phases are used to generate horizontal seismic components in the plane direction; L5. Use vertical target spectrum Fitting generates vertical seismic component; L6. Finally, output three sets of mutually orthogonal seismic components as the ground motion time history for time history analysis.

7. The nuclear fusion device construction seismic spectrum automatic fitting and structure checking system of claim 6, wherein, The iterative algorithm used in the L3 step to perform acceleration spectrum iterative fitting includes: L31. Calculate the response spectrum for the current time history and compare to the target spectrum to obtain an error function with the formula: ; L32. Correct the frequency domain amplitude according to the error distribution and re-update the schedule. L33. Repeat the iteration of the L31-L32 operations until the error meets the tolerance requirement: is a set of minima.

8. The nuclear fusion device construction seismic spectrum automatic fitting and structure checking system of claim 1, wherein, After the spectrum generation and fitting module generates the horizontal seismic influence coefficient curve, vertical seismic influence coefficient curve and frequency acceleration spectrum, it also simultaneously generates the corresponding seismic influence curve and acceleration spectrum for the specified damping ratio, and gives the corresponding envelope spectrum in the output result, which takes the maximum value of all damping ratio results, and finally forms a conservative spectrum, which is used to support automatic docking with the structure analysis software. Through the built-in fixed reference system database in the system, the seismic time history record consistent with the target response spectrum is synthesized.

9. The nuclear fusion device construction seismic spectrum automatic fitting and structure checking system of claim 1, wherein, The data export and interface module is used for converting a period domain spectrum curve into a frequency domain acceleration spectrum, and outputting a standardized data interface file in Excel / CSV / text format as a readable load spectrum file for directly calling structural analysis software such as Abaqus, ANSYS, NASTRAN and the like.

Images