WO2017077609A1 - 構造解析方法、及び構造解析プログラム - Google Patents

構造解析方法、及び構造解析プログラム Download PDF

Info

Publication number
WO2017077609A1
WO2017077609A1 PCT/JP2015/081098 JP2015081098W WO2017077609A1 WO 2017077609 A1 WO2017077609 A1 WO 2017077609A1 JP 2015081098 W JP2015081098 W JP 2015081098W WO 2017077609 A1 WO2017077609 A1 WO 2017077609A1
Authority
WO
WIPO (PCT)
Prior art keywords
model
structural analysis
analysis
strength
boundary condition
Prior art date
Application number
PCT/JP2015/081098
Other languages
English (en)
French (fr)
Inventor
孝典 根岸
康弘 川島
和也 山浦
正雄 福島
玉木　剛
俊雄 荒井
敏靖 大原
Original Assignee
富士通株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 富士通株式会社 filed Critical 富士通株式会社
Priority to JP2017548568A priority Critical patent/JP6645508B2/ja
Priority to PCT/JP2015/081098 priority patent/WO2017077609A1/ja
Priority to TW105135268A priority patent/TW201721462A/zh
Publication of WO2017077609A1 publication Critical patent/WO2017077609A1/ja
Priority to US15/960,744 priority patent/US11003816B2/en

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/292—Extracting wanted echo-signals
    • G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
    • G01S7/2927—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods by deriving and controlling a threshold value
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/32—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device
    • H04N1/32358—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device using picture signal storage, e.g. at transmitter
    • H04N1/32443—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device using picture signal storage, e.g. at transmitter with asynchronous operation of the image input and output devices connected to the memory
    • H04N1/32448—Controlling data flow to or from the memory in relation to the available memory capacity
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/32—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device
    • H04N1/32358—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device using picture signal storage, e.g. at transmitter
    • H04N1/32443—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device using picture signal storage, e.g. at transmitter with asynchronous operation of the image input and output devices connected to the memory
    • H04N1/32454—Controlling data flow to or from the memory in relation to the amount of data, e.g. file size
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/32—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device
    • H04N1/333—Mode signalling or mode changing; Handshaking therefor
    • H04N1/33376—Mode signalling or mode changing; Handshaking therefor according to characteristics or state of one of the communicating parties, e.g. available memory capacity

Definitions

  • the present invention relates to a structure analysis method and a structure analysis program.
  • a technique for performing structural analysis using a finite element method is known.
  • structural analysis using the finite element method for example, displacement and stress of a structure are obtained by solving simultaneous linear equations.
  • a direct method or an iterative method can be given as a large classification.
  • the direct method is a solution method that directly obtains a solution by triangulating a coefficient matrix
  • the iterative method is a solution method that obtains an approximate solution by iterative calculation.
  • Non-Patent Document 1 introduces an algorithm of a direct method suitable for a sparse matrix.
  • Non-Patent Document 2 introduces an iterative algorithm suitable for parallel processing and a calculation example thereof.
  • an object of the present invention is to provide a structure analysis method and a structure analysis program that can shorten analysis time.
  • model information is acquired, based on the acquired model information, the model size of the model to be analyzed is evaluated, and based on the evaluation result of the model size, the finite element of the model
  • a structure analysis method and a structure analysis program for selecting an algorithm for solving simultaneous linear equations of a structure analysis solver using a method as a direct method or an iterative method are proposed.
  • the boundary condition of the model to be analyzed is acquired, the strength of the boundary condition is evaluated based on the acquired boundary condition, and the strength of the boundary condition is evaluated. Based on the results, a structure analysis method and a structure analysis program for selecting whether the algorithm of the simultaneous linear equations of the structure analysis solver using the finite element method of the model is a direct method or an iterative method are proposed.
  • model information and boundary conditions of a model to be analyzed are acquired, the model scale of the model is evaluated based on the acquired model information, and the acquired boundary conditions Based on the above, the strength of the boundary condition is evaluated. Based on the evaluation result of the model scale and the strength evaluation result of the boundary condition, the simultaneous linear equation solution of the structural analysis solver using the finite element method of the model is performed.
  • a structure analysis method and a structure analysis program for selecting whether the algorithm is a direct method or an iterative method are proposed.
  • analysis time can be shortened.
  • FIG. 1 is an explanatory diagram showing an example of the operation of the structural analysis apparatus according to the present invention.
  • FIG. 2 is an explanatory view showing another operation example by the structural analysis apparatus according to the present invention.
  • FIG. 3 is an explanatory diagram showing elements and nodes used in the structural analysis of the finite element method.
  • FIG. 4 is a block diagram illustrating a hardware configuration example of the structural analysis apparatus.
  • FIG. 5 is a block diagram illustrating a functional configuration example of the structural analysis apparatus.
  • FIG. 6 is an explanatory diagram illustrating an example of model information.
  • FIG. 7 is an explanatory diagram illustrating an example of a memory capacity threshold value according to the number of elements.
  • FIG. 1 is an explanatory diagram showing an example of the operation of the structural analysis apparatus according to the present invention.
  • FIG. 2 is an explanatory view showing another operation example by the structural analysis apparatus according to the present invention.
  • FIG. 3 is an explanatory diagram showing elements and nodes used in the structural analysis of the finite element
  • FIG. 8 is a flowchart illustrating an example of a selection processing procedure based on the model scale by the structural analysis apparatus.
  • FIG. 9 is an explanatory diagram illustrating a model example used for deriving a threshold value.
  • FIG. 10 is an explanatory diagram showing an example of deriving a threshold value when the number of elements is 5000.
  • FIG. 11 is an explanatory diagram illustrating an example of deriving a threshold value when the number of elements is 10,000.
  • FIG. 12 is an explanatory diagram illustrating an example of a threshold corresponding to the model scale.
  • FIG. 13 is an explanatory diagram illustrating an example of calculating the strength of the constraint condition and an example of selecting a solution.
  • FIG. 9 is an explanatory diagram illustrating a model example used for deriving a threshold value.
  • FIG. 10 is an explanatory diagram showing an example of deriving a threshold value when the number of elements is 5000.
  • FIG. 11 is an explanatory diagram illustrating an example of
  • FIG. 14 is a flowchart illustrating an example of a selection procedure based on the strength of the constraint condition by the structural analysis apparatus.
  • FIG. 15 is an explanatory diagram illustrating a system application example.
  • FIG. 16 is an explanatory diagram illustrating an example of an analysis database.
  • FIG. 17 is an explanatory diagram illustrating an example of selecting a solution based on the performance in the system.
  • FIG. 18 is an explanatory diagram of an example of managing the server state and the memory amount.
  • FIG. 19 is a flowchart (part 1) illustrating an example of a structure analysis processing procedure by the system.
  • FIG. 20 is a flowchart (part 2) of the structural analysis processing procedure example by the system.
  • FIG. 1 is an explanatory diagram showing an example of the operation of the structural analysis apparatus according to the present invention.
  • the structural analysis apparatus 100 is a computer that performs structural analysis using a finite element method for an FEM (Finite Element Method) model to be analyzed.
  • the finite element method is a numerical analysis method for approximately solving a differential equation.
  • displacement is evaluated as a stiffness index
  • stress is evaluated as a strength index, and the like.
  • the FEM model to be analyzed is hereinafter also simply referred to as a model, an analysis target model, and the like.
  • a region where an equation is defined is divided into small regions, and an equation in each small region is approximated by a common interpolation function.
  • the small area is referred to as an element 102.
  • a node 103 is arranged at each vertex.
  • a region where an equation is defined is a structure to be analyzed.
  • the analysis target model 101 is a state in which structural analysis by FEM is possible by approximating a structure to be analyzed, and is placed in a simulation space.
  • the simulation space is a space set to represent on the computer an analysis region including the inside of the physical structure to be analyzed and the space outside the structure.
  • the structure to be analyzed is divided into elements 102 which are fine small areas, and the entire shape is represented as a collection of elements 102.
  • a collection of elements 102 is also referred to as a mesh.
  • the size of one lattice in the mesh is also referred to as mesh size.
  • the shape of one element 102 is a tetrahedron, a triangle, etc. in the case of two dimensions, and a tetrahedron, a pentahedron, a hexahedron, etc. in the case of three dimensions.
  • a tetrahedral secondary element is taken as an example of the element 102, but an example of the shape of the element 102 and the relationship between the element 102 and the node 103 will be described in detail with reference to FIG. 2.
  • the developer creates, for example, three-dimensional model information of a three-dimensional solid model in which a structure is represented by polygons using CAD (Computer Aided Design).
  • the three-dimensional model information includes polygon coordinate data, for example.
  • CAE Computer Aided Engineering
  • the CAE has, for example, three pieces of software called a preprocessor, a solver, and a post processor.
  • the CAE preprocessor can create model information 111 representing an FEM model by performing mesh division and boundary condition definition.
  • the CAE solver is a structural analysis solver that performs numerical analysis on the FEM model represented by the model information 111.
  • the CAE post processor visualizes the analysis results obtained by the CAE solver.
  • the force balance equation is a simultaneous linear equation.
  • K is a stiffness matrix.
  • u is a vector representing the displacement of each node 103 for each degree of freedom.
  • u is an unknown vector obtained by analysis.
  • f is a vector representing the load.
  • f is a load value corresponding to the degree of freedom of each node 103 when the load condition is set.
  • the direct method is, for example, a solution method that directly obtains a solution by triangulating a coefficient matrix.
  • Examples of the direct method include a multifrontal method.
  • the iterative method is a method for obtaining an approximate solution by iterative calculation.
  • As an iterative method for example, there is a CG method with preprocessing.
  • the direct method generally requires more memory than the iterative method.
  • the direct method has a characteristic of simultaneous linear equations, and if the non-zero region increases, more calculation amount and memory amount are required, and analysis time may be longer.
  • the model is unstable, such as when the constraint on the model 101 by the constraint condition is insufficient, the solution may not easily converge, and the analysis time may be long.
  • the structural analysis apparatus 100 evaluates the model scale based on the information representing the analysis target model 101, and solves the solution by simultaneous linear equations using FEM according to the model scale as a direct method and an iterative method. Select from. As a result, the amount of memory used can be estimated, so that analysis by the direct method is performed without running out of memory, and the analysis time is shortened.
  • the structural analysis apparatus 100 acquires model information 111.
  • the model information 111 includes information representing the shape of the analysis target model 101, for example. A detailed example of the model information 111 will be described with reference to FIG.
  • the structural analysis apparatus 100 evaluates the scale of the analysis target model 101 based on the acquired model information 111.
  • the structural analysis apparatus 100 evaluates the model scale based on, for example, the number of elements and / or the number of nodes and / or the number of degrees of freedom of analysis calculated from the model information 111, the mesh size, and the element type.
  • the structural analysis apparatus 100 selects whether the algorithm for solving the simultaneous linear equations of the structural analysis solver using the finite element method of the model 101 is the direct method or the iterative method based on the evaluated model scale. For example, the structural analysis apparatus 100 selects the iterative method as the model size increases.
  • the structure analysis apparatus 100 selects whether the algorithm is a direct method or an iterative method based on the evaluated model scale and the amount of memory used at the time of analysis.
  • the structural analysis apparatus 100 performs the direct method when, for example, the amount of memory used when analyzing the model scale of the model 101 using the direct method is less than the memory usage that can be used at the time of analysis. select.
  • the structural analysis apparatus 100 selects the iterative method when, for example, the amount of memory used when analyzing the model size of the model 101 using the direct method is greater than or equal to the amount of memory used that can be used at the time of analysis. . In the example of FIG. 1, an iterative method is selected.
  • the estimated amount of memory used when analyzing a model can be calculated based on the model size calculated from the element type, the number of elements, the number of nodes, and the like.
  • the amount of memory to be used for one element may be prepared in advance in a storage device or the like. Then, the structural analysis apparatus 100 may calculate the memory amount to be used according to the number of elements of the analysis target model 101.
  • the direct method analysis is generally more accurate than the iterative method analysis, and it is understood that the direct method analysis will not take time, the direct method is used to improve the accuracy. There is. However, as described above, in the case of analysis by the direct method, if the memory usage is insufficient, the analysis takes time. For this reason, by selecting the analysis method according to the model size, the memory usage can be estimated based on the model size, so the analysis by the direct method can be performed without running out of memory, improving the accuracy and analyzing time. Can be shortened.
  • FIG. 2 is an explanatory diagram showing another example of the operation of the structural analysis apparatus according to the present invention.
  • the structural analysis apparatus 100 is a computer that performs structural analysis using a finite element method on the FEM model to be analyzed.
  • the structural analysis apparatus 100 shown in FIG. 2 and the structural analysis apparatus 100 shown in FIG. 1 may be the same apparatus or different apparatuses, and thus the same reference numerals are given.
  • FIG. 2 and the example shown in FIG. 1 detailed description of the same configuration and function will be omitted.
  • the structural analysis apparatus 100 evaluates the strength of the boundary condition of the model 101, and according to the strength of the boundary condition, solves the simultaneous linear equation using FEM as a direct method and an iterative method. And select from. Thereby, the analysis time can be shortened.
  • the boundary condition is a condition set at the boundary to represent the load state of the structure.
  • boundary conditions There are two types of boundary conditions, for example, constraint conditions and load conditions.
  • the constraint condition is a condition such as which part is fixed how and how much forced displacement is given to the structure to be analyzed.
  • the load condition is a condition such as how much force is applied to which part in which direction. In the case where the boundary condition is strong, for example, the constraint condition is strong, it can be determined whether the calculation model is stable due to the boundary constraint by the constraint condition, and the analysis time can be shortened.
  • the structural analysis apparatus 100 acquires the boundary condition of the analysis target model 101.
  • the structural analysis apparatus 100 evaluates the strength of the boundary condition of the analysis target model 101 based on the acquired boundary condition.
  • the strength of the boundary condition includes, for example, at least one of a constraint condition strength and a load condition weakness.
  • Constraint type includes the type and direction of the constraint, the position and area to be constrained, and the like.
  • the strength of the constraint condition is, for example, the strength with which the analysis target model 101 is constrained according to the constraint condition.
  • the strength of the constraint condition is evaluated based on, for example, the area of a region corresponding to the range constrained on the surface of the analysis target model 101 by the constraint condition as described later.
  • As the load condition the type of load, the magnitude and direction of the load, the position and area where the load is applied, and the like are set.
  • the weakness of the load condition is, for example, the reciprocal of the magnitude of the influence of the load condition on the rigidity.
  • the weakness of the load condition is evaluated based on, for example, the reciprocal of the magnitude of the influence of the load condition on the rigidity.
  • the structural analysis apparatus 100 selects, for example, whether the algorithm for solving the linear equations of the structural analysis solver using the finite element method of the model 101 is a direct method or an iterative method based on the strength of the evaluated boundary condition. To do. For example, the structural analysis apparatus 100 selects the iterative method as the strength of the boundary condition increases.
  • the structural analysis apparatus 100 selects the iterative method when the strength of the constraint condition is greater than or equal to the threshold value, and selects the direct method when the strength of the constraint condition is less than the threshold value.
  • the degree of stability of the boundary can be determined based on the strength of the constraint condition, and the analysis time can be shortened.
  • FIG. 3 is an explanatory diagram showing elements and nodes used for structural analysis of the finite element method.
  • the structure is, for example, a continuous body.
  • the CAE preprocessor creates an FEM model by approximating a three-dimensional model of a structure created in a three-dimensional simulation space.
  • the element types include a three-dimensional solid element, a two-dimensional shell element, a one-dimensional beam element and a truss element, and a zero-dimensional mass element.
  • Solid element shapes include hexahedron, pentahedron, and tetrahedron.
  • the shape of the shell element includes a quadrangle and a triangle.
  • the shape of the beam element or truss element is a line.
  • the shape of the mass element is a point.
  • the element shape is a tetrahedron
  • there are a plurality of node defining methods such as a primary element and a secondary element.
  • the nodes included in the element are the vertices of the tetrahedron.
  • the nodes are the vertices of the tetrahedron and the center point of the line connecting the vertices.
  • the number of nodes constituting an element differs depending on the element shape and the node definition method.
  • description will be made using a tetrahedral secondary element as an element type.
  • the solid element is also referred to as a solid element.
  • the shell element is also referred to as a plate element, for example.
  • the beam element is also referred to as a beam element.
  • the number of analysis degrees of freedom of the nodes described later differs.
  • FIG. 4 is a block diagram illustrating a hardware configuration example of the structural analysis apparatus.
  • the structural analysis apparatus 100 includes a CPU (Central Processing Unit) 401, a ROM (Read Only Memory) 402, a RAM (Random Access Memory) 403, a disk drive 404, and a disk 405.
  • the structural analysis apparatus 100 includes an I / F (Interface) 406, a keyboard 407, a mouse 408, and a display 409. Further, the CPU 401, the ROM 402, the RAM 403, the disk drive 404, the I / F 406, the keyboard 407, the mouse 408, and the display 409 are connected by a bus 400, respectively.
  • the CPU 401 controls the entire structure analysis apparatus 100.
  • the ROM 402 stores programs such as a boot program.
  • the RAM 403 is used as a work area for the CPU 401.
  • the disk drive 404 controls reading / writing of data with respect to the disk 405 according to the control of the CPU 401.
  • the disk 405 stores data written under the control of the disk drive 404. Examples of the disk 405 include a magnetic disk and an optical disk.
  • the I / F 406 is connected to a network 410 such as a LAN (Local Area Network), a WAN (Wide Area Network), or the Internet through a communication line, and is connected to another device via the network 410.
  • the I / F 406 controls an internal interface with the network 410 and controls input / output of data from an external device.
  • a modem or a LAN adapter may be employed as the I / F 406.
  • the keyboard 407 and the mouse 408 are interfaces for inputting various data by user operation.
  • a display 409 is an interface that outputs data in accordance with an instruction from the CPU 401.
  • the structural analysis apparatus 100 may be provided with an input device that captures images and moving images from a camera and an input device that captures audio from a microphone. Although not shown, the structure analysis apparatus 100 may be provided with an output device such as a printer.
  • FIG. 5 is a block diagram illustrating a functional configuration example of the structural analysis apparatus.
  • the structural analysis apparatus 100 includes an acquisition unit 501, a first evaluation unit 502, a second evaluation unit 503, and a selection unit 504.
  • the processing of the control unit from the acquisition unit 501 to the selection unit 504 is coded in a program stored in a storage device such as the ROM 402, the RAM 403, and the disk 405 that can be accessed by the CPU 401 shown in FIG.
  • the CPU 401 reads the program from the storage device and executes the process coded in the program. Thereby, the process of a control part is implement
  • the processing result of the control unit is stored in a storage device such as the ROM 402, the RAM 403, and the disk 405, for example.
  • the acquisition unit 501 acquires the model information 111.
  • the model information 111 is information representing the analysis target model.
  • the model information 111 includes information representing the shape of the analysis target model, for example.
  • Information representing the shape of the model includes element information of elements included in the model, nodes included in the elements, position information of the nodes, and the like.
  • the information representing the shape of the analysis target model may be information such as volume and surface area, for example.
  • the acquired model information 111 may include information indicating the material of the model, for example.
  • the acquired model information 111 may include, for example, information on boundary conditions set in the model at the time of analysis.
  • the acquired model information 111 may include, for example, information on a result output item.
  • the information of the output item of the result is information indicating which analysis result is calculated and output, such as displacement, stress, reaction force, and the like.
  • FIG. 6 is an explanatory diagram showing an example of model information.
  • 6 (1) to 6 (3) show, as model information 111, element information of elements included in the model, nodes included in the elements, position information of the nodes, and the like.
  • model information 111 includes information representing the shape of the model, information representing the material of the model, information representing boundary conditions set in the model at the time of analysis, and the like.
  • the boundary condition is a condition set at the boundary to represent the load state of the structure.
  • boundary conditions There are two types of boundary conditions, for example, constraint conditions and load conditions.
  • the constraint condition is a condition such as which part is fixed how and how much forced displacement is given to the structure to be analyzed.
  • the load condition is a condition such as how much force is applied to which part in which direction.
  • the model information 111 has model ID and element ID fields. Identification information for identifying a model is set in the model ID field. In the element ID field, identification information for identifying an element included in the model to be analyzed is set. In the example of FIG. 6 (1), the model ma has elements en1 to en3, for example.
  • the model information 111 has fields for element ID and node ID. Identification information for identifying the element is set in the element ID field. Identification information for identifying the node included in the element is set in the node ID field.
  • the element en1 has nodes nn1 to nn4.
  • the model information 111 has a field of node ID and node position information. Identification information for identifying the node is set in the node ID field. In the node position information field, node position information is set.
  • the node position information is, for example, the coordinate values of the x-axis, y-axis, and z-axis of the node.
  • the structure analysis apparatus 100 can specify the shape of the model based on the model information 111 shown in FIGS. 6 (1) to 6 (3).
  • the structural analysis apparatus 100 can calculate the volume and surface area of the model.
  • the model information 111 includes model ID, volume, mass, main moment of inertia, number of nodes, element type, mesh size, number of elements, number of analysis degrees of freedom, minimum sphere volume, It has fields such as surface area and boundary conditions.
  • Identification information for identifying a model is set in the model ID field.
  • the volume of the model is set in the volume field.
  • the mass of the model is set in the mass field.
  • the main inertia moment of the model is set in the main inertia moment field.
  • the volume, mass and main moment of inertia are information determined from the material and shape of the model.
  • degrees of freedom for example, the total number of degrees of freedom of each node included in the model is set.
  • the final number of degrees of freedom of analysis cannot be determined without completing the assembly of simultaneous linear equations performed within the CAE solver. For this reason, as the number of degrees of freedom in analysis, for example, an approximate value based on the number of nodes ⁇ the number of degrees of freedom per node may be set.
  • the element length of one representative element for mesh division is set.
  • the case where all the elements are the same element type will be described as an example, and therefore the mesh size is one.
  • the number of nodes included in the model is set in the number of nodes field.
  • the number of elements included in the model is set.
  • a predicted value can be calculated based on the mesh size, the volume of the shape, and the element type.
  • the number of nodes can be calculated as an approximate prediction value based on the number of elements and the element type, and the accuracy of prediction can be increased by taking the model shape into consideration.
  • the number of elements and the number of nodes can be specified based on the model information 111 when the model information 111 shown in FIGS. 6 (1) to (3) already exists.
  • the minimum sphere volume covering the model is set.
  • there is a field of the minimum sphere volume but the present invention is not limited to this.
  • the surface area of the model is set in the surface area field.
  • the boundary condition field includes, for example, a constraint condition, a constraint area, a constraint degree of freedom, and a load condition field.
  • a constraint condition is set in the constraint condition field.
  • the constraint condition include one-side end surface fixing, both-end fixing, and all-around fixing.
  • the constraint area field for example, an area that is constrained in accordance with a constraint condition among the surface areas of the model is set.
  • the number of degrees of freedom to be constrained for example, the number of degrees of freedom to be constrained among the degrees of freedom per node is set. For example, when there are three directions of freedom per node, the number of directions to be restricted among the three directions is the number of degrees of freedom of restriction.
  • the load condition is set in the load condition field. Examples of the load condition include a gravity load in each direction.
  • the structural analysis apparatus 100 evaluates a model scale based on information representing the analysis target model 101, and describes an example in which a solution method using simultaneous linear equations using FEM is selected from a direct method and an iterative method according to the model scale. To do.
  • the first evaluation unit 502 evaluates the model scale of the model based on the model information 111 acquired by the acquisition unit 501.
  • the model size of the model may be, for example, the number of elements, the number of nodes, or the number of degrees of freedom of analysis.
  • the model scale may be a bandwidth, a skyline, a profile, the number of non-zero terms, the number of fill-ins, or the like related to the coefficient matrix, or may be a value calculated by combining these.
  • the first evaluation unit 502 evaluates the model scale based on, for example, the number of elements and / or the number of nodes and / or the number of degrees of freedom of analysis calculated from the model information 111 and the mesh size.
  • the model information 111 includes information representing the shape of the model, for example, the model volume
  • the first evaluation unit 502 calculates the number of elements based on the model volume, the mesh size, and the element type. calculate.
  • the number of elements calculated here is a predicted value.
  • the first evaluation unit 502 calculates the number of nodes based on the model volume, mesh size, and element type, for example.
  • the number of nodes calculated here is a predicted value.
  • the first evaluation unit 502 calculates the number of degrees of freedom in analysis based on, for example, the model volume, mesh size, element type, and constraint conditions.
  • the number of degrees of freedom of analysis calculated here is a predicted value.
  • the first evaluation unit 502 identifies the number of elements, the number of nodes, and the number of degrees of freedom of analysis from the model information 111. May be.
  • the analysis degree of freedom is, for example, the sum of the degrees of freedom of each node included in the model.
  • the number of degrees of freedom in analysis cannot be determined unless the assembly of simultaneous linear equations performed within the CAE solver is completed. Therefore, the number of degrees of freedom in analysis can be calculated as a predicted value by the number of nodes ⁇ the number of degrees of freedom per node.
  • the number of elements will be described as an example of the model scale index value.
  • the number of degrees of freedom per node is determined by the element type, but may be included in the model information 111, for example.
  • the selection unit 504 selects whether the algorithm for solving the simultaneous linear equations of the structural analysis solver using the model finite element method is a direct method or an iterative method based on the evaluation result of the model scale. For example, the selection unit 504 selects the iterative method as the model size increases.
  • the selection unit 504 acquires, for example, the amount of free memory.
  • the free memory amount is the amount of free space in a storage device such as the RAM 403 and the disk 405, for example.
  • the selection unit 504 selects the direct method or the iterative method based on the amount of free memory and the evaluated model size. The larger the model size, the more memory usage.
  • a threshold for the amount of free memory is prepared according to the model size, and the selection unit 504 uses the direct method if the acquired amount of free memory is equal to or greater than the threshold according to the model size.
  • the selection unit 504 uses an iterative method when the acquired amount of free memory is less than a threshold corresponding to the model scale.
  • FIG. 7 is an explanatory diagram showing an example of a memory capacity threshold value according to the number of elements.
  • the threshold value table 700 has, for example, threshold values for the memory amount for each number of elements.
  • the number of elements is, for example, 1 to 3000, 3001 to 5000, 5001 to 7000, 7001 to 10,000. For example, if the model size is 4500, a threshold mt2 corresponding to 3001 to 5000 is selected.
  • the larger the model size the larger the threshold for the amount of memory.
  • the larger the model size the larger the threshold value, so the iterative method is used. For this reason, a shortage of memory during analysis can be avoided, and the analysis time can be shortened.
  • the direct method can be used if the amount of free memory is large, so that the analysis accuracy can be improved.
  • FIG. 8 is a flowchart illustrating an example of a selection processing procedure based on the model scale by the structural analysis apparatus.
  • the structural analysis apparatus 100 acquires model information 111 representing the analysis target model (step S801). Then, the structural analysis apparatus 100 evaluates the model scale (step S802).
  • the structural analysis apparatus 100 confirms the amount of free memory (step S803).
  • the structural analysis apparatus 100 selects a solution according to the model scale and the amount of free memory (step S804). Then, the structural analysis apparatus 100 performs numerical analysis using the selected solution (step S805), and ends a series of processes.
  • the second evaluation unit 503 evaluates the strength of the boundary condition based on the acquired model information 111 and the boundary condition.
  • the boundary condition includes a constraint condition and a load condition.
  • the strength of the boundary condition includes, for example, at least one of a constraint condition strength and a load condition weakness.
  • Constraint conditions include, for example, the type and direction of constraint, the position and region to be restrained, and the like.
  • the strength of the constraint condition is, for example, the strength with which the analysis target model is constrained according to the constraint condition.
  • the strength of the constraint condition is evaluated based on, for example, the area of a region corresponding to a range that is constrained on the surface of the model to be analyzed by the constraint condition as described later.
  • As the load conditions for example, the type of load, the magnitude and direction of the load, and the position and region where the load is applied are set.
  • the weakness of the load condition is evaluated based on, for example, the reciprocal of the magnitude of the influence of the load condition on the rigidity.
  • the second evaluation unit 503 evaluates the strength of the constraint condition based on the constraint area of the model and the constraint freedom.
  • the restrained area is an area of a portion restrained in the model according to the restraint condition. More specifically, the second evaluation unit 503 calculates the strength of the constraint condition using the following formula (1).
  • Strength of restraint condition restraint area / surface area ⁇ (degree of freedom of restraint per node / number of degrees of freedom per node) Formula (1)
  • the sum of these may be used as the strength of the constraint conditions of the entire model.
  • the selection unit 504 selects, based on the evaluation result of the constraint condition strength, whether the algorithm for solving the linear equations of the structural analysis solver using the finite element method of the model is a direct method or an iterative method. For example, the selection unit 504 causes the iterative method to be selected as the constraint condition increases. Specifically, for example, the selection unit 504 selects the iterative method when the strength of the constraint condition is greater than or equal to a threshold value, and selects the direct method when the strength of the constraint condition is less than the threshold value.
  • the threshold value is set in advance and stored in a storage device such as the ROM 402, the RAM 403, and the disk 405.
  • the threshold setting method is described as being performed by the structural analysis apparatus 100.
  • the present invention is not limited to this, and a threshold different from the structural analysis apparatus 100 may be determined and provided to the structural analysis apparatus 100. Good.
  • the structural analysis apparatus 100 prepares a plurality of constraint conditions for a model, and calculates the strength of the constraint conditions for each of the constraint conditions.
  • the structural analysis apparatus 100 performs numerical analysis by a direct method and an iterative method on the model for each of a plurality of constraint conditions.
  • the structural analysis device 100 obtains the strength of the constraint condition when the performance ratio of the analysis time by the direct method and the analysis time by the iterative method is 1, by linear interpolation, and sets it as a threshold value.
  • FIG. 9 is an explanatory diagram showing a model example used for deriving a threshold value.
  • the threshold value derivation will be described by taking the model m1 as an example.
  • the dimension of the model m1 is, for example, a rectangular parallelepiped of 50 ⁇ 100 ⁇ 1000 [mm].
  • constraint condition 1 is cantilever.
  • constraint condition 2 is fixed at both ends.
  • constraint condition 3 is fixed all around.
  • the strength of the constraint condition and the analysis time for each of the case where the model m1 is mesh-divided so that the number of elements is about 5000 and the case where the model is divided so that the number of elements is about 10000
  • the threshold value is set as follows.
  • FIG. 10 is an explanatory diagram showing an example of deriving a threshold value when the number of elements is 5000.
  • (1) to (4) show numerical values used for deriving the strength of the constraint condition
  • (5) and (6) show the analysis time.
  • the structural analysis apparatus 100 determines the strength of the constraint condition by “constraint area / surface area ⁇ (number of degrees of freedom of constraint per node / number of degrees of freedom per node)”. calculate.
  • the restrained area is shown in (3).
  • the surface area is shown in (1).
  • the degree of freedom of constraint per node is shown in (4).
  • the number of degrees of freedom per node is shown in (2).
  • the strength of the constraint condition is “0.016” by “(5000/310000) ⁇ (3/3)”.
  • the strength of the constraint condition is “0.032” by “(10000/310000) ⁇ (3/3)”.
  • the structural analysis apparatus 100 obtains an analysis time, for example, by performing numerical analysis by a direct method and numerical analysis by an iterative method for each constraint condition. Then, the structural analysis apparatus 100 obtains the ratio of the analysis time by the direct method and the analysis time by the iterative method as the performance ratio of the analysis method.
  • the analysis time by the direct method is shown in (5).
  • the analysis time by the iterative method is shown in (6).
  • the structural analysis apparatus 100 obtains the strength of the constraint condition when the performance ratio is 1 by linearly interpolating the constraint condition strength of the constraint condition 2 and the constraint condition strength of the constraint condition 3. Then, the structural analysis apparatus 100 sets the strength of the constraint condition when the performance ratio is 1 as the threshold value. For example, not only linear interpolation but also Lagrange interpolation, spline interpolation, or least square method may be used.
  • the analysis time by the direct method is shorter than the analysis time by the iterative method.
  • the analysis time by the iterative method is shorter than the analysis time by the direct method.
  • the selection unit 504 selects the iterative method if the strength of the constraint condition calculated by the second evaluation unit 503 is equal to or greater than the threshold, and selects the direct method if it is less than the threshold. By doing so, the analysis time can be shortened.
  • FIG. 11 is an explanatory diagram showing an example of deriving a threshold value when the number of elements is 10,000.
  • (1) to (4) show numerical values used for deriving the strength of the constraint condition
  • (5) and (6) show the analysis time.
  • the strength of the constraint condition is the same as in the example of FIG.
  • the performance ratio is 1 between the performance ratio of the constraint condition 2 and the performance ratio of the constraint condition 3.
  • FIG. 11 shows an example in which the threshold value is obtained in the same manner as in FIG.
  • the selection unit 504 shown in FIG. 5 directly calculates an algorithm for solving the linear equations of the structural analysis solver using the model finite element method based on the evaluation result of the model scale and the evaluation result of the strength of the constraint condition. Or iterative method. Since the evaluation of the model scale is the evaluation by the first evaluation unit 502 described above, detailed description thereof is omitted.
  • the selection unit 504 changes the selection criterion based on the model size when selecting whether to use the direct method or the iterative method based on the strength of the constraint condition. Specifically, the selection unit 504 selects, for example, whether the solution method is a direct method or an iterative method based on the model scale, and determines whether the solution method is a direct method or an iterative method based on the strength of the constraint condition. Select whether to use the law. Then, the selection unit 504 selects, based on the selection result based on the model size and the selection result based on the strength of the constraint condition, for example, whether the solution method is to be a direct method or an iterative method in the end. May be.
  • the selection criterion is changed based on the model scale will be described.
  • FIG. 12 is an explanatory diagram showing an example of a threshold value corresponding to the model scale.
  • the threshold table 1200 has a threshold for each number of elements, for example.
  • the number of elements is given as an example of model scale evaluation.
  • the number of elements is not limited to this, and may be the number of nodes or the number of degrees of freedom of analysis as described above.
  • the number of elements is, for example, 3000, 5000, 7500, 10000.
  • the iterative method is selected as the number of elements is larger.
  • the selection unit 504 acquires a threshold value corresponding to the number of elements from the threshold value table 1200, and compares the calculated strength of the constraint condition with the acquired threshold value. Specifically, for example, when the number of elements of the analysis target model is 4000 or less, the selection unit 504 acquires a threshold corresponding to the number of elements of 3000. Specifically, for example, when the number of elements of the analysis target model is greater than 4000 and 6250 or less, the selection unit 504 acquires a threshold corresponding to the number of elements of 5000. Thus, a range of the number of elements may be provided. The selection unit 504 selects the iterative method when the calculated constraint strength is equal to or greater than the acquired threshold, and selects the direct method when the calculated constraint strength is less than the acquired threshold. .
  • FIG. 13 is an explanatory diagram showing an example of calculating the strength of constraint conditions and an example of selecting a solution.
  • the analysis target model has, for example, a size of 50 ⁇ 100 ⁇ 500.
  • the first evaluation unit 502 evaluates the model size of the analysis target model, for example.
  • the number of elements is used as the model scale index value.
  • the second evaluation unit 503 evaluates the strength of the constraint condition based on the constraint condition of the analysis target model, for example.
  • the selection unit 504 acquires a threshold value according to the model scale evaluated from the threshold value table 1200, for example.
  • the selection unit 504 acquires a threshold corresponding to the number of elements of 5000 from the threshold table 1200, for example.
  • the threshold corresponding to 5000 elements is “0.202” using the example of FIG.
  • the selection unit 504 determines that the strength of the constraint condition is less than the threshold value, and selects the direct method. If the iterative method is used when the model is not sufficiently constrained by the constraint conditions, the solution does not converge easily. Therefore, the analysis time can be shortened by selecting the direct method when the strength of the constraint condition of the model is weak such that the strength of the constraint condition is less than the threshold value.
  • FIG. 14 is a flowchart illustrating an example of a selection procedure based on the strength of the constraint condition by the structural analysis apparatus.
  • the structural analysis apparatus 100 acquires model information 111 of the analysis target model (step S1401).
  • the structure analysis apparatus 100 evaluates the strength of the boundary condition (step S1402). Then, the structural analysis apparatus 100 evaluates the model scale (step S1403). Next, the structure analysis apparatus 100 acquires a threshold value corresponding to the model scale (step S1404). The structural analysis apparatus 100 compares the boundary condition strength with a threshold value (step S1405).
  • step S1405 If the strength of the boundary condition is less than the threshold value (step S1405: less than the threshold value), the structural analysis apparatus 100 selects the direct method (step S1406) and proceeds to step S1408. On the other hand, when the strength of the boundary condition is greater than or equal to the threshold (step S1405: greater than or equal to the threshold), the structure analysis apparatus 100 selects an iterative method (step S1407), and proceeds to step S1408.
  • the structural analysis apparatus 100 performs numerical analysis by the selected solution (step S1408) and ends a series of processing.
  • FIG. 15 is an explanatory diagram showing a system application example.
  • the system 1500 includes, for example, a client terminal device 1501 and a data center 1502.
  • the client terminal device 1501 may be plural as shown in FIG.
  • the data center 1502 includes a plurality of servers 1503.
  • the client terminal device 1501 and the data center 1502 are connected via the network 410.
  • the management method of the server 1503 in the data center 1502 is not particularly limited.
  • the structure analysis apparatus 100 may be realized by a client terminal device 1501 and a plurality of servers 1503.
  • the acquisition unit 501, the first evaluation unit 502, and the second evaluation unit 503 may be realized by the client terminal device 1501
  • the selection unit 504 may be realized by the server 1503.
  • the server 1503 executes, for example, CAD and CAE.
  • the client terminal device 1501 acquires the model information 111 created by CAE and evaluates the strength of the constraint condition. Then, the client terminal apparatus 1501 transmits the strength of the evaluated constraint condition to the server 1503. Then, the server 1503 selects one of the iterative method and the direct method based on the strength of the constraint condition, and executes numerical analysis using the selected solution.
  • the server 1503 may register past analysis results in a database. Then, the server 1503 may select a solution based on the strength of the constraint condition received from the client terminal device 1501.
  • FIG. 16 is an explanatory diagram showing an example of an analysis database.
  • the analysis database 1600 includes, for example, fields of model ID, boundary condition strength, model scale, direct method A analysis time, direct method B analysis time, iterative method C analysis time, and iterative method D analysis time. . By setting information in each field, it is stored as records (1601-1, 1601-2).
  • the analysis database 1600 is realized by a storage device such as the ROM 402, the RAM 403, and the disk 405, for example.
  • model ID field identification information for identifying a model is set.
  • boundary condition strength field the boundary condition strength of the model indicated by the model ID is set. For example, the strength of the constraint condition is set.
  • model size field the model size of the model indicated by the model ID is set. For example, the number of elements is set.
  • an analysis time when a numerical analysis is performed using the direct method A is set.
  • an analysis time when a numerical analysis is performed using the direct method B is set.
  • an analysis time when a numerical analysis is performed using the iterative method C is set.
  • an analysis time when a numerical analysis is performed using the iterative method D is set.
  • an analysis time field is prepared for each type of the direct method or the iterative method. Note that the method is not limited to the direct method or the iterative method, and there may be an analysis field of another method.
  • FIG. 17 is an explanatory diagram showing an example of selecting a solution based on the results in the system.
  • the client terminal device 1501 transmits evaluation results such as the strength of constraint conditions and the model scale to the server 1503.
  • the server 1503 searches the analysis database 1600 for a record 1601 having an evaluation result closest to the evaluation result received from the client terminal device 1501. Specifically, the server 1503 compares the evaluation result with the evaluation result in the record 1601 included in the analysis database 1600, scores the similarity, and detects the record 1601 having the closest evaluation result. Good.
  • the server 1503 selects the solution with the shortest analysis time among the analysis times of the solutions included in the detected record 1601 as the solution of the analysis target model. In the example of FIG. 17, iterative method C is selected. For example, the server 1503 determines a calculation resource based on the selected solution and model size. In the example of FIG. 1, the analysis is performed using the server 1503-1 and the server 1503-2. An example of determining a calculation resource will be described with reference to FIG.
  • the server 1503-1 and the server 1503-2 perform numerical analysis on the analysis target model using the selected iteration method C by parallel processing.
  • the server 1503 transmits the analysis result to the client terminal device 1501. Further, the server 1503 registers the evaluation result and the analysis time in the analysis database 1600 as a record 1601.
  • the server 1503 performs numerical analysis by a solution that is not selected as a solution of the analysis target model when computer resources are available, and registers the analysis time in the analysis database 1600.
  • the server 1503 may select a solution by comparing the strength of the boundary condition and the threshold corresponding to the model scale as described above. Good.
  • FIG. 18 is an explanatory diagram showing an example of managing the server state and the memory amount.
  • the table 1800 has information for managing the status and resources of the server 1503 in the data center 1502.
  • the table 1800 has fields such as a server ID, a state, a memory amount, and the number of CPUs.
  • identification information for identifying the server 1503 is set.
  • Information indicating the status of the server 1503 is set in the status field.
  • the state of the server 1503 is represented by “in use” and “unused”. In more detail, the state of the server 1503 may be represented by the number of CPUs in use.
  • the memory amount field the memory amount of the RAM 403 or the disk 405 included in the server 1503 is set.
  • the number of CPUs field for example, the number of processors included in the server 1503 is set.
  • the server 1503-1 to the server 1503-3 have the same amount of memory and the same number of CPUs, but the server 1503-4 has a larger amount of memory and the number of CPUs than the servers 1503-1 to 1503-3. Take the case as an example.
  • the model to be analyzed is from x to z. Assume that the model size of the analysis target is about the same.
  • the constraint condition strength of the analysis target model x is ax
  • the constraint condition strength of the analysis target model y is ay
  • the constraint condition strength of the analysis target model z is az.
  • the relationship between the strengths of the constraint conditions is ax> threshold> ay> az.
  • the iterative method is selected as the solution method of the analysis target model x
  • the direct method is selected as the solution method of the analysis target model y
  • the direct method is selected as the solution method of the analysis target model z.
  • the server 1503 selects a computing resource with a large amount of memory. For example, the server 1503 selects the server 1503-4 as a calculation resource. For example, when the iterative method is selected, the server 1503 selects one of the servers 1503-1 to 1503-3. For example, when a plurality of iterative methods can be selected, the server 1503 performs simultaneous processing by a plurality of servers 1503 among the servers 1503-1 to 1503-3 depending on the type of the iterative method. Also good. Further, the server 1503 may select a plurality of computing resources in the case of a solution that allows parallel processing using a plurality of servers. Further, the server 1503 may estimate in advance a calculation resource such as a memory amount necessary for efficient calculation from the model information, and may select an appropriate calculation resource according to the result.
  • a calculation resource such as a memory amount necessary for efficient calculation from the model information
  • the server 1503-4 is selected as a calculation resource.
  • the server 1503-1 is selected as a calculation resource.
  • the server 1503-2 is selected as a calculation resource.
  • the selection unit 504 selects a calculation resource for each analysis target model.
  • FIGS. 19 and 20 are flowcharts showing an example of the structure analysis processing procedure by the system.
  • the system 1500 includes the client terminal device 1501 and the server 1503 as described above. Therefore, the processing illustrated in FIGS. 19 and 20 may be performed by the server 1503 or the client terminal device 1501.
  • the system 1500 acquires the model information 111 of the analysis target model by the client terminal device 1501 (step S1901). Next, the system 1500 uses the client terminal device 1501 to evaluate the model size of the analysis target model (step S1902). The system 1500 evaluates the strength of the boundary condition by the client terminal device 1501 (step S1903).
  • the system 1500 searches the record 1501 including the evaluation result closest to the evaluation result of the analysis target model by the server 1503 from the analysis database 1600 (step S1904).
  • the system 1500 selects the solution with the shortest analysis time from the record 1601 searched by the server 1503 (step S1905).
  • the system 1500 performs numerical analysis by the solution specified by the server 1503 (step S1906).
  • the system 1500 outputs an analysis result by the server 1503 (step S2001).
  • the server 1503 stores the evaluation result and the analysis time in the analysis database 1600 (step S2002).
  • the system 1500 determines whether or not the load on the computer environment is high (step S2003). When the load on the computer environment is high (step S2003: Yes), the system 1500 returns to step S2003.
  • the high load on the computer environment means, for example, a state where there is no server 1503 without processing.
  • step S2003 determines whether the server 1503 has executed analysis with all the solutions.
  • step S2004 determines whether the server 1503 has executed analysis with all the solutions.
  • step S2004 executes numerical analysis using a solution having the shortest analysis time among the searched records 1601 among the unexecuted solutions.
  • Step S2005 The system 1500 stores the analysis time in the analysis database 1600 by the server 1503 (step S2006), and returns to step S2003.
  • step S2004 when it is determined that the analysis has been executed for all the solutions (step S2004: Yes), the system 1500 ends the series of processes.
  • the structural analysis apparatus 100 evaluates the model size based on the information representing the model to be analyzed, and selects a solution method using simultaneous linear equations using FEM from the direct method and the iterative method according to the model size. To do.
  • the amount of memory used can be estimated based on the model scale, so that the analysis by the direct method can be performed without running out of memory. Therefore, the analysis time can be shortened while improving accuracy.
  • the structural analysis apparatus 100 evaluates the model scale based on the number of elements and / or the number of nodes and / or the number of degrees of freedom of analysis calculated from the model information, the mesh size, and the element type.
  • the model scale can be easily evaluated, and the time required for the evaluation can be shortened.
  • the structural analysis apparatus 100 evaluates the strength of the boundary condition based on the boundary condition of the model, and selects a solution method using simultaneous linear equations using FEM from the direct method and the iterative method according to the strength of the boundary condition. . Thereby, it can be determined whether the boundary constraint by the boundary condition is stable, and a solution method more suitable for the model can be selected. Therefore, the analysis time can be shortened.
  • the structural analysis apparatus 100 performs the evaluation of the strength of the boundary condition based on the constraint area and the constraint freedom included in the constraint condition of the model boundary conditions. Thereby, it is possible to easily determine how much the boundary constraint is stabilized, and it is possible to reduce the time required for the evaluation.
  • the structural analysis apparatus 100 evaluates the model scale based on the model information, evaluates the strength of the boundary condition based on the boundary condition, and uses simultaneous linear equations using FEM according to the model scale and the strength of the boundary condition. Choose between direct and iterative solutions. Thereby, the possibility that the solution does not converge can be reduced, and the solution method can be selected so that the amount of memory used is sufficient, and the analysis time can be shortened.
  • the structural analysis apparatus 100 evaluates the model scale based on the number of elements and / or the number of nodes and / or the number of degrees of freedom of analysis calculated from the model information, the mesh size, and the element type. As a result, the model scale can be easily evaluated, and the time required for the evaluation can be shortened. In addition, the structural analysis apparatus 100 evaluates the strength of the constraint condition based on the constraint area and the constraint freedom included in the model constraint condition. Thereby, it is possible to easily determine how much the boundary constraint is stabilized, and it is possible to reduce the time required for the evaluation.
  • the structural analysis apparatus 100 evaluates the model scale of the model to be analyzed based on the model information, selects a numerical analysis algorithm for structural analysis using the finite element method of the model based on the model scale, Determine the computational resources used in the numerical analysis. As a result, it is possible to select a more appropriate solution and calculation resource in accordance with the characteristics of the numerical analysis of the model, and shorten the analysis time.
  • the structural analysis apparatus 100 evaluates the strength of the boundary condition based on the boundary condition of the model, selects a numerical analysis algorithm for the structural analysis using the finite element method of the model according to the strength of the boundary condition, Determine the computational resources used in the numerical analysis. As a result, it is possible to select a more appropriate solution and calculation resource in accordance with the characteristics of the numerical analysis of the model, and shorten the analysis time.
  • the structural analysis method described in this embodiment can be realized by executing a prepared structural analysis program on a computer such as a personal computer or a workstation.
  • the structural analysis program is recorded on a computer-readable recording medium such as a magnetic disk, an optical disk, or a USB (Universal Serial Bus) flash memory, and is executed by being read from the recording medium by the computer.
  • the structural analysis program may be distributed through a network such as the Internet.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

構造解析装置（１００）は、モデル情報（１１１）を取得する。構造解析装置（１００）は、取得したモデル情報（１１１）に基づき、解析対象モデル（１０１）の規模を評価する。構造解析装置（１００）は、モデル情報（１１１）とメッシュサイズと要素タイプとから算出された要素数及び／又は節点数に基づいてモデル規模の評価を行う。構造解析装置（１００）は、評価したモデル規模に基づいて、モデル（１０１）の有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する。

Description

構造解析方法、及び構造解析プログラム
　本発明は、構造解析方法、及び構造解析プログラムに関する。
　従来、有限要素法を用いて構造解析を行う技術が公知である。有限要素法を用いた構造解析では、例えば、連立一次方程式を解くことにより構造物の変位や応力を求める。この連立一次方程式の求解は、計算量が多く計算に時間がかかるため、計算機を用いて高速に解くための様々な技術がある。連立一次方程式の解法としては、例えば、大きな分類として直接法や反復法などが挙げられる。直接法は、例えば、係数行列を三角化することで解を直接求める解法であるのに対して、反復法は、反復計算により近似解を求める解法である。
　直接法や反復法に代表される連立一次方程式求解のアルゴリズムは、計算機をより効率的に使用して、より高速に解くための様々な技術開発が行われてきている。例えば、非特許文献１には疎行列向きの直接法のアルゴリズムが紹介されている。非特許文献２には並列処理に適した反復法のアルゴリズムとその計算例が紹介されている。
山本有作:"疎行列連立一次方程式の直接解法"，計算工学，Ｖｏｌ．１１，Ｎｏ．４，ｐｐ．１４－１８（２００６） 大原敏靖，川島康弘，藤野清次："Ｈｙｂｒｉｄ反復法ソルバによるＦＥＭ構造解析の処理時間短縮"，計算工学講演会論文集　Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅ　Ｃｏｎｆｅｒｅｎｃｅ　ｏｎ　Ｃｏｍｐｕｔａｔｉｏｎａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　ａｎｄ　Ｓｃｅｉｅｎｃｅ　１９，４ｐ，２０１４－０６
　しかしながら、連立一次方程式の解法には得手不得手があり、全ての問題を確実に、かつ最も早く解く解法は存在しない。したがって、従来、解析対象のモデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを特定の解法として数値解析を行うと、解析対象のモデルによっては他の解法よりも解析時間が長くなるという問題点がある。例えば、直接法は、連立一次方程式の特性により、非ゼロ領域が増えると演算量やメモリ量がより多く必要となり、解析時間が長くなる場合がある。一方、例えば、反復法は、連立一次方程式の特性により解が収束しにくいことがあり、解析時間が長くなる場合がある。
　１つの側面では、本発明は、解析時間の短縮化を図ることができる構造解析方法、及び構造解析プログラムを提供することを目的とする。
　本発明の一側面によれば、モデル情報を取得し、取得した前記モデル情報に基づいて、解析対象のモデルのモデル規模を評価し、前記モデル規模の評価結果に基づいて、前記モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する構造解析方法、及び構造解析プログラムが提案される。
　また、本発明の他の側面によれば、解析対象のモデルの境界条件を取得し、取得した前記境界条件に基づいて、前記境界条件の強さを評価し、前記境界条件の強さの評価結果に基づいて、前記モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式のアルゴリズムを直接法とするか反復法とするか選択する構造解析方法、及び構造解析プログラムが提案される。
　また、本発明の他の側面によれば、モデル情報と解析対象のモデルの境界条件とを取得し、取得した前記モデル情報に基づいて、前記モデルのモデル規模を評価し、取得した前記境界条件に基づき、前記境界条件の強さを評価し、前記モデル規模の評価結果及び前記境界条件の強さの評価結果に基づいて、前記モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する構造解析方法、及び構造解析プログラムが提案される。
　本発明の一態様によれば、解析時間の短縮化を図ることができる。
図１は、本発明にかかる構造解析装置による一の動作例を示す説明図である。 図２は、本発明にかかる構造解析装置による他の動作例を示す説明図である。 図３は、有限要素法の構造解析に用いる要素と節点を示す説明図である。 図４は、構造解析装置のハードウェア構成例を示すブロック図である。 図５は、構造解析装置の機能的構成例を示すブロック図である。 図６は、モデル情報例を示す説明図である。 図７は、要素数に応じたメモリ容量の閾値例を示す説明図である。 図８は、構造解析装置によるモデル規模に基づく選択処理手順例を示すフローチャートである。 図９は、閾値の導出に用いるモデル例を示す説明図である。 図１０は、要素数が５０００の場合の閾値の導出例を示す説明図である。 図１１は、要素数が１００００の場合の閾値の導出例を示す説明図である。 図１２は、モデル規模に応じた閾値例を示す説明図である。 図１３は、拘束条件の強さの算出例と解法の選択例を示す説明図である。 図１４は、構造解析装置による拘束条件の強さに基づく選択手順例を示すフローチャートである。 図１５は、システム適用例を示す説明図である。 図１６は、解析データベース例を示す説明図である。 図１７は、システムにおける実績に基づく解法の選択例を示す説明図である。 図１８は、サーバの状態とメモリ量の管理例を示す説明図である。 図１９は、システムによる構造解析処理手順例を示すフローチャート（その１）である。 図２０は、システムによる構造解析処理手順例を示すフローチャート（その２）である。
　以下に添付図面を参照して、本発明にかかる構造解析方法、及び構造解析プログラムの実施の形態を詳細に説明する。
　図１は、本発明にかかる構造解析装置による一の動作例を示す説明図である。構造解析装置１００は、解析対象のＦＥＭ（Ｆｉｎｉｔｅ　Ｅｌｅｍｅｎｔ　Ｍｅｔｈｏｄ）モデルについて有限要素法を用いた構造解析を行うコンピュータである。ここで、有限要素法とは、微分方程式を近似的に解くための数値解析の方法である。構造解析による評価では、例えば、剛性の指標として変位、強度の指標として応力などが評価される。
　解析対象のＦＥＭモデルは、以降単にモデル、解析対象モデルなどとも称する。有限要素法では、方程式が定義された領域を小領域に分割し、各小領域における方程式を共通な補間関数で近似する方法である。ここで、小領域を要素１０２と称する。要素１０２は、例えば、各頂点などに節点１０３が配置される。有限要素法を用いた構造解析では、方程式が定義された領域が解析対象の構造物である。解析対象モデル１０１は、解析対象の構造物を近似させることによりＦＥＭによる構造解析が可能な状態にしたものであり、シミュレーション空間上に配置されるものである。
　シミュレーション空間とは、解析対象の物理的な構造物の内部とその構造物の外部の空間とを含む解析領域をコンピュータ上に表すために設定された空間である。解析対象モデル１０１では、解析対象の構造物を細かな小領域である要素１０２に分割し、要素１０２の集まりとして全体の形状を表す。要素１０２の集まりをメッシュとも称する。メッシュにおける１つの格子のサイズをメッシュサイズとも称する。
　１つの要素１０２の形状は、２次元の場合に四角形や三角形などであり、３次元の場合に四面体や五面体や六面体などである。図１では、要素１０２として四面体の二次要素を例に挙げるが、要素１０２の形状の例や要素１０２と節点１０３との関係については図２を用いて詳細に説明する。
　開発者は、例えば、ＣＡＤ（Ｃｏｍｐｕｔｅｒ　Ａｉｄｅｄ　Ｄｅｓｉｇｎ）を用いて構造物をポリゴンで表した３次元形状の立体モデルの立体モデル情報を作成する。立体モデル情報は例えばポリゴンの座標データなどを有する。つぎに、開発者は、ＣＡＥ（Ｃｏｍｐｕｔｅｒ　Ａｉｄｅｄ　Ｅｎｇｉｎｅｅｒｉｎｇ）を用いて、設計した構造物が要求性能を満たすかどうかを、実際に構造物を作る前にコンピュータ上でシミュレーションする。ここで、ＣＡＥは、例えば、プリプロセッサとソルバーとポストプロセッサと呼ばれる３つのソフトウェアを有する。
　まず、ＣＡＥのプリプロセッサは、メッシュ分割や境界条件定義を行うことによりＦＥＭモデルを表すモデル情報１１１を作成することができる。つぎに、ＣＡＥのソルバーは、モデル情報１１１が表すＦＥＭモデルに対して数値解析を行う構造解析ソルバーである。最後に、ＣＡＥのポストプロセッサは、ＣＡＥのソルバーによって得られた解析結果の可視化を行う。
　従来、ＦＥＭによる構造解析では、力のつりあい方程式を解くことにより、構造物の変位や応力を求めることができる。力のつりあい方程式は連立一次方程式である。連立一次方程式は「Ｋｕ＝ｆ」とも表す。Ｋは剛性行列である。ｕはそれぞれの節点１０３の自由度ごとの変位を表すベクトルである。ｕは解析で求める未知ベクトルである。ｆは荷重を表すベクトルである。ｆは、荷重条件を設定した場合、各節点１０３の自由度に対応した荷重値である。
　従来、連立一次方程式の求解は、計算量が多く計算に時間がかかるため、計算機を用いてより高速に解くための様々な手法が開発されている。連立一次方程式の解法として、直接法や反復法がある。直接法は、例えば、係数行列を三角化することで解を直接求める解法である。直接法としては、例えば、マルチフロンタル法などがある。反復法は、反復計算により近似解を求める解法である。反復法としては、例えば、前処理付きＣＧ法などがある。
　しかしながら、従来技術では、解析対象モデル１０１の有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを特定の解法として数値解析を行うと、モデルによっては他の解法よりも解析時間が長くなるという問題点がある。例えば、一般的に直接法は反復法よりメモリ量が多く必要である。また、直接法は連立一次方程式の特性により、非ゼロ領域が増えると演算量やメモリ量がより多く必要となり、解析時間が長くなる場合がある。一方、例えば、反復法は、拘束条件によるモデル１０１への拘束が不十分な場合などモデルが不安定であると、解が収束しにくいことがあり、解析時間が長くなる場合がある。
　そこで、本実施の形態では、構造解析装置１００は、解析対象モデル１０１を表す情報に基づきモデル規模を評価し、モデル規模に応じてＦＥＭを用いた連立一次方程式による解法を直接法と反復法とから選択する。これにより、メモリの使用量を推定できるため、メモリを不足させずに直接法による解析を行い、解析時間の短縮化を図る。
　まず、構造解析装置１００は、モデル情報１１１を取得する。モデル情報１１１は、例えば、解析対象モデル１０１の形状を表す情報を含む。モデル情報１１１の詳細な例は、後述する図６を用いて説明する。
　つぎに、構造解析装置１００は、取得したモデル情報１１１に基づき、解析対象モデル１０１の規模を評価する。構造解析装置１００は、例えば、モデル情報１１１とメッシュサイズと要素タイプとから算出された要素数及び／又は節点数及び／又は解析自由度数に基づいてモデル規模の評価を行う。
　そして、構造解析装置１００は、評価したモデル規模に基づいて、モデル１０１の有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する。構造解析装置１００は、例えば、モデル規模が大きいほど、反復法を選択する。
　解析対象モデル１０１のモデル規模が大きいと、解析する節点１０３の数が多いため、メモリ使用量が多くなる。このため、構造解析装置１００は、評価したモデル規模と、解析時に使用可能なメモリ使用量と、に基づいて、アルゴリズムを直接法とするか反復法とするか選択する。
　具体的に、構造解析装置１００は、例えば、直接法を用いてモデル１０１のモデル規模を解析する場合に使用するメモリ量が、解析時に使用可能なメモリ使用量未満である場合に、直接法を選択する。一方、構造解析装置１００は、例えば、直接法を用いてモデル１０１のモデル規模を解析する場合に使用するメモリ量が、解析時に使用可能なメモリ使用量以上である場合に、反復法を選択する。図１の例では、反復法が選択される。
　例えば、モデルを解析する場合に使用するメモリ量は、要素タイプや要素数や節点数などから算出されるモデル規模に基づいて推測値を算出可能である。例えば、直接法を用いる場合に要素１つに対してどの程度のメモリ量を使用するかを予め記憶装置などに用意しておいてもよい。そして、構造解析装置１００は、解析対象モデル１０１の要素数に応じて使用するメモリ量を算出してもよい。
　例えば、一般的に直接法による解析の方が反復法による解析よりも精度が高いため、直接法による解析に時間がかからないことが分かる場合に、精度の向上のために、直接法を利用したい場合がある。しかし、上述したように、直接法による解析の場合、メモリの使用量が不足すると、解析に時間がかかる。このため、モデル規模に応じて解析法を選択することで、モデル規模によってメモリの使用量を推定できるため、メモリを不足させずに直接法による解析を行い、精度の向上を図りつつ、解析時間の短縮化を図ることができる。
　図２は、本発明にかかる構造解析装置による他の動作例を示す説明図である。構造解析装置１００は、解析対象のＦＥＭモデルについて有限要素法を用いた構造解析を行うコンピュータである。図２に示す構造解析装置１００と、図１に示す構造解析装置１００とは同一装置であっても、異なる装置であってもよいため、同一の符号を付す。また、ここで、図２に示す例と、図１に示す例とにおいて、同一の構成や機能についての詳細な説明を省略する。
　図１で説明したように、従来技術では、解析対象モデル１０１の有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを特定の解法として数値解析を行うと、モデルによっては他の解法よりも解析時間が長くなるという問題点がある。直接法は、例えば、解析時にメモリ量を多く使用するため、解析時間が長くなる場合がある。一方、反復法は、例えば、拘束条件によるモデル１０１への拘束が不十分であると、解が収束しない場合がある。このため、解析時間が長くなる場合がある。
　そこで、本実施の形態では、構造解析装置１００は、モデル１０１の境界条件の強さを評価し、境界条件の強さに応じて、ＦＥＭを用いた連立一次方程式による解法を直接法と反復法とから選択する。これにより、解析時間の短縮化を図ることができる。ここで、境界条件とは、構造物の負荷状態を表すために境界に設定する条件である。境界条件には、例えば、拘束条件と荷重条件の２種類がある。拘束条件は、解析対象の構造物に対して、どの部位をどのように固定するか、どのくらい強制的な変位を与えるか、などの条件のことである。荷重条件はどの部位にどのような方向でどのくらい力を加えるかなどの条件である。境界条件の強さが、例えば、拘束条件の強さの場合、拘束条件による境界の拘束により計算モデルが安定しているか否かを判定でき、解析時間の短縮化を図ることができる。
　まず、構造解析装置１００は、解析対象モデル１０１の境界条件を取得する。つぎに、構造解析装置１００は、取得した境界条件に基づき、解析対象モデル１０１の境界条件の強さを評価する。境界条件の強さは、例えば、拘束条件の強さと荷重条件の弱さとの少なくともいずれかを含む。
　拘束条件は、拘束の種類や方向、拘束する位置や領域、などが設定される。拘束条件の強さは、例えば、拘束条件に応じて解析対象モデル１０１が拘束される強さである。拘束条件の強さは、例えば、後述するように拘束条件によって解析対象モデル１０１の表面のうち拘束される範囲に相当する領域の面積に基づいて評価される。荷重条件は、荷重の種類、荷重の大きさや方向、荷重をかける位置や領域、などが設定される。荷重条件の弱さは、例えば、荷重条件が剛性に与える影響の大きさの逆数である。荷重条件の弱さは、例えば、荷重条件が剛性に与える影響の大きさの逆数に基づいて評価される。
　構造解析装置１００は、例えば、評価した境界条件の強さに基づいて、モデル１０１の有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する。構造解析装置１００は、例えば、境界条件の強さが大きいほど、反復法を選択する。
　ここでは、拘束条件の強さを例に挙げて説明する。解析対象モデル１０１の拘束条件の強さが大きいと、拘束条件による境界の拘束により解析モデルが安定するため、反復法において解が収束する可能性が高くなる。一方、解析対象モデル１０１の拘束条件の強さが小さいと、拘束条件による境界の拘束が不十分な可能性があり、反復法において解が収束しない可能性が高くなる。
　そこで、構造解析装置１００は、例えば、拘束条件の強さが閾値以上の場合に、反復法を選択し、拘束条件の強さが閾値未満の場合に、直接法を選択する。このように、拘束条件の強さによって境界の安定度合いを判別でき、解析時間の短縮化を図ることができる。
　図３は、有限要素法の構造解析に用いる要素と節点を示す説明図である。図３（１）に示すように、構造物は、例えば、連続体である。上述したように、ＣＡＥのプリプロセッサは、３次元のシミュレーション空間上に作成された構造物の立体モデルを近似することによりＦＥＭモデルを作成する。
　図３（２）に示すように、要素タイプは、３次元のソリッド要素、２次元のシェル要素、１次元のビーム要素やトラス要素、０次元の質量要素などが挙げられる。ソリッド要素の形状は、六面体、五面体、四面体などがある。シェル要素の形状は、四角形、三角形などがある。また、ビーム要素やトラス要素の形状は、線である。質量要素の形状は、点である。
　図３（３）に示すように、例えば、要素形状が四面体の場合、一次要素と二次要素のように節点の定義方法が複数ある。一次要素の場合、要素に含まれる節点は、四面体の頂点である。二次要素の場合、節点は、四面体の頂点と、頂点間を結ぶ線の中心点と、である。このように、要素形状や節点の定義方法により要素を構成する節点数が異なる。本実施の形態では、要素タイプとして四面体の二次要素を用いて説明する。
　ここで、ソリッド要素は、固体要素とも称する。また、シェル要素は、例えば、板要素とも称する。また、ビーム要素は、梁要素とも称する。これらの要素タイプにより、それぞれ後述する節点の解析自由度数が異なる。
（構造解析装置１００のハードウェア構成例）
　図４は、構造解析装置のハードウェア構成例を示すブロック図である。図４において、構造解析装置１００は、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）４０１と、ＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）４０２と、ＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）４０３と、ディスクドライブ４０４と、ディスク４０５と、を有する。また、構造解析装置１００は、Ｉ／Ｆ（Ｉｎｔｅｒｆａｃｅ）４０６と、キーボード４０７と、マウス４０８と、ディスプレイ４０９と、を有する。また、ＣＰＵ４０１と、ＲＯＭ４０２と、ＲＡＭ４０３と、ディスクドライブ４０４と、Ｉ／Ｆ４０６と、キーボード４０７と、マウス４０８と、ディスプレイ４０９とは、バス４００によってそれぞれ接続される。
　ここで、ＣＰＵ４０１は、構造解析装置１００の全体の制御を司る。ＲＯＭ４０２は、ブートプログラムなどのプログラムを記憶する。ＲＡＭ４０３は、ＣＰＵ４０１のワークエリアとして使用される。ディスクドライブ４０４は、ＣＰＵ４０１の制御にしたがってディスク４０５に対するデータのリード／ライトを制御する。ディスク４０５は、ディスクドライブ４０４の制御で書き込まれたデータを記憶する。ディスク４０５としては、磁気ディスク、光ディスクなどが挙げられる。
　Ｉ／Ｆ４０６は、通信回線を通じてＬＡＮ（Ｌｏｃａｌ　Ａｒｅａ　Ｎｅｔｗｏｒｋ）、ＷＡＮ（Ｗｉｄｅ　Ａｒｅａ　Ｎｅｔｗｏｒｋ）、インターネットなどのネットワーク４１０に接続され、このネットワーク４１０を介して他の装置に接続される。そして、Ｉ／Ｆ４０６は、ネットワーク４１０と内部のインターフェースを司り、外部装置からのデータの入出力を制御する。Ｉ／Ｆ４０６には、例えばモデムやＬＡＮアダプタなどを採用することができる。
　キーボード４０７やマウス４０８は、利用者の操作により、各種データの入力を行うインターフェースである。ディスプレイ４０９は、ＣＰＵ４０１の指示により、データを出力するインターフェースである。
　また、図示を省略するが、構造解析装置１００には、カメラから画像や動画を取り込む入力装置やマイクから音声を取り込む入力装置が設けられていてもよい。また、図示を省略するが、構造解析装置１００には、プリンタなどの出力装置が設けられていてもよい。
（構造解析装置１００の機能的構成例）
　図５は、構造解析装置の機能的構成例を示すブロック図である。構造解析装置１００は、取得部５０１と、第１評価部５０２と、第２評価部５０３と、選択部５０４と、を有する。取得部５０１から選択部５０４までの制御部の処理は、例えば、図４に示すＣＰＵ４０１がアクセス可能なＲＯＭ４０２、ＲＡＭ４０３、ディスク４０５などの記憶装置に記憶されたプログラムにコーディングされている。そして、ＣＰＵ４０１が記憶装置から該プログラムを読み出して、プログラムにコーディングされている処理を実行する。これにより、制御部の処理が実現される。また、制御部の処理結果は、例えば、ＲＯＭ４０２、ＲＡＭ４０３、ディスク４０５などの記憶装置に記憶される。
　取得部５０１は、モデル情報１１１を取得する。モデル情報１１１は、解析対象モデルを表す情報である。モデル情報１１１は、例えば、解析対象モデルの形状を表す情報を含む。また、モデルの形状を表す情報は、モデルに含まれる要素の要素情報、要素に含まれる節点及び節点の位置情報などが挙げられる。解析対象モデルの形状を表す情報は、例えば、体積、表面積、などの情報であってもよい。また、取得されるモデル情報１１１は、例えば、モデルの材質を表す情報を含んでいてもよい。また、取得されるモデル情報１１１は、例えば、解析時にモデルに設定される境界条件の情報を含んでいてもよい。また、取得されるモデル情報１１１は、例えば、結果の出力項目の情報を含んでいてもよい。結果の出力項目の情報は、例えば変位、応力、反力など、どの解析結果を計算出力するかの情報である。
　図６は、モデル情報例を示す説明図である。図６（１）～図６（３）には、モデル情報１１１として、モデルに含まれる要素の要素情報、要素に含まれる節点及び節点の位置情報などを示す。また、図６（４）には、モデル情報１１１として、モデルの形状を表す情報やモデルの材質を表す情報や解析時にモデルに設定される境界条件を表す情報などを示す。ここで、境界条件とは、構造物の負荷状態を表すために境界に設定する条件である。境界条件には、例えば、拘束条件と荷重条件の２種類がある。拘束条件は、解析対象の構造物に対して、どの部位をどのように固定するか、どのくらい強制的な変位を与えるか、などの条件のことである。荷重条件はどの部位にどのような方向でどのくらい力を加えるかなどの条件である。
　図６（１）では、モデル情報１１１は、モデルＩＤ、要素ＩＤのフィールドを有する。モデルＩＤのフィールドには、モデルを識別する識別情報が設定される。要素ＩＤのフィールドには、解析対象のモデルに含まれる要素を識別する識別情報が設定される。図６（１）の例では、モデルｍａは、例えば、要素ｅｎ１～ｅｎ３を有する。
　つぎに、図６（２）の例では、モデル情報１１１は、要素ＩＤ、節点ＩＤのフィールドを有する。要素ＩＤのフィールドには、要素を識別する識別情報が設定される。節点ＩＤのフィールドには、要素に含まれる節点を識別する識別情報が設定される。図６（２）の例では、要素ｅｎ１は、節点ｎｎ１～ｎｎ４を有する。
　そして、図６（３）の例では、モデル情報１１１は、節点ＩＤ、節点位置情報のフィールドを有する。節点ＩＤのフィールドには、節点を識別する識別情報が設定される。節点位置情報のフィールドには、節点の位置情報が設定される。節点の位置情報は、例えば、節点のｘ軸、ｙ軸、ｚ軸の座標値である。
　図６（１）～（３）に示すモデル情報１１１によって、構造解析装置１００は、モデルの形状を特定可能である。例えば、構造解析装置１００は、モデルの体積、表面積などを算出可能である。
　つぎに、図６（４）の例では、モデル情報１１１は、モデルＩＤ、体積、質量、主慣性モーメント、節点数、要素タイプ、メッシュサイズ、要素数、解析自由度数、最小の球の体積、表面積、境界条件などのフィールドを有する。モデルＩＤのフィールドには、モデルを識別する識別情報が設定される。体積のフィールドには、モデルの体積が設定される。質量のフィールドには、モデルの質量が設定される。主慣性モーメントのフィールドには、モデルの主慣性モーメントが設定される。ここで、体積や質量や主慣性モーメントは、モデルの材質や形状から決定される情報である。
　要素タイプのフィールドには、図２に示したような要素タイプを識別する情報が設定される。また、解析自由度数のフィールドには、例えば、モデルに含まれる各節点の自由度数の総和が設定される。最終的な解析自由度数は、ＣＡＥのソルバー内部で行われる連立一次方程式の組立が完了しないと判明しない。このため、解析自由度数は、例えば、節点数×１節点あたりの自由度数などによる概算値が設定されてもよい。
　メッシュサイズのフィールドには、メッシュ分割する際の代表的な１要素の要素長が設定される。本実施の形態では、すべての要素が同一の要素タイプである場合を例に挙げて説明するため、メッシュサイズは一通りである。
　節点数のフィールドには、モデルに含まれる節点の数が設定される。また、要素数のフィールドには、モデルに含まれる要素の数が設定される。要素数については、メッシュサイズと、形状の体積と、要素タイプと、によって予測値を算出することが可能である。また、節点数は、要素数と、要素タイプとに基づいておおよその予測値を算出することが可能であり、さらにモデル形状を考慮することで予測の精度を高めることも可能である。また、要素数と節点数とについては、図６（１）～（３）に示すモデル情報１１１がすでにある場合、当該モデル情報１１１に基づいて特定可能である。
　最小の球の体積のフィールドには、モデルを覆う最小の球の体積が設定される。ここでは、最小の球の体積のフィールドがあるが、これに限らず、例えば、モデルを覆う最小の立方体の体積のフィールドがあってもよい。
　表面積のフィールドには、モデルの表面積が設定される。境界条件のフィールドは、例えば、拘束条件、拘束面積、拘束自由度数、荷重条件のフィールドを有する。そして、拘束条件のフィールドには、例えば、拘束条件が設定される。拘束条件としては、例えば、片側端面固定、両端固定、全周固定などが挙げられる。拘束面積のフィールドには、例えば、モデルの表面積のうち、拘束条件に応じて拘束される面積が設定される。拘束自由度数のフィールドには、例えば、１節点当たりの自由度数のうち拘束される自由度数が設定される。例えば、１節点当たりの自由度が３方向である場合に、３方向のうち何方向を拘束するかが拘束自由度数である。荷重条件のフィールドには、荷重条件が設定される。荷重条件としては、例えば、各方向における重力荷重などが挙げられる。
＜モデル規模に応じた選択＞
　まず、構造解析装置１００は、解析対象モデル１０１を表す情報に基づきモデル規模を評価し、モデル規模に応じてＦＥＭを用いた連立一次方程式による解法を直接法と反復法とから選択する例を説明する。
　第１評価部５０２は、取得部５０１によって取得されたモデル情報１１１に基づいて、モデルのモデル規模を評価する。モデルのモデル規模は、例えば、要素数、節点数、解析自由度数であってもよい。また、モデル規模は、係数行列に関する、バンド幅、スカイライン、プロフィル、非ゼロ項数、フィルインの数などであってもよく、これらを組み合わせて算出した値であってもよい。
　具体的に、第１評価部５０２は、例えば、モデル情報１１１とメッシュサイズとから算出された要素数及び／又は節点数及び／又は解析自由度数に基づいてモデル規模の評価を行う。例えば、モデル情報１１１が、例えば、モデルの体積などのモデルの形状を表す情報を含む場合、第１評価部５０２は、モデルの体積と、メッシュサイズと、要素タイプと、に基づいて要素数を算出する。ここで算出される要素数は予測値である。また、第１評価部５０２は、例えば、モデルの体積と、メッシュサイズと、要素タイプと、に基づいて節点数を算出する。ここで算出される節点数は、予測値となる。また、第１評価部５０２は、例えば、モデルの体積と、メッシュサイズと、要素タイプと、拘束条件と、に基づいて解析自由度数を算出する。ここで算出される解析自由度数は、予測値となる。
　また、モデル情報１１１にモデルに含まれる要素や節点や解析自由度数の情報などがモデル情報１１１に含まれる場合、第１評価部５０２はモデル情報１１１から要素数や節点数や解析自由度数を特定してもよい。
　また、上述したように、解析自由度数は、例えば、モデルに含まれる各節点の自由度数の総和である。解析自由度数は、ＣＡＥのソルバー内部で行われる連立一次方程式の組立が完了しないと判明しない。このため、解析自由度数は、節点数×１節点あたりの自由度数によって予測値を算出することが可能である。ここでは、モデル規模の指標値として要素数を例に挙げて説明する。１節点あたりの自由度数は要素タイプにより決定されるが、例えば、モデル情報１１１に含まれていてもよい。
　つぎに、選択部５０４は、モデル規模の評価結果に基づいて、モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する。選択部５０４は、例えば、モデル規模が大きいほど、反復法が選択されるようにする。
　具体的に、選択部５０４は、例えば、空きメモリ量を取得する。空きメモリ量とは、例えば、ＲＡＭ４０３やディスク４０５などの記憶装置のうちの空き領域の量である。選択部５０４は、空きメモリ量と、評価したモデル規模と、に基づいて、直接法とするか反復法とするか選択する。モデル規模が大きいほど、メモリの使用量が多くなる。
　例えば、モデル規模に応じて空きメモリ量についての閾値を用意しておき、選択部５０４は、取得した空きメモリ量が、モデル規模に応じた閾値以上であれば、直接法を用いる。選択部５０４は、取得した空きメモリ量が、モデル規模に応じた閾値未満であれば、反復法を用いる。
　図７は、要素数に応じたメモリ容量の閾値例を示す説明図である。閾値テーブル７００は、例えば、要素数ごとのメモリ量の閾値を有する。要素数は、例えば、１～３０００、３００１～５０００、５００１～７０００、７００１～１００００である。例えば、モデル規模が４５００であれば、３００１～５０００に対応する閾値ｍｔ２が選ばれる。
　例えば、モデル規模が大きいほど、メモリ量についての閾値は、大きくなる。例えば、閾値ｍｔ１＜閾値ｍｔ２＜閾値ｍｔ３＜閾値ｍｔ４である。これにより、モデル規模が大きいほど、閾値が大きいため、反復法を用いることになる。このため、解析時のメモリ量の不足を回避でき、解析時間の短縮化を図ることができる。また、モデル規模が大きくとも、空きメモリ量が多ければ、直接法を用いることができるため、解析精度の向上を図ることができる。
（構造解析装置１００によるモデル規模に基づく選択処理手順例）
　図８は、構造解析装置によるモデル規模に基づく選択処理手順例を示すフローチャートである。構造解析装置１００は、解析対象モデルを表すモデル情報１１１を取得する（ステップＳ８０１）。そして、構造解析装置１００は、モデル規模を評価する（ステップＳ８０２）。
　つぎに、構造解析装置１００は、空きメモリ量を確認する（ステップＳ８０３）。構造解析装置１００は、モデル規模と空きメモリ量に応じて解法を選択する（ステップＳ８０４）。そして、構造解析装置１００は、選択した解法により数値解析を実行し（ステップＳ８０５）、一連の処理を終了する。
＜境界条件の強さとモデル規模に応じた選択＞
　つぎに、境界条件の強さに基づいて解析法を選択する例と、境界条件の強さとモデル規模とに基づいて解析法を選択する例と、について説明する。
　第２評価部５０３は、取得したモデル情報１１１と境界条件とに基づいて、境界条件の強さを評価する。図５に示すモデル情報１１１において説明したように、境界条件には、拘束条件と、荷重条件とがある。境界条件の強さは、例えば、拘束条件の強さと荷重条件の弱さとの少なくともいずれかを含む。
　拘束条件は、例えば、拘束の種類や方向、拘束する位置や領域等が設定される。拘束条件の強さは、例えば、拘束条件に応じて解析対象モデルが拘束される強さである。拘束条件の強さは、例えば、後述するように拘束条件によって解析対象モデルの表面のうち拘束される範囲に相当する領域の面積に基づいて評価される。荷重条件は、例えば、荷重の種類、荷重の大きさや方向、荷重をかける位置や領域、が設定される。荷重条件の弱さは、例えば、荷重条件が剛性に与える影響の大きさの逆数に基づいて評価される。
　ここでは、境界条件に含まれる拘束条件を例に挙げ、境界条件の強さとして拘束条件の強さを例に挙げて説明する。第２評価部５０３は、モデルの拘束面積、拘束自由度に基づいて、拘束条件の強さを評価する。拘束面積は、拘束条件に応じてモデル内で拘束される部分の面積である。より具体的に、第２評価部５０３は、以下式（１）によって拘束条件の強さを算出する。
　拘束条件の強さ＝拘束面積／表面積×（１節点当たりの拘束自由度数／１節点当たりの自由度数）・・・式（１）
　モデルに複数の拘束条件が設定されている場合は、これらの総和を、モデル全体の拘束条件の強さとしてもよい。
　選択部５０４は、拘束条件の強さの評価結果に基づいて、モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する。選択部５０４は、例えば、拘束条件の強さが大きいほど、反復法が選択されるようにする。具体的に、選択部５０４は、例えば、拘束条件の強さが閾値以上である場合に、反復法を選択し、拘束条件の強さが閾値未満である場合に、直接法を選択する。ここで、閾値は、予め設定されてＲＯＭ４０２、ＲＡＭ４０３、ディスク４０５などの記憶装置に記憶される。
　つぎに、選択部５０４による詳細な選択例の説明の前に、選択部５０４による選択時の閾値の設定方法について説明する。ここでは、閾値の設定方法を構造解析装置１００が行うとして説明するが、これに限らず、構造解析装置１００と異なる装置が閾値を決定しておき、構造解析装置１００に閾値を提供してもよい。
　まず、構造解析装置１００は、例えば、モデルに対して複数の拘束条件を用意し、拘束条件の各々について拘束条件の強さを算出する。また、構造解析装置１００は、複数の拘束条件の各々について、モデルに対して、直接法及び反復法による数値解析を実行する。構造解析装置１００は、直接法による解析時間と、反復法による解析時間との性能比が１となる場合の拘束条件の強さを線形補間により求めて閾値とする。
　図９は、閾値の導出に用いるモデル例を示す説明図である。モデルｍ１を例に挙げて閾値の導出について説明する。モデルｍ１の寸法は、例えば、５０×１００×１０００［ｍｍ］の直方体である。例えば、拘束条件１は、片持ちである。例えば、拘束条件２は、両端固定である。例えば、拘束条件３は、全周固定である。
　本実施の形態では、モデルｍ１を要素数が５０００程度になるようにメッシュ分割した場合と、要素数が１００００程度となるようにメッシュ分割した場合と、の各々について、拘束条件の強さと解析時間とによって閾値を設定する例を示す。
　図１０は、要素数が５０００の場合の閾値の導出例を示す説明図である。図１０において（１）～（４）には拘束条件の強さの導出に用いる各数値と、（５）と（６）には解析時間を示す。
　構造解析装置１００は、例えば、上述した式（１）に示したように、「拘束面積／表面積×（１節点当たりの拘束自由度数／１節点当たりの自由度数）」によって拘束条件の強さを算出する。拘束面積は（３）に示す。表面積は（１）に示す。１節点当たりの拘束自由度数は（４）に示す。１節点当たりの自由度数は（２）に示す。
　拘束条件１の片持ちの場合、拘束条件の強さは、「（５０００／３１００００）×（３／３）」により「０．０１６」となる。
　拘束条件２の両端固定の場合、拘束条件の強さは、「（１００００／３１００００）×（３／３）」により「０．０３２」となる。
　拘束条件３の全周固定の場合、拘束条件の強さは、「（１１００００／３１００００）×（３／３）」により「０．３５５」となる。
　構造解析装置１００は、例えば、拘束条件の各々について、直接法による数値解析と、反復法による数値解析と、を行うことにより解析時間を求める。そして、構造解析装置１００は、直接法による解析時間と反復法による解析時間との比を解析方法の性能比として求める。直接法による解析時間は（５）に示す。反復法による解析時間は（６）に示す。
　図１０の例では、拘束条件２の性能比と拘束条件３の性能比との間に性能比１がある。そこで、構造解析装置１００は、拘束条件２の拘束条件の強さと、拘束条件３の拘束条件の強さと、を線形補間して性能比が１となる場合の拘束条件の強さを求める。そして、構造解析装置１００は、性能比が１となる場合の拘束条件の強さを閾値として設定する。例えば、線形補間に限らず、ラグランジュ補間、スプライン補間、最小二乗法によって求めてもよい。
　閾値＝拘束条件２の強さ＋（拘束条件３の強さ－拘束条件２の強さ）×（１－拘束条件２の性能比）／（拘束条件３の性能比－拘束条件２の性能比）
　＝０．０３２＋（０．３５５－０．０３２）×（１－０．３７）／（１．５７－０．３７）
　＝０．２０２
　拘束条件２の場合、反復法による解析時間よりも直接法による解析時間が短い。これに対して、拘束条件３の場合、直接法による解析時間よりも反復法による解析時間が短い。このため、上述したように、選択部５０４は、第２評価部５０３により算出された拘束条件の強さが閾値以上であれば、反復法を選択し、閾値未満であれば、直接法を選択することにより、解析時間の短縮化を図ることができる。
　図１１は、要素数が１００００の場合の閾値の導出例を示す説明図である。図１１において（１）～（４）には拘束条件の強さの導出に用いる各数値と、（５）と（６）には解析時間を示す。拘束条件の強さについては、図１０の例と同じである。図１１の例では性能比が１となるのは、拘束条件２の性能比と拘束条件３の性能比との間である。このため、図１１では、図１０と同様に閾値を求める例であるため、詳細な説明を省略するが、以下のように求まる。
　閾値＝拘束条件２の強さ＋（拘束条件３の強さ－拘束条件２の強さ）×（１－拘束条件２の性能比）／（拘束条件３の性能比－拘束条件２の性能比）
　＝０．０３２＋（０．３５５－０．０３２）×（１－０．４７）／（２．１６－０．４７）
　＝０．１３３
　また、図５に示す選択部５０４は、モデル規模の評価結果及び拘束条件の強さの評価結果に基づいて、モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する。モデル規模の評価については、上述した第１評価部５０２による評価であるため、詳細な説明を省略する。
　具体的に、選択部５０４は、例えば、拘束条件の強さに基づいて、解法を直接法とするか反復法とするか選択する際に、モデル規模に基づいて選択基準を変更する。また、具体的に、選択部５０４は、例えば、モデル規模に基づいて解法を直接法とするか反復法とするかを選択し、拘束条件の強さに基づいて解法を直接法とするか反復法とするかを選択する。そして、選択部５０４は、例えば、モデル規模に基づく選択結果と、拘束条件の強さに基づく選択結果と、に基づいて、最終的に解法を直接法とするか反復法とするかを選択してもよい。ここでは、モデル規模に基づいて選択基準を変更する例について説明する。
　図１２は、モデル規模に応じた閾値例を示す説明図である。閾値テーブル１２００は、例えば、要素数ごとに閾値を有する。図１２では、モデル規模の評価の一例として、要素数を挙げるが、これに限らず、上述したように節点数や解析自由度数などであってもよい。
　閾値テーブル１２００では、要素数は、例えば、３０００、５０００、７５００、１００００である。要素数が３０００に対応する閾値ｔｈ１＞要素数が５０００に対応する閾値ｔｈ２＞要素数が７５００に対応する閾値ｔｈ３＞要素数が１００００に対応する閾値ｔｈ４である。このように、要素数が小さいほど、閾値が大きくなり、要素数が大きいほど、閾値が小さくなる。これにより、要素数が大きいほど、閾値が小さいため、要素数が大きいほど、反復法が選択されるようになる。
　選択部５０４は、閾値テーブル１２００から要素数に応じた閾値を取得して、算出した拘束条件の強さと取得した閾値とを比較する。具体的に、選択部５０４は、例えば、解析対象モデルの要素数が４０００以下であれば、要素数が３０００に対応する閾値を取得する。具体的に、選択部５０４は、例えば、解析対象モデルの要素数が４０００より大きく６２５０以下であれば、要素数が５０００に対応する閾値を取得する。このように、要素数の範囲を設けてもよい。そして、選択部５０４は、算出した拘束条件の強さが取得した閾値以上である場合に反復法を選択し、算出した拘束条件の強さが取得した閾値未満である場合、直接法を選択する。
　つぎに、具体例を用いて拘束条件の強さの算出例と拘束条件の強さによる解法の選択例について説明する。
　図１３は、拘束条件の強さの算出例と解法の選択例を示す説明図である。解析対象モデルは、例えば、寸法が５０×１００×５００である。
　第１評価部５０２は、例えば、解析対象モデルのモデル規模を評価する。ここでは、モデル規模の指標値として要素数を用いる。つぎに、第２評価部５０３は、例えば、解析対象モデルの拘束条件に基づいて拘束条件の強さを評価する。
　拘束条件の強さ＝拘束面積／表面積×（１節点当たりの拘束自由度数／１節点当たりの自由度数）
　＝（１０００／１６０００）×（３／３）
　＝０．０６２５
　選択部５０４は、例えば、閾値テーブル１２００から評価したモデル規模に応じて閾値を取得する。図１３の例では、要素数が４７２４であるため、選択部５０４は、例えば、閾値テーブル１２００から、要素数が５０００に対応する閾値を取得する。ここで、要素数が５０００に対応する閾値として、図１０の例を用い、「０．２０２」である。
　このため、選択部５０４は、拘束条件の強さが閾値未満であると判定し、直接法を選択する。拘束条件によるモデルの拘束が不十分な場合に、反復法を用いると、解がなかなか収束しない。そのため、拘束条件の強さが閾値未満のようにモデルの拘束条件の強さが弱い場合に、直接法が選択されることにより、解析時間の短縮化を図ることができる。
（構造解析装置１００による拘束条件の強さに基づく選択処理手順例）
　図１４は、構造解析装置による拘束条件の強さに基づく選択手順例を示すフローチャートである。構造解析装置１００は、解析対象モデルのモデル情報１１１を取得する（ステップＳ１４０１）。
　構造解析装置１００は、境界条件の強さを評価する（ステップＳ１４０２）。そして、構造解析装置１００は、モデル規模を評価する（ステップＳ１４０３）。つぎに、構造解析装置１００は、モデル規模に応じた閾値を取得する（ステップＳ１４０４）。構造解析装置１００は、境界条件の強さと閾値とを比較する（ステップＳ１４０５）。
　境界条件の強さが閾値未満である場合（ステップＳ１４０５：閾値未満）、構造解析装置１００は、直接法を選択し（ステップＳ１４０６）、ステップＳ１４０８へ移行する。一方、境界条件の強さが閾値以上である場合（ステップＳ１４０５：閾値以上）、構造解析装置１００は、反復法を選択し（ステップＳ１４０７）、ステップＳ１４０８へ移行する。
　そして、構造解析装置１００は、選択した解法により数値解析を実行し（ステップＳ１４０８）、一連の処理を終了する。
　図１５は、システム適用例を示す説明図である。システム１５００は、例えば、クライアント端末装置１５０１と、データセンタ１５０２と、を有する。クライアント端末装置１５０１は、図１５に示すように複数であってもよい。データセンタ１５０２は、例えば、複数のサーバ１５０３を有する。クライアント端末装置１５０１とデータセンタ１５０２とはネットワーク４１０を介して接続される。データセンタ１５０２におけるサーバ１５０３の管理方法については特に限定しない。
　構造解析装置１００は、クライアント端末装置１５０１と、複数のサーバ１５０３と、によって実現されてもよい。例えば、取得部５０１と、第１評価部５０２と第２評価部５０３についてはクライアント端末装置１５０１によって実現され、選択部５０４についてはサーバ１５０３によって実現されてもよい。ここで、データセンタ１５０２のいずれのサーバ１５０３が選択部５０４による処理を行うかを予め決定していてもよい。
　サーバ１５０３は、例えば、ＣＡＤとＣＡＥを実行する。そして、クライアント端末装置１５０１は、例えば、ＣＡＥによって作成されたモデル情報１１１を取得して拘束条件の強さを評価する。そして、クライアント端末装置１５０１は、評価した拘束条件の強さをサーバ１５０３へ送信する。そして、サーバ１５０３は、拘束条件の強さに基づいて反復法と直接法とのいずれかの解法を選択し、選択した解法により数値解析を実行する。ここで、サーバ１５０３は、過去の解析の実績をデータベースに登録してもよい。そして、サーバ１５０３は、クライアント端末装置１５０１から受信した拘束条件の強さに基づいて解法を選択してもよい。
　図１６は、解析データベース例を示す説明図である。解析データベース１６００は、例えば、モデルＩＤ、境界条件の強さ、モデル規模、直接法Ａの解析時間、直接法Ｂの解析時間、反復法Ｃの解析時間、反復法Ｄの解析時間のフィールドを有する。各フィールドに情報が設定されることによりレコード（１６０１－１，１６０１－２）として記憶される。解析データベース１６００は、例えば、ＲＯＭ４０２、ＲＡＭ４０３、ディスク４０５などの記憶装置によって実現される。
　モデルＩＤのフィールドには、モデルを識別する識別情報が設定される。境界条件の強さのフィールドには、モデルＩＤが示すモデルの境界条件の強さが設定される。例えば、拘束条件の強さが設定される。モデル規模のフィールドには、モデルＩＤが示すモデルのモデル規模が設定される。例えば、要素数が設定される。
　直接法Ａの解析時間のフィールドには、直接法Ａを用いて数値解析した場合の解析時間が設定される。直接法Ｂの解析時間のフィールドには、直接法Ｂを用いて数値解析した場合の解析時間が設定される。反復法Ｃの解析時間のフィールドには、反復法Ｃを用いて数値解析した場合の解析時間が設定される。反復法Ｄの解析時間のフィールドには、反復法Ｄを用いて数値解析した場合の解析時間が設定される。このように、解法には様々な方法があるため、直接法や反復法の種類ごとに解析時間のフィールドが用意される。なお、直接法や反復法に限られるものではなく、他の方法の解析フィールドがあってもよい。
　また、図１６に示すフィールドに限らず、モデルの寸法、体積、メッシュサイズ、要素タイプなども記憶されていてもよい。
　図１７は、システムにおける実績に基づく解法の選択例を示す説明図である。クライアント端末装置１５０１は、例えば、拘束条件の強さ、モデル規模などの評価結果をサーバ１５０３へ送信する。つぎに、サーバ１５０３は、クライアント端末装置１５０１から受信した評価結果に最も近い評価結果を有するレコード１６０１を解析データベース１６００から検索する。具体的に、サーバ１５０３は、評価結果と、解析データベース１６００に含まれるレコード１６０１内の評価結果と、を比較して、類似度をスコア化して最も近い評価結果を有するレコード１６０１を検出してもよい。
　サーバ１５０３は、例えば、検出したレコード１６０１に含まれる解法の解析時間のうち最も短い解析時間の解法を、解析対象モデルの解法として選択する。図１７の例では、反復法Ｃが選択される。サーバ１５０３は、例えば、選択した解法やモデル規模に基づいて、計算資源を決定する。図１の例では、サーバ１５０３－１とサーバ１５０３－２とを用いて解析が行われる。計算資源の決定例については、例えば、図１８を用いて説明する。
　サーバ１５０３－１とサーバ１５０３－２とが、並列処理により、選択した反復法Ｃを用いて解析対象モデルに対して数値解析を行う。
　つぎに、サーバ１５０３は、解析結果をクライアント端末装置１５０１へ送信する。また、サーバ１５０３は、評価結果と、解析時間と、をレコード１６０１として解析データベース１６００へ登録する。
　また、サーバ１５０３は、計算機資源が空いている場合に、解析対象モデルの解法として選択しなかった解法による数値解析を行い、解析時間を解析データベース１６００へ登録する。
　また、検出された最も近い評価結果について、各類似度が低い場合には、サーバ１５０３は、上述したように境界条件の強さとモデル規模に応じた閾値とを比較して解法を選択してもよい。
　図１８は、サーバの状態とメモリ量の管理例を示す説明図である。テーブル１８００は、データセンタ１５０２内にあるサーバ１５０３の状態と資源を管理するための情報を有する。テーブル１８００は、例えば、サーバＩＤ、状態、メモリ量、ＣＰＵ数などのフィールドを有する。
　サーバＩＤのフィールドには、サーバ１５０３を識別する識別情報が設定される。状態のフィールドには、サーバ１５０３の状態を示す情報が設定される。図１８の例では、「使用中」と「未使用」によってサーバ１５０３の状態を表す。また、より詳細に、使用中のＣＰＵ数などによってサーバ１５０３の状態を表してもよい。
　メモリ量のフィールドには、サーバ１５０３が有するＲＡＭ４０３やディスク４０５のメモリ量が設定される。ＣＰＵ数のフィールドには、例えば、サーバ１５０３が有するプロセッサの数が設定される。
　ここで、サーバ１５０３－１～サーバ１５０３－３については、メモリ量とＣＰＵ数が同じであるが、サーバ１５０３－４は、サーバ１５０３－１～サーバ１５０３－３よりもメモリ量及びＣＰＵ数が多い場合を例に挙げる。
　解析対象モデルがｘ～ｚまでの３つの例を挙げる。解析対象のモデル規模は同程度であるとする。また、解析対象モデルｘの拘束条件の強さがａｘであり、解析対象モデルｙの拘束条件の強さがａｙであり、解析対象モデルｚの拘束条件の強さがａｚとする。拘束条件の強さの関係はａｘ＞閾値＞ａｙ＞ａｚである。解析対象モデルｘの解法は反復法が選択され、解析対象モデルｙの解法は直接法が選択され、解析対象モデルｚの解法は直接法が選択される。
　ここで、サーバ１５０３は、直接法を選択した場合に、メモリ量が多い計算資源を選択する。サーバ１５０３は、例えば、サーバ１５０３－４を計算資源として選択する。サーバ１５０３は、例えば、反復法を選択した場合、サーバ１５０３－１～サーバ１５０３－３のいずれかを選択する。また、サーバ１５０３は、例えば、複数の反復法が選択可能な場合、反復法の種類に応じて、サーバ１５０３－１～サーバ１５０３－３のうちの複数のサーバ１５０３によって同時処理を行うようにしてもよい。また、サーバ１５０３は、複数のサーバを用いた並列処理が可能な解法の場合には、複数の計算資源を選択してもよい。さらに、サーバ１５０３は、モデル情報から、効率的な計算に必要なメモリ量等の計算資源を予め見積もり、その結果に応じて適切な計算資源を選択してもよい。
　例えば、解析対象モデルｘについては、例えば、計算資源としてサーバ１５０３－４が選択される。解析対象モデルｙについては、例えば、計算資源としてサーバ１５０３－１が選択される。解析対象モデルｚについては、例えば、計算資源としてサーバ１５０３－２が選択される。このように、選択部５０４は、それぞれの解析対象モデルについての計算資源を選択する。
　また、計算資源の選択のために、テーブル１８００を用いる例を挙げたが、特に限定しない。
（システム１５００による構造解析処理手順例を示すフローチャート）
　図１９及び図２０は、システムによる構造解析処理手順例を示すフローチャートである。ここで、システム１５００は、上述したように、クライアント端末装置１５０１と、サーバ１５０３とを有する。そのため、図１９及び図２０に示す処理については、サーバ１５０３が行ってもよいし、クライアント端末装置１５０１が行ってもよい。
　システム１５００は、クライアント端末装置１５０１によって解析対象モデルのモデル情報１１１を取得する（ステップＳ１９０１）。つぎに、システム１５００は、クライアント端末装置１５０１によって解析対象モデルのモデル規模を評価する（ステップＳ１９０２）。システム１５００は、クライアント端末装置１５０１によって境界条件の強さを評価する（ステップＳ１９０３）。
　つぎに、システム１５００は、解析データベース１６００から、サーバ１５０３によって、解析対象モデルの評価結果に最も近い評価結果を含むレコード１６０１を検索する（ステップＳ１９０４）。システム１５００は、サーバ１５０３によって検索したレコード１６０１から最も解析時間の短い解法を選択する（ステップＳ１９０５）。つぎに、システム１５００は、サーバ１５０３によって特定した解法により数値解析を実行する（ステップＳ１９０６）。
　つぎに、システム１５００は、サーバ１５０３によって解析結果を出力する（ステップＳ２００１）。そして、システム１５００は、サーバ１５０３によって評価結果と解析時間を解析データベース１６００に格納する（ステップＳ２００２）。システム１５００は、計算機環境の負荷が高いか否かを判定する（ステップＳ２００３）。計算機環境の負荷が高い場合（ステップＳ２００３：Ｙｅｓ）、システム１５００は、ステップＳ２００３へ戻る。計算機環境の負荷が高いとは、例えば、処理のないサーバ１５０３がない状態である。
　計算機環境の負荷が低い場合（ステップＳ２００３：Ｎｏ）、システム１５００は、サーバ１５０３によってすべての解法で解析実行したか否かを判断する（ステップＳ２００４）。すべての解法で解析実行していないと判断された場合（ステップＳ２００４：Ｎｏ）、システム１５００は、未実行の解法のうち検索したレコード１６０１の中で最も解析時間が短い解法により数値解析を実行する（ステップＳ２００５）。そして、システム１５００は、サーバ１５０３によって解析時間を解析データベース１６００に格納し（ステップＳ２００６）、ステップＳ２００３へ戻る。ステップＳ２００４において、すべての解法で解析実行したと判断された場合（ステップＳ２００４：Ｙｅｓ）、システム１５００は、一連の処理を終了する。
　以上説明したように、構造解析装置１００は、解析対象のモデルを表す情報に基づきモデル規模を評価し、モデル規模に応じてＦＥＭを用いた連立一次方程式による解法を直接法と反復法とから選択する。これにより、モデル規模によってメモリの使用量を推定できるため、メモリを不足させずに直接法による解析を行うことができる。したがって、精度の向上を図りつつ、解析時間の短縮化を図ることができる。
　また、構造解析装置１００は、モデル規模の評価を、モデル情報とメッシュサイズと要素タイプとから算出された要素数及び／又は節点数及び／又は解析自由度数などに基づいて行う。これにより、モデル規模を簡単に評価することができ、評価にかかる時間の短縮化を図ることができる。
　また、構造解析装置１００は、モデルの境界条件に基づき境界条件の強さを評価し、境界条件の強さに応じてＦＥＭを用いた連立一次方程式による解法を直接法と反復法とから選択する。これにより、境界条件による境界の拘束が安定であるかを判定でき、よりモデルに適した解法が選択できる。したがって、解析時間の短縮化を図ることができる。
　また、構造解析装置１００は、境界条件の強さの評価をモデルの境界条件のうちの拘束条件に含まれる拘束面積、拘束自由度に基づいて行う。これにより、境界の拘束がどの程度で安定するかを簡単に判定することができ、評価にかかる時間の短縮化を図ることができる。
　また、構造解析装置１００は、モデル情報に基づきモデル規模を評価し、境界条件に基づき境界条件の強さを評価し、モデル規模及び境界条件の強さに応じてＦＥＭを用いた連立一次方程式による解法を直接法と反復法とから選択する。これにより、解が収束しない可能性を減らすことができ、かつメモリの使用量が足りるように解法が選択可能となり、解析時間の短縮化を図ることができる。
　また、構造解析装置１００は、モデル規模の評価を、モデル情報とメッシュサイズと要素タイプとから算出された要素数及び／又は節点数及び／又は解析自由度数などに基づいて行う。これにより、モデル規模を簡単に評価することができ、評価にかかる時間の短縮化を図ることができる。また、構造解析装置１００は、拘束条件の強さの評価を、モデルの拘束条件に含まれる拘束面積、拘束自由度に基づいて行う。これにより、境界の拘束がどの程度で安定するかを簡単に判定することができ、評価にかかる時間の短縮化を図ることができる。
　また、構造解析装置１００は、モデル情報に基づいて、解析対象のモデルのモデル規模を評価し、モデル規模に基づいて、モデルの有限要素法を用いた構造解析の数値解析アルゴリズムを選択するとともに、数値解析で用いる計算資源を決定する。これにより、当該モデルの数値解析の特性に応じた、より適した解法と計算資源の選択が可能となり、解析時間の短縮化を図ることができる。
　また、構造解析装置１００は、モデルの境界条件に基づき境界条件の強さを評価し、境界条件の強さに応じてモデルの有限要素法を用いた構造解析の数値解析アルゴリズムを選択するとともに、数値解析で用いる計算資源を決定する。これにより、当該モデルの数値解析の特性に応じた、より適した解法と計算資源の選択が可能となり、解析時間の短縮化を図ることができる。
　なお、本実施の形態で説明した構造解析方法は、予め用意された構造解析プログラムをパーソナル・コンピュータやワークステーション等のコンピュータで実行することにより実現することができる。本構造解析プログラムは、磁気ディスク、光ディスク、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ　Ｓｅｒｉａｌ　Ｂｕｓ）フラッシュメモリなどのコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行される。また、構造解析プログラムは、インターネット等のネットワークを介して配布してもよい。
　１００　構造解析装置
　１０１　解析対象モデル
　１０２　要素
　１０３　節点
　１１１　モデル情報
　５０１　取得部
　５０２　第１評価部
　５０３　第２評価部
　５０５　選択部
　ｍ１　モデル

Claims (11)

  1. 　モデル情報を取得し、
    　取得した前記モデル情報に基づいて、解析対象のモデルのモデル規模を評価し、
    　前記モデル規模の評価結果に基づいて、前記モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する、
    　処理をコンピュータが実行することを特徴とする構造解析方法。
  2. 　前記モデル規模の評価は、前記モデル情報とメッシュサイズと要素タイプとから算出された要素数及び／又は節点数に基づいて行われる、ことを特徴とする請求項１に記載の構造解析方法。
  3. 　解析対象のモデルの境界条件を取得し、
    　取得した前記境界条件に基づいて、前記境界条件の強さを評価し、
    　前記境界条件の強さの評価結果に基づいて、前記モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する、
    　処理をコンピュータが実行することを特徴とする構造解析方法。
  4. 　前記境界条件の強さの評価は、前記モデルの境界条件のうちの拘束条件に含まれる拘束面積、拘束自由度に基づいて行われる、ことを特徴とする請求項３に記載の構造解析方法。
  5. 　モデル情報と解析対象のモデルの境界条件とを取得し、
    　取得した前記モデル情報に基づいて、前記モデルのモデル規模を評価し、
    　取得した前記境界条件に基づき、前記境界条件の強さを評価し、
    　前記モデル規模の評価結果及び前記境界条件の強さの評価結果に基づいて、前記モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する、
    　処理をコンピュータが実行することを特徴とする構造解析方法。
  6. 　前記モデル規模の評価は、前記モデル情報とメッシュサイズと要素タイプとから算出された要素数及び／又は節点数に基づいて行われ、
    　前記境界条件の強さの評価は、前記境界条件のうちの拘束条件に含まれる拘束面積、拘束自由度に基づいて行われる、
    　ことを特徴とする請求項５に記載の構造解析方法。
  7. 　モデル情報を取得し、
    　取得した前記モデル情報に基づいて、解析対象のモデルのモデル規模を評価し、
    　前記モデル規模の評価結果に基づいて、前記モデルの有限要素法を用いた構造解析の数値解析アルゴリズムを選択するとともに、数値解析で用いる計算資源を決定する、
    　処理をコンピュータが実行することを特徴とする構造解析方法。
  8. 　解析対象のモデルの境界条件を取得し、
    　取得した前記境界条件に基づき、前記境界条件の強さを評価し、
    　前記境界条件の強さの評価結果に基づいて、前記モデルの有限要素法を用いた構造解析の数値解析アルゴリズムを選択するとともに、数値解析で用いる計算資源を決定する、
    　処理をコンピュータが実行することを特徴とする構造解析方法。
  9. 　モデル情報を取得し、
    　取得した前記モデル情報に基づいて、解析対象のモデルのモデル規模を評価し、
    　前記モデル規模の評価結果に基づいて、前記モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する、
    　処理をコンピュータに実行させることを特徴とする構造解析プログラム。
  10. 　解析対象のモデルの境界条件を取得し、
    　取得した前記境界条件に基づいて、前記境界条件の強さを評価し、
    　前記境界条件の強さの評価結果に基づいて、前記モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する、
    　処理をコンピュータに実行させることを特徴とする構造解析プログラム。
  11. 　モデル情報と解析対象のモデルの境界条件とを取得し、
    　取得した前記モデル情報に基づいて、前記モデルのモデル規模を評価し、
    　取得した前記境界条件に基づき、前記境界条件の強さを評価し、
    　前記モデル規模の評価結果及び前記境界条件の強さの評価結果に基づいて、前記モデルの有限要素法を用いた構造解析ソルバーの連立一次方程式求解のアルゴリズムを直接法とするか反復法とするか選択する、
    　処理をコンピュータに実行させることを特徴とする構造解析プログラム。
PCT/JP2015/081098 2015-11-04 2015-11-04 構造解析方法、及び構造解析プログラム WO2017077609A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2017548568A JP6645508B2 (ja) 2015-11-04 2015-11-04 構造解析方法、及び構造解析プログラム
PCT/JP2015/081098 WO2017077609A1 (ja) 2015-11-04 2015-11-04 構造解析方法、及び構造解析プログラム
TW105135268A TW201721462A (zh) 2015-11-04 2016-10-31 構造解析方法及構造解析程式
US15/960,744 US11003816B2 (en) 2015-11-04 2018-04-24 Structure analysis device and structure analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2015/081098 WO2017077609A1 (ja) 2015-11-04 2015-11-04 構造解析方法、及び構造解析プログラム

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US15/960,744 Continuation US11003816B2 (en) 2015-11-04 2018-04-24 Structure analysis device and structure analysis method

Publications (1)

Publication Number Publication Date
WO2017077609A1 true WO2017077609A1 (ja) 2017-05-11

Family

ID=58661642

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/081098 WO2017077609A1 (ja) 2015-11-04 2015-11-04 構造解析方法、及び構造解析プログラム

Country Status (4)

Country Link
US (1) US11003816B2 (ja)
JP (1) JP6645508B2 (ja)
TW (1) TW201721462A (ja)
WO (1) WO2017077609A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11061738B2 (en) * 2019-02-28 2021-07-13 Movidius Limited Methods and apparatus to store and access multi dimensional data
CA3157588A1 (en) * 2019-11-06 2021-05-14 Kyle Kosic Methods and systems for the estimation of the computational cost of simulation
TWI764312B (zh) * 2020-10-08 2022-05-11 國立中央大學 基於等效節點割線質量逼近之結構體分析方法、裝置與電腦程式產品

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0573527A (ja) * 1991-09-12 1993-03-26 Nippon Steel Corp 電磁場解析装置

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01145723A (ja) * 1987-08-28 1989-06-07 Hitachi Ltd プログラム生成方法
JP3190697B2 (ja) 1991-02-28 2001-07-23 日本電気株式会社 有限要素法による構造解析システム
TW388921B (en) * 1997-11-28 2000-05-01 Nippon Electric Co Semiconductor process device simulation method and storage medium storing simulation program
US6313490B1 (en) * 1999-03-16 2001-11-06 Micron Technology, Inc. Base current reversal SRAM memory cell and method
JP3648089B2 (ja) * 1999-03-18 2005-05-18 富士通株式会社 設計システムおよび記録媒体
JP4629255B2 (ja) * 2001-03-27 2011-02-09 株式会社富士通長野システムエンジニアリング 構造解析処理方法および構造解析処理プログラム
US7657412B2 (en) * 2001-03-29 2010-02-02 Honda Giken Kogyo Kabushiki Kaisha Method of and apparatus for displaying structure optimizing result
US7289635B2 (en) * 2001-07-20 2007-10-30 Edag Ag Structural noise source predictor
US6813749B2 (en) * 2002-02-12 2004-11-02 The Boeing Company Method, system and computer program product for multidisciplinary design analysis of structural components
JP4583698B2 (ja) * 2002-08-01 2010-11-17 トヨタ自動車株式会社 解析対象物の部品間境界条件の同定装置
CN100429654C (zh) * 2002-11-06 2008-10-29 甲骨文国际公司 用于管理来自单个接口的多层次数据的技术
US7319943B2 (en) * 2002-12-16 2008-01-15 Canon Kabushiki Kaisha Optimum design method, and apparatus, and program for the same
JP4381743B2 (ja) * 2003-07-16 2009-12-09 独立行政法人理化学研究所 境界表現データからボリュームデータを生成する方法及びそのプログラム
JPWO2007086193A1 (ja) * 2006-01-27 2009-06-18 国立大学法人名古屋大学 有限要素法による構造解析方法及びプログラム
JP2008165488A (ja) * 2006-12-28 2008-07-17 Fujitsu Ltd 組立操作性評価装置、方法及びプログラム
JP4973957B2 (ja) 2009-11-02 2012-07-11 独立行政法人科学技術振興機構 有限要素法解析方法、有限要素法解析装置及び有限要素法解析プログラム
US8306793B2 (en) * 2010-06-04 2012-11-06 Livermore Software Technology Corporation Systems and methods of performing vibro-acoustic analysis of a structure
GB201104413D0 (en) * 2011-03-16 2011-04-27 Airbus Operations Ltd Stresses induced by random loading
US20140149090A1 (en) * 2012-11-27 2014-05-29 Autodesk, Inc. Distributed Analysis of Structures
US10008036B2 (en) * 2012-12-10 2018-06-26 Ansys, Inc. System and method for generating a mesh
US11205807B2 (en) * 2014-01-22 2021-12-21 Samsung Electronics Co., Ltd. Computationally screening the stability of battery electrode materials with mixtures of redox couple elements
CN107077518B (zh) * 2014-09-18 2020-10-20 日本制铁株式会社 结构体的部件的变形模式分析方法、采用了变形模式分析方法的结构体的加强方法及记录介质

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0573527A (ja) * 1991-09-12 1993-03-26 Nippon Steel Corp 電磁場解析装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"3-jigen Netsu Kaiseki XCOOL Gijutsu Setsumeisho (Settenho Kaiseki Module Hen", 1 October 1997 (1997-10-01), pages 13, 94 *
YOSHIAKI HARADA: "Yugen Yosoho Yomoyama Banashi Dai 61 Wa Hanpukuho no Fantasy", FEMINGWAY, 30 October 2014 (2014-10-30), Retrieved from the Internet <URL:http://web.archive.org/web/20141030130722/http://femingway.com/?p=2394> [retrieved on 20160115] *

Also Published As

Publication number Publication date
TW201721462A (zh) 2017-06-16
JP6645508B2 (ja) 2020-02-14
US11003816B2 (en) 2021-05-11
JPWO2017077609A1 (ja) 2018-07-26
US20180239854A1 (en) 2018-08-23

Similar Documents

Publication Publication Date Title
Shojaei et al. An adaptive multi-grid peridynamic method for dynamic fracture analysis
JP6143893B2 (ja) ３ｄ点の集合にプリミティブ形状をフィッティングする方法
Georgiou et al. Learning fluid flows
Trobec et al. Parallel scientific computing: theory, algorithms, and applications of mesh based and meshless methods
JP6645508B2 (ja) 構造解析方法、及び構造解析プログラム
US20110144928A1 (en) Vortex extraction apparatus, method, program storage medium, and display system
US11295050B2 (en) Structural analysis method and structural analysis apparatus
US10068037B2 (en) Simulation method and simulation apparatus for continuum motion analysis using a particle method
Benaissa et al. Crack identification using model reduction based on proper orthogonal decomposition coupled with radial basis functions
JP2020184337A (ja) 変形モードについての記述子を使用した物体の構造最適化のための方法
Li et al. Predicting inference latency of neural architectures on mobile devices
Li et al. A MapReduce algorithm to create contiguity weights for spatial analysis of big data
Avril et al. Dynamic adaptation of broad phase collision detection algorithms
Li et al. Flotation froth image texture extraction method based on deterministic tourist walks
US20220067246A1 (en) Simulation method and system
Lu et al. Parallel contact-aware simulations of deformable particles in 3D Stokes flow
JP4830094B2 (ja) 計算機
Kracker et al. Automatic analysis of crash simulations with dimensionality reduction algorithms such as PCA and t-SNE
Li et al. Inference latency prediction at the edge
JP7056452B2 (ja) シミュレーション装置、シミュレーション方法およびシミュレーションプログラム
KR20130101332A (ko) 삼각형 메쉬로 표현된 3차원 물체를 ＤｏＳｕｒｆａｃｅ 표현 방법으로 변환하는 시스템 및 그 방법
JP2006313400A5 (ja)
Zapata et al. A GPU parallel finite volume method for a 3D Poisson equation on arbitrary geometries
Valero-Lara MRF satellite image classification on GPU
US20170091353A1 (en) Computer-readable recording medium, computing machine resource allocation method, and particle simulation apparatus

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15907795

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2017548568

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15907795

Country of ref document: EP

Kind code of ref document: A1