Patent Application
20160196378
The present invention generally relates to computer-aided engineering analysis, more particularly to methods and systems for numerically simulating bi-phase material that changes phase after crossing a directional spatial boundary.
This section is for the purpose of summarizing some aspects of the present invention and to briefly introduce some preferred embodiments. Simplifications or omissions in this section as well as in the abstract and the title herein may be made to avoid obscuring the purpose of the section. Such simplifications or omissions are not intended to limit the scope of the present invention.
Systems, methods and software product for numerically simulating bi-phase material that changes phase after crossing a directional spatial boundary are disclosed. Definition of a directional spatial boundary of a bi-phase material and a finite element analysis (FEA) model containing a plurality of finite elements for representing the bi-phase material are received in a computer system having at least one application module installed thereon. Each of the finite elements is associated with a material identifier that contains first and second sets of material properties corresponding to respective first and second phases of the bi-phase material. The bi-phase material changes from the first phase to the second phase after crossing the directional spatial boundary.
Each of the finite elements is initially assigned with the first set of material properties. A time-marching simulation to obtain numerically-simulated structural behaviors of the bi-phase material moving in the material flow direction using the FEA model is conducted. At each of a plurality of solution cycles during the time-marching simulation, the second set of material properties under the same material identifier is assigned to those of the finite elements determined to have moved across the direction spatial boundary for a directional spatial boundary having instant phase change. Material properties of a finite element located in the transition region are calculated by interpolating first and second set of material properties under the same material identifier for a directional spatial boundary having gradual phase transition.
Numerically-simulated structural behaviors are calculated with the finite elements that are grouped together in accordance with the same material identifier.
In one embodiment, the directional spatial boundary comprises a plane that partitions the three-dimensional space to two regions. Bi-phase material changes phase instantly after crossing the plane. The plane is defined by first and second nodes with the first node located on the plane and the second node being used to define the material flow direction.
In another embodiment, the directional spatial boundary comprises first and second parallel planes for beginning and ending of the phase transition of the bi-phase material, respectively. Two parallel planes are defined by first and second nodes with the first node located on the first plane and the second node located on the second plane. The material flow direction is defined as a direction from the first node to the second node.
One of the advantages of the present invention is related to computational efficiency. In explicit solution technique used in a time-marching simulation, the required smallest time step size for the first phase and the second phase of the bi-phase material can be drastically different. Material properties of the two phases are defined under the same material identifier, which allows the finite elements to be grouped together to employ certain numerical scheme to speed up the computation, for example, subcycling.
Other objectives, features, and advantages of the present invention will become apparent upon examining the following detailed description of an embodiment thereof, taken in conjunction with the attached drawings.
These and other features, aspects, and advantages of the present invention will be better understood with regard to the following description, appended claims, and accompanying drawings as follows:
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will become obvious to those skilled in the art that the present invention may be practiced without these specific details. The descriptions and representations herein are the common means used by those experienced or skilled in the art to most effectively convey the substance of their work to others skilled in the art. In other instances, well-known methods, procedures, and components have not been described in detail to avoid unnecessarily obscuring aspects of the present invention.
Reference herein to “one embodiment” or “an embodiment” means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the invention. The appearances of the phrase “in one embodiment” in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. Further, the order of blocks in process flowcharts or diagrams representing one or more embodiments of the invention do not inherently indicate any particular order nor imply any limitations in the invention.
Embodiments of the present invention are discussed herein with reference to
Process 100 starts, at action 102, by receiving a definition of a directional spatial boundary of a bi-phase material represented by a finite element analysis (FEA) model in a computer system (e.g., computer system 600 of
The FEA model containing a plurality of finite elements (e.g., solid elements, plate elements, beam elements). Example finite elements include, but are not limited to, quadrilateral element, triangular element, hexahedral element and tetrahedral element. Each of the finite elements is associated with a material identifier that contains first and second set of material properties corresponding to respective first and second phases of the bi-phase material.
Next, at action 103, the material flow direction and type of the directional spatial boundary are determined from the received definition. There are two different types of directional spatial boundary. The first type (an example shown in
In one embodiment,
The first example directional spatial boundary comprises a plane 400 as shown in
In another embodiment, a second example directional spatial boundary in the material flow direction 310 is shown in
One example representation of the second example directional spatial boundary is shown in
At action 104, process 100 assigns the first set of material properties to all of the finite elements in the FEA model. In other words, every finite element in the FEA model that represents the bi-phase material is initialized as first phase.
Next, at decision 105, it is determined what type of directional spatial boundary is used for a time-marching simulation. When instant phase change type is determined, at action 106, a time-marching simulation is conducted to obtain numerically-simulated structural behaviors of the bi-phase material moving in the material flow direction using the FEA model.
The time-marching simulation contains a number of solution cycles (i.e., time steps). At each solution cycle, finite elements already assigned with the first set of material properties are checked to determine whether any of them has moved across the location of the directional spatial boundary (e.g., first node of the example shown in
To any finite element so determined, the second set of material properties under the same material identifier is assigned. The numerical structural behaviors of the bi-phase material are calculated. The calculations are performed with finite elements that are grouped together according to respective material identifiers. The finite elements having the same material identifier are gathered in one group. Since bi-phase material is defined under the same material identifier with two different material properties, the calculations can be performed with high efficiency, for example, keeping vectorization.
When gradual phase transition is determined at decision 105, process 100 moves to action 108. A time-marching simulation is conducted to obtain simulated structural behaviors of the bi-phase material moving in the material flow direction. At each solution cycle during the time-marching simulation, each finite element is checked whether it is in the transition region (e.g., the transition region between the first node 320a and the second node 320b shown in
Simulated structural behaviors can then be calculated using all of the finite elements in the FEA model. Since the finite elements under same material identifier are grouped together, computation efficiency can be achieved both for performing the interpolation and for calculating structural behaviors.
According to one aspect, the present invention is directed towards one or more computer systems capable of carrying out the functionality described herein. An example of a computer system 600 is shown in
Computer system 600 also includes a main memory 608, preferably random access memory (RAM), and may also include a secondary memory 610. The secondary memory 610 may include, for example, one or more hard disk drives 612 and/or one or more removable storage drives 614, representing a floppy disk drive, a magnetic tape drive, an optical disk drive, etc. The removable storage drive 614 reads from and/or writes to a removable storage unit 618 in a well-known manner. Removable storage unit 618, represents a floppy disk, magnetic tape, optical disk, etc. which is read by and written to by removable storage drive 614. As will be appreciated, the removable storage unit 618 includes a computer readable storage medium having stored therein computer software and/or data.
In alternative embodiments, secondary memory 610 may include other similar means for allowing computer programs or other instructions to be loaded into computer system 600. Such means may include, for example, a removable storage unit 622 and an interface 620. Examples of such may include a program cartridge and cartridge interface (such as that found in video game devices), a removable memory chip (such as an Erasable Programmable Read-Only Memory (EPROM), Universal Serial Bus (USB) flash memory, or PROM) and associated socket, and other removable storage units 622 and interfaces 620 which allow software and data to be transferred from the removable storage unit 622 to computer system 600. In general, Computer system 600 is controlled and coordinated by operating system (OS) software, which performs tasks such as process scheduling, memory management, networking and I/O services.
There may also be a communications interface 624 connecting to the bus 602. Communications interface 624 allows software and data to be transferred between computer system 600 and external devices. Examples of communications interface 624 may include a modem, a network interface (such as an Ethernet card), a communications port, a Personal Computer Memory Card International Association (PCMCIA) slot and card, etc. The computer 600 communicates with other computing devices over a data network based on a special set of rules (i.e., a protocol). One of the common protocols is TCP/IP (Transmission Control Protocol/Internet Protocol) commonly used in the Internet. In general, the communication interface 624 manages the assembling of a data file into smaller packets that are transmitted over the data network or reassembles received packets into the original data file. In addition, the communication interface 624 handles the address part of each packet so that it gets to the right destination or intercepts packets destined for the computer 600. In this document, the terms “computer program medium” and “computer usable medium” are used to generally refer to media such as removable storage drive 614, and/or a hard disk installed in hard disk drive 612. These computer program products are means for providing software to computer system 600. The invention is directed to such computer program products.
The computer system 600 may also include an input/output (I/O) interface 630, which provides the computer system 600 to access monitor, keyboard, mouse, printer, scanner, plotter, and alike.
Computer programs (also called computer control logic) are stored as application modules 606 in main memory 608 and/or secondary memory 610. Computer programs may also be received via communications interface 624. Such computer programs, when executed, enable the computer system 600 to perform the features of the present invention as discussed herein. In particular, the computer programs, when executed, enable the processor 604 to perform features of the present invention. Accordingly, such computer programs represent controllers of the computer system 600.
In an embodiment where the invention is implemented using software, the software may be stored in a computer program product and loaded into computer system 600 using removable storage drive 614, hard drive 612, or communications interface 624. The application module 606, when executed by the processor 604, causes the processor 604 to perform the functions of the invention as described herein.
The main memory 608 may be loaded with one or more application modules 606 that can be executed by one or more processors 604 with or without a user input through the I/O interface 630 to achieve desired tasks. In operation, when at least one processor 604 executes one of the application modules 606, the results are computed and stored in the secondary memory 610 (i.e., hard disk drive 612). The status of the finite element analysis is reported to the user via the I/O interface 630 either in a text or in a graphical representation.
Although the present invention has been described with reference to specific embodiments thereof, these embodiments are merely illustrative, and not restrictive of, the present invention. Various modifications or changes to the specifically disclosed exemplary embodiments will be suggested to persons skilled in the art. For example, whereas the directional spatial boundary has been shown and described as a plane, other forms of physical objects can be used instead, for example, a box. In summary, the scope of the invention should not be restricted to the specific exemplary embodiments disclosed herein, and all modifications that are readily suggested to those of ordinary skill in the art should be included within the spirit and purview of this application and scope of the appended claims.
