Automatic building seismic spectrum fitting and structure checking system for nuclear fusion device

Abstract

The invention discloses an automatic building seismic spectrum fitting and structure checking system for a nuclear fusion device, and the system comprises a GB parameter library and index module which is used for carrying out the retrieval through an index key, and returning a parameter set composed of a specified data field; the spectrum generation and fitting module is used for calling a GB database to extract parameters required by an anti-seismic design according to the internal data of the parameter set, automatically cutting off negative values, controlling the change of a spectrum tail to be consistent, and generating a fitted horizontal earthquake influence coefficient curve, a vertical earthquake influence coefficient curve and a corresponding frequency acceleration spectrum; the static equivalent and FRS judgment module is used for obtaining a first inherent frequency through modal analysis in fitting data, judging according to a preset threshold value, and respectively and correspondingly adopting static equivalent acceleration and response spectrum analysis operation; and the data export and interface module is used for outputting the readable load spectrum file. According to the invention, automatic generation and export operation of an earthquake influence curve and a frequency acceleration spectrum can be realized, and anti-seismic spectrum analysis and static strength check of a nuclear fusion device are supported.

Description

Technical Field

[0001] This invention relates to the field of seismic resistance technology for nuclear fusion device buildings, specifically to a Matlab-based automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings. Background Technology

[0002] With the increasing global demand for clean energy, nuclear fusion has emerged as a potential high-efficiency and clean energy source. A nuclear fusion device, as a novel clean energy device, typically comprises a vacuum chamber, a superconducting magnet, a central support system, and surrounding buildings. These structures must withstand complex static and dynamic loads throughout their service life, with seismic forces being a crucial factor that must be considered. To ensure the safety and reliability of the device and its critical components under seismic conditions, seismic response spectrum fitting and structural strength verification must be performed in accordance with relevant national standards.

[0003] In current engineering practice, seismic action analysis and structural verification for nuclear fusion devices primarily rely on manual lookup and fitting methods based on GB standards. Engineers first manually consult appendices and influence coefficient curves in the *Code for Seismic Design of Buildings* and relevant nuclear power plant standards, based on the device's geographical location, structural type, and design seismic group. Then, they adjust these curves using damping ratios. Afterward, they manually or with the aid of general software perform curve fitting on this data and plot the seismic influence curve. These standards provide seismic influence coefficient curves and calculation methods under different seismic fortification intensities, site categories, and structural damping ratios.

[0004] However, the following problems still exist in actual use: (1) The data interface is not uniform, which increases the secondary conversion work and causes problems in the automation of seismic parameter acquisition. Existing methods require manual review of GB standards and manual input of parameters, which is inefficient and prone to errors. This invention aims to automate the extraction of parameters such as seismic fortification intensity, design seismic group, site category, structural damping ratio, maximum value of horizontal seismic influence coefficient, and characteristic period by establishing a GB database interface.

[0005] (2) Manual table lookup and curve fitting are inefficient and redundant, leading to longer project cycles and reduced efficiency. Engineers often need to manually consult tables and curves in the appendices of GB standards based on the location and design grouping of the device, and then combine them with correction formulas such as damping ratio to draw seismic impact curves. This process is tedious and prone to human error, making it difficult to meet the requirements of efficiency and accuracy for the verification of complex structures in nuclear fusion devices.

[0006] (3) Seismic spectrum generation and structural analysis are disconnected, making it difficult to achieve full-process automation and error control. In the existing process, there is no direct interface between the generated seismic impact curves and structural analysis software. Data typically requires manual secondary processing to obtain seismic acceleration spectra suitable for spectral analysis. This approach is not only inefficient but may also introduce data bias.

[0007] (4) Static equivalent analysis lacks automated support For different components whose first natural frequency is higher or lower than 30Hz in the modal analysis results, the standard requires the use of different seismic action calculation methods. However, the existing process relies heavily on manual judgment and calculation, and lacks a unified and automated tool to quickly select the static equivalent acceleration or call the floor response spectrum (FRS) for spectral analysis.

[0008] (5) Existing technologies are not sufficiently adaptable to the special characteristics of nuclear fusion devices. Nuclear fusion devices involve various components such as vacuum chambers, superconducting magnets, and support systems. Their operating environment is extreme, demanding higher accuracy and automation in seismic analysis. However, existing technologies are primarily geared towards general buildings or equipment, lacking effective support for such diverse structures with varied parameters.

[0009] To address these technical problems, this application proposes an automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings. Summary of the Invention

[0010] The main objective of this invention is to provide an automatic seismic spectrum fitting and structural verification system for nuclear fusion devices. Based on specified data standards and combined with the device's geographical location and structural type, the system can automatically retrieve seismic parameters, automatically generate seismic influence curves, and automatically export frequency acceleration spectra. Furthermore, the system has built-in logic for automatically determining static equivalent acceleration and calling FRS (floor response spectrum) spectrum analysis for modal frequency thresholds (30Hz), thereby supporting the seismic spectrum analysis and static strength verification of nuclear fusion devices and solving the technical problems mentioned in the background art.

[0011] The present invention solves the above-mentioned technical problems by adopting the following technical solutions: An automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings, implemented on the Matlab platform, includes: (1) GB parameter library and index module, which stores specified data fields (including seismic fortification intensity, design seismic group, site soil type, structural damping ratio). Maximum value of horizontal earthquake influence coefficient Characteristic Period , curve tail end (etc.) and index keys (including province / city / county (or latitude and longitude mapping), site category, device structure type, etc.), used to automatically retrieve and return a parameter set consisting of a specified data domain using the index keys based on location, site, and structure type conditions {intensity, grouping, site category, structural damping ratio, etc.}. Maximum value of horizontal earthquake influence coefficient Characteristic Period }; (2) The spectrum generation and fitting module has built-in damping correction function and piecewise function. It is used to extract the parameters required for seismic design from the GB database based on the device location, site category and structural type data in the parameter set, automatically truncate negative values ​​and control the change trend of the end area of ​​the seismic spectrum to keep it stable and consistent, and finally generate the fitted horizontal seismic influence coefficient curve. Vertical seismic influence coefficient curve and the corresponding frequency acceleration spectrum , ; (3) Static equivalence and FRS determination module: The first natural frequency is obtained from the fitting data of the spectrum generation and fitting module through modal analysis. ,like For frequencies ≥30Hz, the corresponding static equivalent acceleration is directly used; if If the frequency is less than 30Hz, the FRS is invoked and the response spectrum analysis is performed. At the same time, a set of frequency banded static equivalent parameters are provided for verification. (4) Data export and interface module, outputting earthquake impact curves ( ), frequency acceleration spectrum ( The system can read load spectrum files (Excel / CSV) and transmit data through a spectrum analysis / static analysis software interface (text / table format).

[0012] Preferably, the specific calculation process for 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: S1. Maximize the horizontal seismic influence coefficient in the parameter set. According to the actual structural damping ratio After correction, the maximum value of the actual horizontal earthquake influence coefficient is obtained. And calculate the control coefficients for the corresponding spectral shape, using the following formula:

[0013]

[0014]

[0015]

[0016] in, The structural damping ratio, The decay index of the descent segment. The slope coefficient of the descending segment. This is the damping correction factor. This represents the maximum influence coefficient of a horizontal earthquake, as shown in the table. This represents the maximum influence coefficient of a horizontal earthquake after damping correction. S2. The horizontal seismic influence coefficient curve is calculated using a set of preset piecewise functions. The piecewise function includes:

[0017] At this time, a non-negative truncation operation is performed simultaneously, and the calculation formula is as follows:

[0018] in, For the oscillation period, For characteristic period, This is the duration of the spectral ending, with a default value of 6; S3. Calculate the vertical seismic influence coefficient curve using the project default values. ,have:

[0019] Among them, the default value of the project Consistent with the rules commonly used in nuclear-grade engineering; S4. Let the frequency vector be... (Unit: Hz) This is the default value. Hz, for ,Pick Combined with the horizontal earthquake influence coefficient curve Vertical Seismic Influence Coefficient Curve The generated frequency acceleration spectrum is calculated using the following formula:

[0020]

[0021] in, For frequency, Let the acceleration be weight, taken as 9.81 m / s². 2 , The horizontal earthquake frequency acceleration spectrum, This represents the vertical seismic frequency acceleration spectrum.

[0022] Preferably, in the static equivalent and FRS determination module, for the first natural frequency... ,like For frequencies ≥30Hz, the corresponding static equivalent acceleration is directly used; if If the frequency is less than 30Hz, the fixed reference system (FRS) is invoked and response spectrum analysis is performed. At the same time, a set of banded static equivalent parameters are provided for verification.

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

[0024] in, It is a static equivalent acceleration; After the system generates a static load case file (such as Abaqus's *CLOAD card, ANSYS's F file, or NASTRAN's SPCD load card), the equivalent acceleration is directly written and applied to the component mass node or reference point.

[0025] Preferably, the specific operation procedure for the response spectrum analysis includes: By using the system's built-in floor response spectrum database, appropriate FRS parameters are preset in the database according to different parts of the nuclear fusion device (such as the central column support, vacuum chamber, superconducting magnet support structure, etc.). When performing low-frequency component analysis, the corresponding FRS parameters are automatically matched and the response spectrum analysis method is called in the background to generate a set of static equivalent acceleration parameters that are adapted to the natural frequency of the component.

[0026] 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 a time history synthesis algorithm (such as weighted inverse Fourier transform or iterative adjustment method) to synthesize a seismic acceleration time history record with a spectral shape. The specific calculation process of calling the time history synthesis algorithm includes: L1. Select from the target seismic acceleration spectrum , Define the target spectrum and determine the required time history length for synthesis. (typically 20–40 seconds) and sampling step size (Generally taken as 0.005–0.01s), calculate the total number of sampling points. for:

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

[0028] in,

[0029] For the first One frequency component, The initial amplitude is set according to the target spectral envelope. The phase angle is random, and its value range is [value range missing]. ; L3. Iteratively fit the current time history using the acceleration spectrum until the preset requirements are met. The iterative algorithms used for acceleration spectrum iterative fitting at this time include: L31. Calculate the response spectrum at the current time history. and the target spectrum By comparison, the error function is obtained. The calculation formula is as follows: ; L32. Correct the frequency domain amplitude according to the error distribution. And update the schedule again; L33. Repeat the L31-L32 operations iteratively until the error meets the tolerance requirements: , which is a set of local minima, and is generally taken as =0.05, meaning the error is within 5%.

[0030] L4. Horizontal seismic components in the plane (X and Y directions) are generated using mutually independent random phases; L5. Using vertical target spectrum Vertical seismic components are generated through fitting. L6. Finally, three sets of mutually orthogonal seismic components are output as the seismic motion time histories (X, Y, Z) for time history analysis.

[0031] Preferably, the spectrum generation and fitting module, after generating the horizontal seismic influence coefficient curve, the vertical seismic influence coefficient curve, and the frequency acceleration spectrum, also simultaneously targets... Common damping ratios such as 0.02, 0.03, 0.04, 0.05, and 0.07 are used to generate corresponding seismic influence curves and acceleration spectra, and the corresponding envelope spectrum is given in the output results. The envelope spectrum is the maximum value of all damping ratio results, and finally a conservative spectrum is formed to support automatic docking with structural analysis software. Through the system's built-in Fixed Reference System (FRS) database, seismic time history records consistent with the target response spectrum are synthesized.

[0032] Preferably, after the seismic spectrum generation and structural verification are completed, the data export and interface module also automatically outputs a complete report containing input parameters, calculation process, result curves, equivalent loads and logs, ensuring that the process is auditable and traceable.

[0033] Preferably, the data export and interface module is used to convert the periodic domain spectrum curve into a frequency domain acceleration spectrum and output standardized data interface files in Excel / CSV / text formats as readable load spectrum files for direct use by structural analysis software such as Abaqus, ANSYS, and NASTRAN.

[0034] In another aspect, the present invention 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 method described above.

[0035] In another aspect, the present invention also discloses a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor performs the steps of the method described above.

[0036] As can be seen from the above technical solution, the present invention provides an automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings. Compared with the prior art, the present invention has the following advantages: 1. This invention, by establishing a GB standard database and an automatic retrieval mechanism, replaces the traditional methods of manual table lookup and manual fitting. It can quickly and accurately obtain key parameters such as seismic fortification intensity, design seismic group, site category, damping ratio, and characteristic period, significantly improving the efficiency of the seismic design process and reducing the risk of errors caused by manual operation, thereby increasing efficiency and reducing human error.

[0037] 2. By setting damping correction formulas and standard piecewise functions in the spectrum generation and fitting module, this invention can automatically generate horizontal and vertical seismic influence coefficient curves and frequency acceleration spectra that meet the requirements of specified building seismic design codes, thereby achieving automatic fitting of seismic spectra and avoiding the uncertainties that may be caused by traditional manual fitting.

[0038] 3. This invention can automatically convert the generated seismic frequency acceleration spectrum into a load input file that can be directly recognized by mainstream structural analysis software such as Abaqus, ANSYS, and NASTRAN, eliminating the need for manual format conversion and data entry, reducing deviations introduced by secondary processing, and achieving seamless integration with structural analysis software, thereby improving the automation level of seismic verification analysis.

[0039] 4. This invention, through its built-in 30Hz judgment logic, can automatically select different seismic loading methods based on modal analysis results, ensuring that components in different frequency ranges can meet the specifications, reducing manual judgment, improving the reliability of calculations, and realizing intelligent determination of modal frequency thresholds and loading methods.

[0040] 5. This invention can automatically generate a report after the analysis is completed, ensuring that the entire seismic design process is auditable and traceable, and realizing full-process traceability and automatic report generation.

[0041] 6. This invention can realize automatic retrieval of earthquake parameters, automatic generation of earthquake influence curves, and automatic export of frequency acceleration spectrum, and supports seismic spectrum analysis and static strength verification of nuclear fusion devices.

[0042] It should be understood that the descriptions in this section are not intended to identify key or essential features of embodiments of the invention, nor are they intended to limit the scope of the invention. Other features of the invention will become readily apparent from the following description. Of course, implementing any product of the invention does not necessarily require achieving all of the advantages described above simultaneously. Attached Figure Description

[0043] The accompanying drawings, which form part of this application, are used to provide a further understanding of the invention. The illustrative embodiments of the invention and their descriptions are used to explain the invention and do not constitute an undue limitation of the invention. In the drawings: Figure 1 This is a schematic diagram of the automatic seismic spectrum fitting and structural verification process for the nuclear fusion device building of the present invention. Figure 2 This is a schematic diagram of the earthquake influence coefficient curve provided in an embodiment of the present invention; Figure 3 This is a schematic diagram of the seismic frequency acceleration spectrum provided in an embodiment of the present invention. Detailed Implementation

[0044] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of them. Unless otherwise specified, the embodiments and features in the embodiments of this application can be combined with each other. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0045] For details in the embodiments, please refer to Figures 1 to 3 .

[0046] This invention addresses the problems of low efficiency, insufficient automation, inconsistent data interfaces, and lack of targeted support in the seismic design and structural verification of existing nuclear fusion devices. It proposes an automatic seismic spectrum fitting and structural verification system for nuclear fusion devices, implemented on the Matlab platform. The specific system modules include: (A) GB parameter library and index module, storing specified data fields (including seismic fortification intensity, design seismic group, site soil type, structural damping ratio). Maximum value of horizontal earthquake influence coefficient Characteristic Period , curve tail end (etc.) and index keys (including province / city / county (or latitude and longitude mapping), site category, device structure type, etc.), used to automatically retrieve and return a parameter set consisting of a specified data domain using the index keys based on location, site, and structure type conditions {intensity, grouping, site category, structural damping ratio, etc.}. Maximum value of horizontal earthquake influence coefficient Characteristic Period }; (B) The spectrum generation and fitting module has built-in damping correction and piecewise functions. It is used to extract the parameters required for seismic design from the GB database based on the device location, site category and structural type data in the parameter set, automatically truncate negative values ​​and control the consistency of spectrum tail changes, and finally generate the fitted horizontal seismic influence coefficient curve. Vertical seismic influence coefficient curve and the corresponding frequency acceleration spectrum , ; (C) Static equivalence and FRS determination module: The first natural frequency is obtained from the fitted data of the spectrum generation and fitting module through modal analysis. ,like For frequencies ≥30Hz, the corresponding static equivalent acceleration is directly used; if If the frequency is less than 30Hz, the FRS is invoked and the response spectrum analysis is performed. At the same time, a set of frequency banded static equivalent parameters are provided for verification. (D) Data export and interface module, outputting earthquake impact curves ( ), frequency acceleration spectrum ( It can read load spectrum files (Excel / CSV) including those for spectral analysis and static analysis software, and transmit data through the interface (text / table format). It can also automatically interface with various analysis software.

[0047] Based on the above system, this invention also proposes Embodiment 1 in the specific data processing process, which is a data processing method performed based on the above system, such as... Figure 1 As shown, the processing flow includes: STEP 1. Selecting Location and Structural Conditions Users can select or input conditions such as the construction location (administrative region or latitude and longitude), design seismic group (GroupID), site class (SiteClass), and device structure type on the interface.

[0048] STEP2. Automatically retrieve GB parameter sets The system establishes a database of national standards (GB50011, GB50267, etc.) and can automatically extract seismic fortification intensity, design group, site category, damping ratio, etc., based on the device location, site category, and structural type. , Key parameters are controlled to automate the entire process, avoiding manual table lookups and errors.

[0049] The system retrieves the corresponding parameter set from the built-in GB database, including: seismic fortification intensity (unitless), design seismic group (unitless), site soil type (I / II / III / IV), and structural damping ratio. (Unitless, typically 0.02–0.08), maximum value of horizontal seismic influence coefficient (No unit, in) (Based on) characteristic period (Unit: s) and duration of the spectral tail (Unit: s, default is 6).

[0050] At this point, by establishing a GB standard database and an automatic retrieval mechanism, the traditional methods of manual table lookup and manual fitting are replaced. This allows for the rapid and accurate acquisition of key parameters such as seismic fortification intensity, design seismic group, site category, damping ratio, and characteristic period, significantly improving the efficiency of the seismic design process and reducing the risk of errors caused by manual operation.

[0051] STEP 3. Damping Correction and Spectral Amplitude Calculation The maximum value of the horizontal seismic influence coefficient in the parameter set. According to the actual structural damping ratio The maximum value of the actual horizontal earthquake influence coefficient was obtained after correction. And calculate the spectral control coefficients. , , :

[0052]

[0053]

[0054]

[0055] in, The structural damping ratio, The decay index of the descent segment. The slope coefficient of the descending segment. This is the damping correction factor. This represents the maximum influence coefficient of a horizontal earthquake, as shown in the table. This represents the maximum influence coefficient of a horizontal earthquake after damping correction.

[0056] STEP 4. Horizontal Seismic Influence Coefficient Curve

[0057] Set period (s), define a piecewise function:

[0058] And perform non-negative truncation:

[0059] in, For the oscillation period, For characteristic period, This is the duration of the spectral ending, which defaults to 6.

[0060] This system uses Matlab programming to automatically generate standard spectral segmentation functions (ascending segment, plateau segment, descent segment, and linear descent segment), supports damping correction, non-negative truncation, and vertical spectral correction, and ensures that the spectral curve is strictly consistent with the GB standard.

[0061] At this point, the system has built-in damping correction formulas and standard piecewise functions in the Matlab platform, which can automatically generate horizontal and vertical seismic influence coefficient curves and frequency acceleration spectra that meet the requirements of the "Code for Seismic Design of Buildings" (GB50011) and the "Standard for Seismic Design of Nuclear Power Plants" (GB50267), ensuring that the results are strictly consistent with national standards and avoiding the uncertainties that may be caused by traditional manual fitting.

[0062] STEP 5. Vertical Seismic Influence Coefficient Curve

[0063] The default value for this project is 2 / 3 (this value is consistent with the common rules for nuclear-grade projects):

[0064]

[0065] STEP 6. Frequency Acceleration Spectrum Generation

[0066] Let the frequency vector (Unit: Hz, default) Hz). For ,Pick:

[0067] Combination , ,get:

[0068]

[0069] in, For frequency, Let the acceleration be weight, taken as 9.81 m / s². 2 , The horizontal earthquake frequency acceleration spectrum, This represents the vertical seismic frequency acceleration spectrum.

[0070] At this point, the system automatically generates a frequency acceleration spectrum in the range of 0–100 Hz based on the periodic domain spectrum. , It outputs standardized data interfaces that can be directly called by structural analysis software such as Excel / CSV, enabling seamless integration of spectrum generation and structural analysis.

[0071] STEP 7. Data Export Export a tabular file (Excel / CSV) containing Frequency_Hz, A_horizontal_mps2, and A_vertical_mps2, which can be directly read by spectrum analysis or static analysis software.

[0072] STEP 8. Automatic Determination of Static Equivalent Acceleration (1) If ≥30Hz: Static equivalent acceleration :

[0073] The system will automatically As a static equivalent load input for static structural analysis.

[0074] (2) If <30Hz: The system will automatically perform the following operations: (2a) Select the appropriate FRS (by location / hierarchy / device area mapping) (2b) Perform response spectrum analysis (for components / connectors); (2c) Based on the natural frequency and FRS characteristics, output a set of static equivalent acceleration parameters (which can be weighted by frequency band or main mode participation coefficient to form an envelope) for structural verification. This automatic switching logic reduces manual judgment and ensures compliance with nuclear engineering specifications.

[0075] STEP 9. Reporting and Traceability Generate spectrum ( , Automated reporting of key parameters, judgment process, and load list ensures the auditability of the verification process.

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

[0077] Based on the method mentioned in System Embodiment 1 above, the following data processing scheme is further adopted in a specific embodiment: Example 2: Automatic Seismic Spectrum Fitting and Acceleration Spectrum Generation Method Based on Matlab Based on the above method, after inputting specified conditions and sequentially executing database retrieval, damping correction calculation, and seismic influence coefficient curve generation, as follows: Figure 2 and Figure 3 As shown, the system outputs the following results: (1) Curve graph (horizontal and vertical graphs are the same); (2) Acceleration spectrum (horizontal and vertical are the same); (3) Excel / CSV file (column names are Frequency_Hz, A_horizontal_mps2, A_vertical_mps2).

[0078] The system will then automatically generate a report, which includes: (a) Input parameter table (intensity, grouping, site category, , , wait); (b) Damping correction calculation process; (c) Curve graphs and frequency spectrum; (d) Export file path; (e) System version number and specification version number.

[0079] Example 3: Automatic Integration with Structural Analysis Software and Equivalent Load Generation Method Based on the generated seismic acceleration spectrum, this embodiment can also automatically interface with mainstream structural analysis software (such as Abaqus, ANSYS, and NASTRAN), automatically converting it into input files for the structural analysis software (such as Abaqus, ANSYS, and NASTRAN). After generating the frequency acceleration spectrum, the system can automatically convert it into the load input format required by the target software, and combine it with the first natural frequency. Automatically determine the loading method, reduce manual processing steps, and improve the efficiency and accuracy of structural verification.

[0080] Its specific implementation steps include: (1) Importing modal analysis results (1a) The system reads the modal results files exported by the structural analysis software, such as OP2 (NASTRAN), ODB (Abaqus), and RST (ANSYS).

[0081] (1b) Automatically extract the first natural frequency (Hz), and stored in the database.

[0082] (2) Using the same static equivalent acceleration automatic determination rule mentioned in Example 2 above, if If the frequency is ≥30Hz, the system determines that the component is a high-frequency component and uses the static equivalent acceleration method. If the frequency is <30Hz, the system determines that the component is a low-frequency connection component, uses response spectrum analysis, and calls the FRS response spectrum for spectrum input. The results are then imported into the structural software according to component classification, reducing manual data conversion steps and the probability of errors. Specifically: (2a) The system directly takes the frequency based on the acceleration spectrum generated in Example 2. Corresponding horizontal acceleration values:

[0083] The system generates static load case files (such as Abaqus' *CLOAD card, ANSYS' F file, and NASTRAN's SPCD load card), and directly writes the equivalent acceleration, which is applied to the component mass nodes or reference points. (2b) The system is based on the generated acceleration spectrum , Combined with FRS characteristics, a complete reaction spectrum input file is generated, with a file format that supports mainstream solvers, such as: The *RESPONSESPECTRUM statement in Abaqus; The SPECTRUM command in ANSYS; RLOAD and DLOAD cards in NASTRAN.

[0084] At this point, the system can write the spectral values ​​within the frequency range in segments, or generate an envelope spectrum as input, to ensure computational security.

[0085] Here, the system's built-in 30Hz judgment logic can automatically select different seismic loading methods based on modal analysis results, ensuring that components in different frequency ranges can meet the specifications, reducing manual judgment and improving the reliability of calculations.

[0086] (3) Automatically generate interface files: The system generates standardized input files based on the user's selected software environment. Abaqus: .inp files containing either CLOAD or RESPONSESPECTRUM paragraphs; ANSYS: Command stream file (.cdb / .txt) containing F and SPECTRUM commands; NASTRAN: .bdf file, containing SPCD, RLOAD, and DLOAD cards.

[0087] (4) When exporting files, the system automatically generates a log file, which includes: modal frequencies. The determination result, the corresponding calculation method (static equivalent or response spectrum), the output file path and name, and the key data of load value or spectrum value.

[0088] At this point, the generated seismic 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. This eliminates the need for manual format conversion and data entry, reduces deviations introduced by secondary processing, and improves the automation level of seismic verification analysis.

[0089] Example 4: Low-frequency component response analysis method based on floor response spectrum (FRS) For some connecting components in nuclear fusion devices with natural frequencies below 30Hz, existing standards require that their seismic action be analyzed using floor response spectrum (FRS).

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

[0091] Therefore, this embodiment can also incorporate a floor response spectrum (FRS) database. The system presets appropriate FRS parameters in the database based on different parts of the nuclear fusion device (such as the central column support, vacuum chamber, 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 equivalent acceleration parameters adapted to the component's natural frequency. This function not only improves analysis efficiency but also ensures that the processing method conforms to specifications and the actual situation of the device. Its specific implementation steps include: (1) Component frequency identification (1a) The system automatically reads the modal analysis results and selects the first natural frequency. Components with a frequency of <30Hz.

[0092] (1b) Mark these components as “low-frequency components” and enter them into the FRS analysis process.

[0093] (2) FRS database call The system has a built-in floor response spectrum (FRS) database, which can automatically match the low-frequency response requirements of different equipment locations, including: (2a) Typical device areas (such as central column support, vacuum chamber, magnet support structure, hoisting beam, etc.); (2b) Reference seismic spectrum or acceleration parameters for each region (defined according to standards such as GB50011 and GB50267); (2c) Scalable user-defined FRS for engineering projects.

[0094] Once the component location is identified, the system automatically matches the corresponding FRS entry.

[0095] (3) Selection of response spectrum loading method The system is based on the natural frequency of the components. Based on its relationship with FRS bands, the system automatically selects the appropriate loading method. like If the frequency falls within the main frequency band of the FRS spectrum: directly use the FRS spectrum value as input; like Located outside the spectrum: The system performs interpolation or extrapolation to obtain the corrected spectral value.

[0096] (4) Generation of structural analysis files The system will automatically generate a load file from the calculated response spectrum that can be directly imported into structural analysis software. Abaqus: Loads in *CLOAD or *DLOAD format; ANSYS: Enter using the F command or the SPECTRUM keyword; NASTRAN: Implemented using RLOAD / DLOAD cards.

[0097] (5) Results Output and Reporting The system automatically generates an "FRS Analysis Report," which includes: component name and its inherent frequency. The matched FRS entries and corresponding parameters, the seismic response spectrum file used for structural analysis, the path to the automatically generated structural software input file, and the complete log of the judgment and calculation process.

[0098] Example 5: Method for synthesizing seismic time histories consistent with response spectra In certain complex operating conditions, relying solely on seismic influence coefficient curves or frequency acceleration spectra for response analysis is insufficient to meet engineering requirements. For example, during hoisting, transportation, or extreme conditions, nuclear fusion devices require nonlinear time history analysis to assess the dynamic response of the structure under seismic loading. Therefore, this embodiment provides a method for synthesizing seismic time history records consistent with the response spectrum based on the generated target acceleration spectrum. After generating the seismic acceleration spectrum, the system can invoke time history synthesis algorithms (such as weighted inverse Fourier transform or iterative adjustment methods) to synthesize seismic acceleration time history records consistent with the spectrum shape, ensuring consistency and traceability between time history analysis and spectral analysis. This function is suitable for complex operating conditions requiring time history analysis. Users can directly import the generated time history curves into structural dynamic analysis, avoiding reliance on seismic records from uncertain external sources, thereby improving the reliability of dynamic response analysis of nuclear fusion devices.

[0099] Its specific implementation steps include: (1) Definition of target spectrum (1a) Select the target seismic acceleration spectrum generated by Example 2 , As the target spectrum; (1b) Determine the required time history length for synthesis (typically 20–40 seconds), and sampling step size (Typically, the value is 0.005–0.01 s). (1c) Calculate the total number of sampling points :

[0100] (2) Initial time history generation The system uses the inverse Fourier transform method to generate the initial time history:

[0101]

[0102] in, For the first One frequency component, The initial amplitude is set according to the target spectral envelope. The phase angle is random, and its range is [range missing]. .

[0103] (3) Iterative adjustment process An iterative algorithm for fitting the acceleration spectrum is used: (3a) Calculate the response spectrum at the current time history ; (3b) and the target spectrum By comparison, the error function is obtained:

[0104] (3c) Correct the frequency domain amplitude according to the error distribution. Update the schedule; (3d) Repeat the iterations until the error meets the tolerance. (generally taken) =0.05, meaning the error is within 5%.

[0105] (4) Combination of horizontal and vertical components For the horizontal seismic components (X and Y directions), an independent random phase generation method can be used to ensure the independence between directions; For the vertical seismic component, according to the vertical target spectrum Generate by fitting; The final output consists of three mutually orthogonal ground motion time histories (X, Y, Z), which meet the input requirements for time history analysis.

[0106] (5) Results Export The system output includes: time history curves (X, Y, Z directions), a comparison of the spectra before and after fitting (target spectrum vs. synthesized spectrum), time history data files (supporting TXT / CSV format, including time and acceleration columns), and a report file explaining the number of iterations, convergence error, and final spectrum consistency test results.

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

[0108] Example 6: Seismic Spectrum Envelope Generation Method under Multiple Damping Ratio Conditions In the seismic design and structural verification of nuclear fusion devices, the structural damping ratio ( The damping ratio is often not a unique value; different components and different operating conditions may have different damping ratios. Common damping ratio ranges from 0.02 to 0.08. Existing methods generally only generate spectral curves for a single damping ratio, which cannot fully cover the possibility of multiple operating conditions. This embodiment provides a method for automatically generating seismic spectra under multiple damping ratio conditions and forming conservative spectra through envelope processing for structural strength and stability verification.

[0109] Therefore, this embodiment can also generate seismic spectra under different damping ratios in batches. The system can simultaneously perform calculations for... The system generates corresponding seismic influence curves and acceleration spectra for common damping ratios such as 0.02, 0.03, 0.04, 0.05, and 0.07, and provides an "envelope spectrum" in the output. The envelope spectrum is the maximum value of all damping ratio results, serving as a reference input for conservative design. This provides multi-condition coverage for engineering verification, improving the safety margin and completeness of the design. Its specific implementation steps include: (1) Set of input damping ratios Users can manually input or the system can automatically retrieve a set of common damping ratios:

[0110] The system will process each Seismic influence coefficient curves and acceleration spectra are generated respectively.

[0111] (2) Damping correction calculation For each The maximum influence coefficient is adjusted according to the following formula:

[0112]

[0113]

[0114]

[0115] in, The structural damping ratio, The decay index of the descent segment. The slope coefficient of the descending segment. This is the damping correction factor. This represents the maximum influence coefficient of a horizontal earthquake, as shown in the table. This represents the maximum influence coefficient of a horizontal earthquake after damping correction.

[0116] (3) Generation of spectral curves For each damping ratio Generate horizontal / vertical seismic influence coefficient curves. , .

[0117] Calculate the frequency acceleration spectrum:

[0118]

[0119] (4) Envelope spectrum calculation The system is for all The spectrum is compared point by point, and the maximum value is taken as the envelope spectrum. The envelope spectrum ensures that the calculation result is not less than the envelope value under any damping ratio, thus meeting the conservative design requirements.

[0120] (5) Results Export System output: spectral curves for each damping ratio, a comparison chart of the envelope spectrum and the spectrum for a single damping ratio, and envelope spectrum data files (Frequency_Hz, A_horizontal_env_mps2, A_vertical_env_mps2).

[0121] In summary, this method can also generate spectra under different damping ratios in batches and output envelope spectra, meeting the requirements of conservative design.

[0122] In another aspect, the present invention 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 method described above.

[0123] In another aspect, the present invention also discloses a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor performs the steps of the method described above.

[0124] In another embodiment provided in this application, a computer program product containing instructions is also provided, which, when run on a computer, causes the computer to execute any of the above embodiments of the automatic fitting and structural verification system for seismic spectrum of nuclear fusion device building.

[0125] It is understood that the system provided in the embodiments of the present invention corresponds to the method provided in the embodiments of the present invention, and the explanation, examples and beneficial effects of the relevant content can be referred to the corresponding parts of the above methods.

[0126] This application also provides an electronic device, including a processor, a communication interface, a memory, and a communication bus, wherein the processor, communication interface, and memory communicate with each other via the communication bus. Memory, used to store computer programs; The processor, when executing a program stored in memory, performs the above.

[0127] The communication bus mentioned in the above-mentioned electronic devices can be a standard bus for interconnecting peripheral components or an extended industrial standard structure bus, etc. This communication bus can be divided into address bus, data bus, control bus, etc.

[0128] The communication interface is used for communication between the aforementioned electronic devices and other devices.

[0129] The memory may include random access memory or non-volatile memory, such as at least one disk storage device. Optionally, the memory may also be at least one storage device located remotely from the aforementioned processor.

[0130] The processors mentioned above can be general-purpose processors, including central processing units, network processors, etc.; they can also be digital signal processors, application-specific integrated circuits, field-programmable gate arrays or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0131] It should also be noted that electronic devices include terminal devices, which can also be called terminals, user equipment, mobile stations, mobile terminals, etc. Terminal devices can be mobile phones, smart TVs, wearable devices, tablets, computers with wireless transceiver capabilities, virtual reality terminal devices, augmented reality terminal devices, wireless terminals in industrial control, wireless terminals in autonomous driving, wireless terminals in remote surgery, wireless terminals in smart grids, wireless terminals in transportation safety, wireless terminals in smart cities, wireless terminals in smart homes, and so on. The embodiments of this application do not limit the specific technologies or device forms used in the terminal devices.

[0132] In the above embodiments, implementation can be achieved, in whole or in part, through software, hardware, firmware, or any combination thereof. When implemented using software, it can be implemented, in whole or in part, as a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, all or part of the processes or functions described in the embodiments of this application are generated.

[0133] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

[0134] Furthermore, it should be noted that if any directional indication (such as up, down, left, right, front, back, etc.) is involved in the embodiments of the present invention, the directional indication is only used to explain the relative positional relationship and movement of each component in a specific posture. If the specific posture changes, the directional indication will also change accordingly.

[0135] Furthermore, if the embodiments of this invention involve descriptions such as "first" or "second," these descriptions are for descriptive purposes only and should not be construed as indicating or implying their relative importance or implicitly specifying the number of technical features indicated. Therefore, a feature defined with "first" or "second" may explicitly or implicitly include at least one of those features. Additionally, the meaning of "and / or" throughout the text includes three parallel solutions; for example, "A and / or B" includes solution A, solution B, or a solution where both A and B are satisfied simultaneously. Furthermore, in the embodiments of this invention, "multiple" refers to two or more. Moreover, the technical solutions of the various embodiments can be combined with each other, but this must be based on the ability of those skilled in the art to implement them. When the combination of technical solutions is contradictory or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection claimed by this invention.

Claims

1. An automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings, implemented on the Matlab platform, characterized in that, include: The GB parameter library and index module stores specified data fields and index keys, which are used to automatically retrieve and return parameter sets composed of specified data fields using the index keys according to specified conditions. The spectrum generation and fitting module is used to extract the parameters required for seismic design from the GB database based on the device location, site category and structure type data in the parameter set, automatically truncate negative values ​​and control the change trend of the end area of ​​the seismic spectrum to keep it stable 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 fitted data of the spectrum generation and fitting module through modal analysis and determines it according to the preset threshold, and respectively adopts static equivalent acceleration and response spectrum analysis operations. The data export and interface module outputs readable load spectrum files, including seismic influence curves and frequency acceleration spectra, and transmits the data through a spectrum analysis / static analysis software interface.

2. The automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings as described in claim 1, characterized in that, The specific calculation process for 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. Maximize the horizontal seismic influence coefficient in the parameter set. According to the actual structural damping ratio After correction, the maximum value of the actual horizontal earthquake influence coefficient is obtained. And calculate the control coefficients for the corresponding spectral shape, using the following formula: in, The structural damping ratio, The decay index of the descent segment. The slope coefficient of the descending segment. This is the damping correction factor. This represents the maximum influence coefficient of a horizontal earthquake, as shown in the table. This represents the maximum influence coefficient of a horizontal earthquake after damping correction. S2. The horizontal seismic influence coefficient curve is calculated using a set of preset piecewise functions. ; S3. Calculate the vertical seismic influence coefficient curve using the project default values. ,have: Among them, the default value of the project Consistent with the rules commonly used in nuclear-grade engineering; S4. Let the frequency vector be... , As a preset value, for ,Pick Combined with the horizontal earthquake influence coefficient curve Vertical Seismic Influence Coefficient Curve The generated frequency acceleration spectrum is calculated using the following formula: in, For frequency, For weight acceleration, The horizontal earthquake frequency acceleration spectrum, This represents the vertical seismic frequency acceleration spectrum.

3. The automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings as described in claim 2, characterized in that, The piecewise function includes: At this time, a non-negative truncation operation is performed simultaneously, and the calculation formula is as follows: in, For the oscillation period, For characteristic period, This is the duration of the end of the melody.

4. The automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings as described in claim 2, characterized in that, In the static equivalent and FRS determination module, for the first natural frequency ,like For frequencies ≥30Hz, the corresponding static equivalent acceleration is directly used; if If the frequency is less than 30Hz, the fixed reference system (FRS) is invoked and response spectrum analysis is performed. At the same time, a set of banded static equivalent parameters are provided for verification.

5. The automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings as described in claim 4, characterized in that, The formula for calculating the static equivalent acceleration is as follows: in, It is a static equivalent acceleration; After the system generates the static condition file, the equivalent acceleration is directly written and applied to the component mass node or reference point.

6. The automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings as described in claim 5, characterized in that, The specific operational procedures for the response spectrum analysis include: By using the system's built-in floor response spectrum database, appropriate FRS parameters are preset in the database according to different parts of the nuclear fusion device; When performing low-frequency component analysis, the corresponding FRS parameters are automatically matched and the response spectrum analysis method is called in the background to generate a set of static equivalent acceleration parameters that are adapted to the natural frequency of the component.

7. The automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings as described in claim 2, characterized in that, 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 a time history synthesis algorithm to output seismic acceleration time history records consistent with the spectrum shape. The specific calculation process of calling the time history synthesis algorithm includes: L1. Select from the target seismic acceleration spectrum , Define the target spectrum and determine the required time history length for synthesis. and sampling step size Calculate the total number of sampling points for: L2. The system uses the inverse Fourier transform method to generate the initial time history, and the calculation formula is as follows: in, For the first One frequency component, The initial amplitude is set according to the target spectral envelope. The phase angle is random, and its value range is [value range missing]. ; L3. Iteratively fit the current time history using the acceleration spectrum until the preset requirements are met; L4. Horizontal seismic components in the plane direction are generated using mutually independent random phases; L5. Using vertical target spectrum Vertical seismic components are generated through fitting. L6. Finally, three sets of mutually orthogonal seismic components are output as the seismic motion time histories for time history analysis.

8. The automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings as described in claim 7, characterized in that, The iterative algorithm used in step L3 to perform the acceleration spectrum iterative fitting includes: L31. Calculate the response spectrum at the current time history. and the target spectrum By comparison, the error function is obtained. The calculation formula is as follows: ; L32. Correct the frequency domain amplitude according to the error distribution. And update the schedule again; L33. Repeat the L31-L32 operations iteratively until the error meets the tolerance requirements: , which is a set of minimum values.

9. The automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings as described in claim 1, characterized in that, The spectrum generation and fitting module generates horizontal seismic influence coefficient curves, vertical seismic influence coefficient curves, and frequency acceleration spectra. It also generates corresponding seismic influence curves and acceleration spectra for specified damping ratios and provides corresponding envelope spectra in the output results. The envelope spectrum takes the maximum value of all damping ratio results to form a conservative spectrum, which is used to support automatic docking with structural analysis software. Through the system's built-in fixed reference system database, it synthesizes seismic time history records consistent with the target response spectrum.

10. The automatic seismic spectrum fitting and structural verification system for nuclear fusion device buildings as described in claim 1, characterized in that, The data export and interface module is used to convert periodic domain spectral curves into frequency domain acceleration spectra and output standardized data interface files in Excel / CSV / text formats as readable load spectrum files for direct use by structural analysis software such as Abaqus, ANSYS, and NASTRAN.

Images