Cycle monitoring method and apparatus.

BACKGROUND OF THE INVENTION

Field of the Invention

0001This invention is directed to a method and apparatus for monitoring and analyzing, on an on-line basis, the fatigue caused cumulatively in components and systems subjected to fluid flow, theremal and pressure transients, particularly nuclear power plant components and systems. More particularly, the invention is directed to presenting to plant engineers a represention of the margin between current fatigue accumulation and a value of and a value of fatigue accumulation that indicates need for repair, replacement or refubrishment of the component or system is necessary. The invention also encompasses techniques for providing accurate determination of stresses at plant locations ciritical to fatigue, the determination being made through measurements of plant process and response parameters.

Description of the Related Art

0002Components in a nuclear power plant are subjected throughout their 40 year nominal life to stresses induced by changes in process fluid temperature and pressure. It is necessary for the operators and engineers of a nuclear power plant to be constantly aware of the status of criti­cal plant components and systems with respect to fatigue. Presently, this is done by the following means: measure­ments of plant process and response parameters (tempera­ture, pressure, flow level, valve status, setpoint status, for example) are measured and the data is recorded by a plant computer. This data is used by plant engineers to make subjective judgements regarding certain types of transients. The types of transients that must be recorded are specified in the plant technical specifications, which are governed by regulations, e.g., by the U.S. Nuclear regulatory Commission (NRC) if in the United States. The technical specifications allow a specified number of each type of transient during the 40-year life. The prior art methods do not perform fatigue calculations to determine changes in life expectancy, but make educated guesses concerning when a component should be removed from service, based on the conservative transient counts allowed by regulators and subjective determinations, by engineers, as to whether a particular transient should be counted. The transients in a nuclear power plant fall into two general categories; (i) normal transients which occur during plant start-up fand plant power output demand changes and (ii) abnormal transients which occur when a plant component or system is functioning abnormally. If the allowable number for one or other type of transients is exceeded, justifi­cation must be provided to the Regulating Agency as to why the plant is safe for continued operation. This may include revision of the stress analysis which formed the original design basis for the plant.

SUMMARY OF THE INVENTION

0003It is an object of the present invention to facilitate decisions regarding electric power plant life extension.

0004The present invention provides a systematic method and apparatus for accurately measuring fatigue in critical plant components due to all significant thermal and pressure transients; advantageously, the present invention compiles a database that can be used in support of utility requests to regulatory authorities for permis­sion to extend the operating life of electric power plants or used to eliminate the requirements for counting plant cycles; preferably the present invention provides a system­atic method and apparatus for monitoring critical plant parameters and for compiling a database for maintenance, repair and replacement decisions for critical plant compo­nents; optionally the present invention provides a system­atic method and apparatus for monitoring critical plant parameters and for compiling a database that can be used to justify non-design basis operating modes in order to provide maximum plant flexibility and optimization of plant operation.

0005To achieve the above objects, the present inven­tion provides a system and method that acquires, prepares and analyzes analog and/or digital signals from component sensors in a process control plant such as a nuclear power plant. The system continuously monitors the sensors and records steady-state and transient phenomena. Compression of the steady state data is performed to reduce the storage load. Transient data is recorded at a higher sampling rate, so that the parameters of the transients, such as the maximum value, can be determined for stress analysis. The stress analysis determines pressure and thermal stresses which are used to calculate a usage factor which indicates the aging due to stress and thus to indicate the remaining life of a component.

0006The invention in its broad form resides in an on-line monitoring method of assessing cumulative fatigue on a component in use by monitoring cycles using a sensor which produces steady state and transient signals, said method comprising the steps of: (a) sampling a sensor signal from the sensor; the steps of time compressing sampled steady state signals; (b) producing a compressed transient from the sensor signal when the sensor signal includes a transient; and (c) computing fatigue on the component from the compressed steady state and transient, thereby to indicate ahead of time from an expected fatigue tolerance of the component, as to when the component might need refurbishing.

BRIEF DESCRIPTION OF THE DRAWING

0007A more detailed understanding of the invention may be had from the following description of a preferred embodiment, given by way of example and to be understood in conjunction with the accompanying drawing wherein: <ul id="ul0001" list-style="none"><li>Fig. 1, including Figs. 1A and 1B, is a function­al hardware and software block diagram of a preferred embodiment of the present invention;</li><li>Fig. 2, including Figs. 2A and 2B, is a flowchart showing the relationship between the software modules of the system.</li><li>Fig. 3, including Figs. 3A and 3B, is a flowchart of the data acquisition module 30;</li><li>Fig. 4, including Figs. 4A-4C, is a flowchart for the data integrity check 86;</li><li>Fig. 5 is a sketch of a typical component showing monitored and unmonitored/undermonitored locations for a typical transient;</li><li>Fig. 6 is an example of a temperature transient;</li><li>Fig. 7, including Figs. 7A-7C, is a flowchart of the process for determination of transient recording duration;</li><li>Fig. 8 is a flowchart of the transient parameter generation module;</li><li>Fig. 9, including Figs. 9A-9C, is a flowchart of a typical analysis module;</li><li>Fig. 10 is another sample transient;</li><li>Fig. 11 is a flowchart of the method for deter­mining stress age or remaining life;</li><li>Fig. 12 is a sample fatigue curve;</li><li>Fig. 13 illustrates the process of obtaining transient thermal stress;</li><li>Fig. 14 is an example component for transient thermal stress calculation; and</li><li>Fig. 15 is a transient example.</li></ul>

DESCRIPTION OF THE PREFERRED EMBODIMENT

0008The present invention preferably uses a dedicated minicomputer for all data acquisition, computation and analysis activities. However, the system could be imple­mented on an already existing computer at the plant, if that computer had adequate data acquisition and computation capability. The description is directed to calculating fatigue at specific locations in a nuclear power plant, however, this method and system can be used for calculation of fatigue crack growth. The changes necessary to incorpo­rate fatigue crack growth can be performed by one of ordinary skill in the art, based on Hertzberg, R. W., <u style="single">Deformation</u><u style="single">and</u><u style="single">Fracture</u><u style="single">Mechanics</u><u style="single">of</u><u style="single">Engineering</u><u style="single">Materi­als</u>, John Wiley and Sons, 1976. The apparatus and methods described herein are not limited to nuclear power plants, but can be applied to any process plant, or even individual components or systems subject to fatigue loadings due to fluid flow, thermal and pressure transients.

0009Fig. 1 is a functional block diagram of the preferred embodiment of the integrated hardware and soft­ware of the invention. The hardware portion of the system consists of a signal input section 10, a sampling and analysis computer 12 and a set of output peripherals 14. It is possible to substitute a slower speed computer for the sampling and analysis computer 12 if only sampling is performed and to provide a separate remotely located analysis computer. The input section 10 provides pressure, temperature and fluid flow data produced by plant analog sensors 16 and/or digital sensors 18, typically found in nuclear power plants. These sensors 16 and 18 read the process signals and response parameters for the components and locations being monitored. Analog signals read by the analog sensors 16 are processed through signal conditioners 20, which include a sensor signal multiplexer for selecting a single sensor signal at a time. The signal conditioners 20 are also typically found in nuclear power plants, and convert the sensor signals to the 0 to 10 volt signals required by digitizers 22 and 24. The analog signals proceed either to the low speed digitizer 22 or the high speed digitizer 24, depending on the frequency of the signals. Suitable digitizers 22 and 24 are available from Tustin as model number X1500. The decision to use the low spped digitizer 16 or the high speed digitizer 14 for a given input signal is made by the design engineer when the system is installed depending on the highest frequency expected from the sensor. The high speed digitizer 24 has a variable capacity with a typical sampling rate of about 200,000 samples per second. Digital signals are produced by the digital sensors 18 and processed through a buffer and counter unit 26 such as a digital input/output board, available as a discrete isolated input card model 1585 from Tustin. The function of the buffer and counter unit is to determine the location of the signal and whether it is on or off (one or zero). Once the signals are processed through the buffer and counter unit 26 in the case of digital signals, or through the low speed digitizer 22 or high speed digitizer 24 in the case of analog signals, the signals are properly digitized and ready for input into the computer 12. The host computer 12 is preferably a Data General computer model MV-4000 having at least two mega­bytes of main memory and a running real time virtual operating system 28 such as Data General's OSVS operating system. This sytem will handle up to 100 channels of sensor data, if more channels are needed, faster and bigger computer will be necessary. The real time operating system 28 acts as the manager and causes loading of appro­priate software modules 30-38. The real time operating system 28 also controls the link between the computer 12 and the output peripherals 14. The output peripherals 14 include typical units associated with minicomputers and are available from Data General such as a 300 megabyte disk drive 40 for short term data storage, cathode ray tube 42 or a printer 44 for producing results of fatigue analysis, a plotter 46 for showing changes in component life expec­tancy and a magnetic tape unit 48 for archival data storage.

0010The software is modular and expandable such that other component tracker modules 38 can be added as required. Other component tracker modules can be based on fatigue crack growth, or any other technically valid theoretical construct for determination of accumulated fatigue, including automated technical specification transient cycle counting. The preferred embodiment illlustrates application of the systeme using fatigue usage factor methodology. A flowchart showing the relationship between the software modules 30-36, is shown in Fig. 2. The sampling module 30 can be considered the manager of the input section 10 in Fig. 1. The sampling module 30, discussed in detail with respect to Fifs. 3-7, combines an input data base 60 (residing in computer 12) with the input signals into a data structure such as array, then performs a data integrity check. the input database 60 contains: engineering units for converting the sensor samples into pressure, temperature, etc.; threshold limits that determine when a transient has occurred, requiring fatigue analysis; component identification information, including sensor locations, ranges, etc.; loads, such as pipe dead weight and seismic loads; fatigue curves depicting the effect of stress levels on fatigue; stress kernels which describe the stress response of a component to a unit charge in a process parameter such as temperature or pressure; and materials properties. The sampling module 30 stores the data in a predetermined size transient and steady state time history database 62. The transient and steady state time history database 62, when full, is used as input to the tranisent analysis module 32, discussed in detail with respect to Figs. 9 and 10. The module 32 computes transient parameters that result in compilation of a transient parameter database 64. The transient parameter database 64 provides the input required for the fatigue calculation software module 34, discussed in detail with respect to Figs. 11-15. Module 34 which computes the component fatigue usage factor for each component and location being monitored and provides the relevant fatigue data to an output database 66. The display software module 36 provides provisions for operator display. The display module 36 outputs component names, locations and remaining life information using tabular representations, graphs and other data representation schemes useful to plant operators. The display module 36 will not be discussed in detail herein, since the production of a suitable display depends on operator preference, plant design, etc., and can be provided by one of ordinary skill in the art.

0011Modules 50 and 52 are other possible modules that cna count transients by merely keeping a count of the transients that a component experiences or simply logging the transients and operating information for operator and engineering reference.

0012The details of operation of the sampling software module 30 will be explained with respect to Figs. 3, 4 and 5. Referring to Fig. 3, the input database 60 contains an array indicating sensor type. This array associates the many input signals channels with particular sensors in the plant. Sensor signals for pressure, temperature, flowrate and level are recorded along with control signals such as valve position and pump status (on or off) and response signals such as strain or displacement. The recording operation is preferably performed one sensor at a time, however, it is possible to sample all sensors 16 and 18 and then perform integrity checks on the samples. First, if the signal sampled is an analog sensor, a sensor channel is selected 82 by incrementing the input selected by the multiplexer in the signal conditioning unit 20. This performed by incrementing a channel pointer and outputting the pointer to the multiplexer. The appropriate digitizer is then activated 84 to sample the signal and the computer 12 stores the sample in temporary memory. If the sensor is a digital signal, the computer merely reads from the appropriate location. In the data integrity check 86, temperature, pressure, flowrate, and level are compared with the engineering units and allowable sensor ranges of the input database 60 to assure the sensors are working correctly. The allowable ranges are determined by the laws of nature and the capacities of the plant. For example, referring to Fig. 4 which provides the details for representative integrity checks, using the established system of units, temperatures in a power plant would never be expected to be below 0°F. If a sensor indicated a measured temperature 110 below 0°F (118), the data integrity check would cause the operator to be advised 88 of the anomaly. Similarly, if a sensor produced a temperature signal 112, indicating 120 a temperature above a certain value (114) which would be specific to a certain plant, the data integrity check 86 would again cause the operator to be advised 88 of the anomaly. The same procedure is applicable to pressure 128 and 130, flowrate 142 and level 144 sensor signals. Each sensor of each type can, of course, have its own allowable range. Data integrity checks can also be designed for control 146 and response 148 signals. These are not as widely used and are based on binary decisions rather than ranges. As an example, if a signal 128 for a pressure indicates a value within 134 an allowable range, but a control signal indicates 140 that the valve providing that pressure is closed, an inconsistency exists and the data integrity check 86 would alert 88 the operator to the problem. The data integrity check would alter 88 the operator to the problem including any short circuits or if data is lost 160.

0013After the integrity check is performed, a determination is made 89 whether data from the sensor for a predetermined time window has been collected. If so, the samples are compared 90 (Fig. 3) to a threshold level that determines whether a transient has occurred. If the sensor values are within the threshold range, rotating storage 92 is updated 94, from which steady state statistics are computed 96. If the sensor values are outside the transient threshold range, the duration of the recording period for the transient is determined and the transient is recorded 98. Both the steady state statistics and transients are stored in the transient and steady state time history database as, for example, a list of temperatures and the times at which the temperatures were recorded.

0014To discuss the steps 90-96 in greater detail, a typical piping section 160 and a typical transient 162 are shown in Figs. 5 and 6. The pipe 160 has two critical locations 164 and 166 about which fatigue information is desired. One of these 164 is monitored for temperature, pressure, flowrate, level and control and response signals, as previously described. The other locaiton 166 is not monitored, or is insufficiently monitored, that is, only pressure and level are available at location 166. The typical transient 162, shown in Fig. 6, represents the change in temperature over time for a transient detected at a monitored critical location 164. Signals for pressure, level and flowrate also vary at the monitored critical location 164, but use of the temperature signal is used as an illustration, while substantially identical operations are performed for the other types of signals. Assuming that the input signal at the monitored location pass the data integrity check 86 (Fig. 3), the next step is to check 90 for threshold level violations. Data is sampled and collected for a time window of preferably 20 seconds for each sensor. As long as the values for pressure, temperature, flowrate and level for all monitored locations stay within threshold levels indicated in the input database 60 during that time period, all that is done is to update 94 the rotating storage 90 and compute steady state statistics 96. The rotating storage 90 is preferably a data log of the most recent 1/2 hour of steady state data.

0015Steady state statistics can consist of means, maximum, minimum, and other statistical data for the signals for the duration of the time window. This steady state statistical data becomes part of the transient and steady state time history database 62. This database can be kept for the life of the plant, if desired on magnetic tape.

0016If any monitored temperature, pressure, flowrate, or level signal in the entire plant exceeds the threshold values stored in the input database 60, transient recording starts for all monitored locations in the plant and continues until a steady state is reestablished for a predetermined period. The threshold values of the input database 60 are determined by engineering calculation and judgments made by the plant designers and would be values at which fluctuations will cause negligible fatigue. For example, in a typical plant, temperatures may fluctuate by plus or minus "T" degrees Fahrenheit and cause negligible fatigue while temperature fluctuations larger than that range could induce stresses of magnitudes that could begin to cause fatigue.

0017Fig. 7 illustrates the steps 89-98 in greater detail when a particular channel or channels are designated as trigger channels for triggering recording. Not all channels are designated as triggering channels and the designations would be made by the plant designer based on failure likelihoods or regulatory requirements. First, the mean value for a process parameter (for example, temperature) is calculated 170 over a given time window, typically 20 seconds. If the difference between this calculated mean and the mean established for this signal after the end of the previous transient exceeds the threshold level 172, transient recording begins 174 for all channels. If the threshold is not exceeded for the process parameter (temperature) being checked, additional signals are read 176 (for example, pressure) until all process parameter signals have been read for that time window. If none of the parameters (temperature, pressure, etc.) exceeds the threshold, rotating storage is updated 90, as described earlier, and a new time window is read 170.

0018Once a temperature transient 162, or any other temperature, pressure, flowrate or level signal has exceeded 172 the threshold, the operation switches from the steady state recording already described to transient recording. Transient recording retains a great deal more data than steady state recording. In transient recording, all sampled data is retained during each 20 second time window until compressed in the transient analysis module 32 which will be explained later. Thus, rather than 4 or 5 datapoints being kept as in steady state recording, hundreds of data points are kept for each sensor signal. If a sensor signal is sampled at 10 samples per second, 200 samples will be kept for that parameter in a 20 second time window. Transient recording continues until the algorithm to determine recording duration indicates that steady state recording can begin again.

0019After transient recording begins for all channels 174, the base mean is calculated 176 from the first 20 seconds of data designated as the first window. The mean of the signal for current time window for pressure, temperature, flowrate and level signals is calculated 178 and compared 180 to the value for the corresponding signal during the previous time window. When the mean for the current time window is different from the mean for the previous time window by less than the threshold for a given channel, the mean for the previous time window becomes the base mean for comparison for that channel. That is, the base mean for the previous time window is not changed. Comparisons for future time windows of that channel are compared 180 to the base mean for that channel.

0020If the difference between a current time window mean and the base mean equals or exceeds 182 the threshold the previous base mean is invalid and a new one must be established 176. After a threshold check has ocurred, a determination 184 is made concerning whether all parameters for this window have been compared. When base mean values can be established and held 186 for all pressure, temperature, flowrate and level signals simultaneously for a given period, 10 minutes, for example, transient data collection stops and the system is returned to the steady state mode 190.

0021An example transient signal will clarify the operations described above. In Fig. 6, the transient signal 162 stays within the threshold until it violates the threshold at point 192. At point 192 high speed recording begins. That is, rotating storage is not updated and all samples get stored. Between point 192 and period 194, the signal is continually violating the threshold around each new base mean value. During period 194, the signal stays within the threshold for approximately 2 1/2 minutes until point 196. Since this is less than the 10 minute threshold re-establishment period, high speed recording continues. Between point 196 and period 198 once again, the threshold is continually violated. During period 198 the signal stays within the threshold range for the 10 minute period and high speed recording would stop at time 27 minutes and resumption of the updating of rotating storage would occur. After steady state recording is again initialized, if the threshold around the base mean value for one channel is violated, the system again switches to transient mode high speed recording for all channels. While in the transient mode, data is stored in the steady state and time history database 62 at the higher speed sampling rate established for the system.

0022The operation transient parameter analysis module 52 will now be explained with reference to Fig. 7. Using data from the input database 60, and the real time data being provided in the transient and steady state time history database 62, each channel is examined to determine 200 the signal type. The signal type is indicated in the array associated with the sensor identified. The signal types are pressure (P), temperature (T), flowrate (Q), level (L), control (CS) and response (RS). Examples of control signals (CS) are valve position (percent open), or pump status (percent of full load) while examples of response signals (RS) are strain or displacement. The values for the control signal and response signals simply pass 202 and 204 to the transient parameter database 64 for future use by the plant operators or engineers, or are used in other component tracker modules. If the signal is pressure, temperature, flow, or level, a data compression process takes place 206.

0023During data acquisition by the sampling software module 36 (Fig. 2) data is sampled at rates high enough to assure gathering of all relevant data. This process results in accumulation of more data than is needed at some times, but these time periods cannot be determined until after the data is collected. The compression module 206 keeps relevant data and eliminates excessive data. For each type of signal (pressure, temperature, flowrate, level) this is done in a different way. For temperature signals in nuclear power plants, a data point (including temperature and time) is typically kept for every 5°F change in temperature. For pressure, a datapoint is typically kept for about 100 PSI change in pressure. Similarly, based on plant requirements for flowrate signals, a datapoint can be kept for each change in flowrate of some specified magnitude, and for level signals a datapoint can be kept for a certain specified change in level. Other types of processes would require different compresion values. Since pressure and temperature changes are often the most important for fatigue determination, it is also possible to associate each flowrate or level signal with a certain pressure or temperature signal. In this case, each time a datapoint is saved for a given temperature or pressure signal, the corresponding datapoints for any associated flowrate or level signal, are also saved. In addition, during the compression process, statistical data for the recorded transient are calculated.

0024The reduced set of pressure, temperature, flowrate and level signal data points and control and response signal data points and statistical data, form the transient parameter database 64.

0025Fig. 9 illustrates how a temperature signal from one channel is processed through the transient parameter generation module 32 which performs the data compression. Other channel temperature signals and other process signals (pressure, level, etc.), are handled similarly. While the sampling module 30 is a real time program, the transient analysis module 32 is an on-line module. It runs automatically when the transient and steady state time history database 62 contains a predetermined number of samples that is, when it is full. The first transient sample of temperature data is first read and stored 220 as the base transient sample.

0026A transient sample with respect to Fig. 9 will now be defined which is different from the standard definition well known to those skilled in the art. If data is sampled (using the traditional definition of a sample which is a single digitization of an analog value at a time instant) at 10 samples per second, a transient sample in Fig. 9 could be of only one of the 10 points of temperature data in the one second time window. That is, a transient sample would be one of the smaples in a ten sample time window. A transient sample also contains the associated time at which the temperature was read, along with the value of any flowrate or level data which has been linked to the temperature data through the input database 60 sampled at the same time as the temperature data point. The one out of ten samples that becomes the transient sample is chosen by the process of Fig. 9. This definition of a transient sample will be used throughout only Fig. 9.

0027After the first sample is read and stored 220 as the base transient sample, a next or current sample is read and three things can happen. First, the temperature in the current transient sample can equal 224 the temperature of the base transient sample. If this happens, the data of the current sample is discarded and the next or a new current sample 138 is read. Second, the temperature of the current sample can be greater 226 than the value of the temperature for the base transient sample. Third, the temperature of the current sample can be lower 228 than the temperature of the base sample. The process occurring after the second and third possibilities are very similar and the process occurring after the second possibility will now be explained in detail while the process occurring after the third possibility will not be explained in detail, for simplicity purposes.

0028When the value of temperature for the current transient sample is greater than the base transient sample 226, two possibilities exist. First, the current sample can be larger by less than some value "T", for example, 5°F. If this is the case, the data for the subsequent sample is discarded and the next or a new sample is read 222. The second possibility is that the temperature for the current transient sample is greater than the value for the first transient sample by greater than or equal to 5°F. If this is the case, the sample is stored in the transient parameter database as a transient sample and becomes the base transient sample. Then, the next or a new current sample is read 232. If the temperature of this latest current sample is equal to the value of the sample just stored, the current sample is discarded 234 and a new current sample is read 232. If the value of the subsequent sample 146 is less than 5°F greater 236 than the stored sample, the data for this sample is discarded and a new current sample is read 232. If the value of the subsequent sample is greater than or equal to the value fo the stored sample plus 5°F, the subsequent sample is stored 232 as a transient sample and becomes the base transient sample for comparison with future samples. If the value for the current sample is less than 238 the value for the previous stored transient sample or base transient sample, samples after the previous stored sample and before the current sample 146 in question are reviewed 240, by using the transient and steady state time history database 62 and reading starts again at one point past the maximum. That is, the time associated with previous base transient sample is used as a start point and the time of the current transient sample is used as the end point of a time period used for scanning the data in the transient and time history database 62. The maximum temperature in this time period is found 240 and this sample is stored in the transient parameter database 64 as a transient sample.

0029If any attempt to read a subsequent sample anywhere in the above-described process fails 242 or 244 this indicates that the data in the transient and steady state time history database 62 has all been read. At this point, the desired statistics 246 are calculated on the transient parameter database 64 which now contains all data stored from the processes in 136, 166 and 168. The statistical data is then added 248 to the transient parameter database 64.

0030For example, assume that the transient and steady state time history database 62 includes temperature points (transient samples), illustrated in Fig. 10 for a positive going transient. The first transient sample 270 is the base transient sample and is stored (5°F) in the transient parameter database 64 along with its time of occurrence (2 seconds). The subsequent transient samples are then compared with the base transient sample until the 5°F threshold is attained at point 272 with the intervening transient samples 274 being discarded. The temperature (10°F) and the time (20 seconds) are stored in the database 64 and point 272 becomes the new base transient sample. The process continues with point 276 becoming the base transient sample. When the point 278 is reached where the temperature is less than that at point 276, the point 278 is stored and the scanning of the transient and steady state time history database 62 begins at time equal to 32 seconds (point 276) and ends at time 44 seconds (point 278). The scan looks for the maximum transient sample (point 280) by comparing the temperature values and stores the temperature (17°F) and time (36 seconds) in the transient parameter data base 64. When all data is exhausted, the statistics for the transient are determined. That is, the minimum is indicated as point 270 (temperature equal 5°F at 2 seconds), the maximum is indicated as point 280, the initial value as point 270, the final value as point 278, the time t as 42 seconds and the number of cycles as one. Similar storage and analysis occurs if the transient is a negative going transient.

0031Operation of the fatigue calculation software module 34 will be explained with respect to Fig. 11. The input database 60 contains information on each channel of pressure data, each channel of temperature data, each channel of flowrate data, each channel of level data and each control and response signal contained in the transient parameter database 60, pertaining to global to local transformations. First, a global to local transformation is performed 290 between monitored and unmonitored points (See Fig. 5). In this example, there are two critical location possibilities, monitored 164 and unmonitored 166. A monitored critical location has monitoring equipment attached which monitors all relevant process, control, and response parameters. At any monitored critical location 164, any of pressure, temperature, flowrate, level, control signals and response signals may be monitored as necessary. Not all monitored critical locations 164 require that all parameters be monitored. A monitored critical location 164 is considered so if all parameters relevant to fatigue calculations at that location are monitored. At the unmonitored locations 166, either not enough parameters are monitored to fully characterize the fatigue response of the system, or no monitoring is done at all. The input database 60 physically correlates (distance, etc.) of the critical locations 164 with the process and response channels of the transient parameter database 60. For each critical location 164 and 166 (Fig. 5), the input database 60 contains the following information: (1) the parameters monitored at that location, and (2) whether the location is monitored or unmonitored as defined earlier. For unmonitored critical locations 166, a transformation of global transient parameters from a nearby monitored location 164 must be applied to obtain the information necessary to perform fatigue calculations at the unmonitored location 166. In the example, pressure and level are monitored at the unmonitored location 166, but temperature and flowrate are also needed. The transformation 290 of Fig. 9 uses the well known laws of fluid mechanics and heat transfer to extrapolate the values of temperature and flowrate at the unmonitored location 166 in Fig. 5 using the known values of these parameters at the monitored location 164.

0032For example, assume the distance between location 164 and 166 in Fig. 5 is 10 meters, and the pipe diameter is 10 millimeters (values stored in input database 60). These dimensions are smaller than might actually be the case in a plant situation, but are chosen for illustration. We assume that the pressure is monitored at location 164 but not monitored at location 166 and that the flow is directed from 164 to 166. Because of friction losses, the pressure at 166 is less than at 164, and the equation governing the pressure drop is: (P164/ρ + α164 V164²/2 +gz164)-(P166/ρ + α166 V166/2 + gz164) = h₁ + h_1m      (1)<maths id="math0001" num=""><img file="EP0280489A2_D0001.tif" /></maths> For incompressible flow, fluid density (ρ) is constant, and since the pipe area is constant, the flow velocity (V) is constant. The pipe horizontal so the vertical height (z) is constant. Assuming fully developed flow in a straight section of pipe, h1_m = 0. The above equation thus reduces to:<maths id="math0002" num=""><img file="EP0280489A2_D0002.tif" /></maths> where: P164 = Pressure at Location 164     P166 = Pressure at Location 166     f = Friction Factor (from Moody chart)     L = Pipe Length     D = Pipe Diameter     V = Fluid Velocity also: f = f(Re, <img file="EP0280489A2_D0003.tif" />) and Re = <img file="EP0280489A2_D0004.tif" /> where e/D is a measure of pipe surface roughness and is taken as 0.0002 in this example. For water at 20°C, ρ = 999Kf/m³ and µ = 1.0 x 10⁻³ Kf/m(sec). Assume flow rate Q = 0.1m³/sec, so:<maths id="math0003" num=""><img file="EP0280489A2_D0005.tif" /></maths> For Re = 8.48 x 10⁵ and e/D = 0.0002, the Moody chart, which indicates a pipe friction factor for a surface roughness (e/D) and Reynolds number (Re), a dimensionless constant that is a measure of flow turbulence, gives f = 0.0149, then,<maths id="math0004" num=""><img file="EP0280489A2_D0006.tif" /></maths> P₁₆₄ - P₁₆₆ = 15.9 KPa (kilo Pascals)      (6) Since P₁₆₄ would be monitored, the value of P₁₆₆ could be calculated using the above equation. Similar relations exist for heat and fluid flow transfer. References for fluid mechanics and heat transfer local to global transfor­mations that can be used to produce transformations for more complex structures include: Fox and McDonald, <u style="single">Intro­duction</u><u style="single">to</u><u style="single">Fluid</u><u style="single">Mechanics</u>, Second Edition, John Wiley and Sons, New York, 1978; and Ozisik, M. Necati, <u style="single">Basic</u><u style="single">Heat</u><u style="single">Transfer</u>, McGraw Hill, New York 1977. One of ordinary skill in the art could provide suitable transformation equations from the teachings of these references.

0033After the completion of the transformation 290 (Fig. 9), all critical locations have all the process, control, and response parameters required for fatigue calculation either from direct monitoring or extrapolation. At this point, stresses and fatigue are calculated for each component desired (by selecting a component 292), each of which may have multiple critical locations (and by select­ ing a location 294) where stress and fatigue information is desired. The process now described takes place for each critical location and each component.

0034First, stresses due to pressure loads are calcu­lated 296. This is done as follows. Pressure stresses are essentially linear with increasing pressure. Stresses for conservative pressure loads at critical locations are calculated in the design stress report prepared at plant design time by the plant design engineers. The value of this stress, and its associated pressure are scaled to a unit value and stored in the input database 60. Then, given the actual monitored pressure recorded and stored in the transient parameter database, the value of unit stress in the input database 60 is scaled accordingly. For example, suppose the pressure used in the design basis stress report is 20,685 kPa (3,000 psi). The stress at a critical location consists of 6 components, 3 normal stresses (σ_x σ_y σ_z) and 3 shear stresses (σ_xy σ_yz σ_zx). For simplicity, only one component is discussed since the methodology is identical for the others. Suppose the stress from the stress report is 82,740 kPa (12,000 psi) for the 20,685 kPa (3,000 psi) pressure at the location in question. That is, at a fluid pressure of 20,685 kPa, 82,740 kPa of stress is developed. The stress and associ­ated pressure that would be stored permanently in the input database 56 would be 27,580 kPa (4,000 psi) for a 6895 kPa (1,000 psi) load. Suppose now the cycle monitoring system records a pressure of 13,790 kPa (2,000 psi) at the loca­tion in question at a given time. The pressure stress computer 296 for that time would be 55,160 kPa (8,000 psi). Additional information on computing pressure stresses can be obtained from Shigley, <u style="single">Mechanical</u><u style="single">Engineering</u><u style="single">Design,</u> Third Edition, McGraw-Hill, New York, 1977.

0035Next, stresses due to general and thermal tran­sient stress are calculated 298. The method of calculation of transient thermal stresses will be discussed later with respect to Figs. 13-16. Stresses due to general thermal loads also exist. General thermal loads, to distinguish from the thermal transient loads, to distinguish from the thermal transient loads, are loads due to expansion or compression of materials due to, for example, global temperature changes in the temperature of fluid flowing through a pipe. General thermal loads are calculated at the time of plant design. The methodology to caclulate stresses for the general thermal loads is idential to that just described for pressure loads, except temperature is the monitored parameter. Additional information on general thermal stress calculations can be obtained from Roark and Young, <u style="single">Formulas</u><u style="single">for</u><u style="single">Stress</u><u style="single">and</u><u style="single">Strain</u>, Fifth Edition, McGraw-Hill, New York, 1975. Seismic loads and dead weight stress loads are fixed at the design time of the plant and stresses (general stresses) due to these loads are also calculated at plant design time and stored in the input database 60. Stresses due to these loads are read 300 from the input database 60. No scaling is necessary for these stress loads as is done for pressure and general thermal loads and the conservative design stresses are used. The stresses are combined 302 by simple addition of each stress component from each load condition.

0036The next step is to determine 304 the stress intensity ranges. This is done as follows. Each transient experienced by the plant and monitored by the system has one or more times where high or low peaks in the process parameters (temperature, pressure, etc.) appear in the transients. The stress peaks for pressure are coincident with the pressure peaks; however, temperature stress peaks lag behind the temperature transient peak. The procedure discussed hereinafter concerning the Duhamel integration produces the stress peaks. These are stored in the tran­sient parameter database 64. Once the stress peak(s) for the new transient are determined, the stress ranges are updated. For example, assume the present transient being considered has one peak stress set at the location in question (a spot on a pipe wall). Let the stresses for this transient be denoted by:     σx₁, σy₁, σz₁, σxy₁, σyz₁, σzx₁      (7) Assume also that two transients have occurred previously, each with one stress peak at the location desired. The stresses here are     σx₂, σy₂, σz₂, σxy₂, σyz₂, σzx₂      (8)     σx₃, σy₃, σz₃, σxy₃, σyz₃, σzx₃      (9) Note that all of these are peak stresses at a single point, the point of maximum stress through the wall. This point will not usually change with different transients experienced by nuclear power plants, and is usually the inside surface of the pipe, nozzle or component in question. We now make the notation more compact and denote the stresses from the three transients by:     (σ_ij)₁ (σ_ij)₂ and (σ_ij)₃      (10) and create three stress ranges:     (σ_ij)_A = (σ_ij)₁ - (σ_ij)₂      (11)     (σ_ij)_B = (σ_ij)₂ - (σ_ij)₃      (12)     (σ_ij)_C = (σ_ij)₃ - (σ_ij)₁      (13) Then, each of the above three sets of normal and shear stress ranges can be turned into principal stress ranges (σ₁, σ₂, σ₃) consisting only of normal stresses by solving the following cubic equation for each stress range above:     σ³ - Dσ² + Eσ - F = 0      (14) where     D = σ_x + σ_y + σ_z      (15)     E = σ_xσ_y + σ_yσ_z + σ_zσ_x - σ_xy² - σ_yz² - σ_zx²      (16)     F = σ_x σ_y σ_z + 2σ_xy σ_yz σ_zx - σ_xy²σ_z - σ_yz²σ_x - σ_zx²σ_y      (17) Solving the equation for each of the three stress ranges above results in the following principal stress ranges:     (σ₁ σ₂ σ₃)_A      (18)     (σ₁ σ₂ σ₃)_B      (19)     (σ₁ σ₂ σ₃)_C      (20) And the following stress intensity ranges     (σ₁₂ σ₂₃ σ₃₁)_A      (21)     (σ₁₂ σ₂₃ σ₃₁)_B      (22)     (σ₁₂ σ₂₃ σ₃₁)_C      (23) where, for ech of A, B, and C:     σ₁₂ = |σ₁ - σ₂|      (24)     σ₂₃ = |σ₂ - _σ3|      (25)     σ₃₁ = |σ₃ - σ₁|      (26) These are the stress intensity ranges 304 calculated for each transient. Additional information on stress peak determination, stress range and stress intensity range updating can be found in Shames, I.H., <u style="single">Introduction to Solid</u><u style="single">Mechanics</u>, Prentice Hall, Englewood Cliffs, N.J., 1975.

0037Next, the usage factor is determined 306 using the following procedure. The stress intensity ranges have been calculated for each transient peak as just described. For each transient peak, one stress intensity range is a maximum. Let us assume σ₁₂ is maximum for each transient. Each value of σ₁₂, that is (σ₁₂) _A (σ₁₂) _B, (σ₁₂) _C is compared to the fatigue curve for that location contained in the input database 60. The fatigue curve is determined at component design time and an example is illustrated in Fig. 12. The ordinate of the fatigue curve is stress intensity range and the abscissa is the allowable number of cycles on a logarithmic scale. Let us say that 100 cycles of (σ₁₂) _A are allowed, 50 cycles of (σ₁₂) _B are allowed, and 200 cycles of (σ₁₂) _C are allowed. Then, assume that one cycle of each stress intensity range occurs, the total fatigue usage at the point is:<maths id="math0005" num=""><img file="EP0280489A2_D0007.tif" /></maths> Where a value of 1.0 indicates end of component life. the new usage factor is used to update the output database 66 which is used for output and display. The above method of computation of stress intensity range 304 and fatigue usage 306 is typical and could be simplified depending on the situation at any given plant. Additional information on stress intensity range calculation and fatigue calculation can be found in the American Society of Mechanical Engineers Boiler and Pressure Vessel Code Section III, Subsection NB, 1983 Edition.

0038The stress amplitudes and usage factors, along with transient parameter information after global to local transformation, are stored in an output database 62 which can be accessed by personnel at the plant, using peripheral equipment 14, such as the cathode ray tube 42, printer 44, or plotter 46. The engineers can now use the information to make decisions regarding plant component maintenance, repair, replacement and life extension.

0039To determine transient thermal stress, the superposition principle is used which states that the stress acting on a component is the sum of the individual stresses acting on the component separately. For example, consider a beam subjected to two mechanical loads. At any location on the beam, the displacement, and indirectly the stress, is the sum of the displacements (stresses) caused by each of the loads separately. The same principle is being applied in the present invention except that temperature distributions are the load mechanisms and the location of the stress and the time are independent variables. For example, consider initially a plate with zero stress at an ambient temperature T. Assume that at the beginning of an analysis, a temperature change of Δ T, is applied to the plate for a period of time Δt₁. Stress will be induced in the plate as a result of the temperature change. If, in a subsequent time interval Δt₂, the plate is given another temperature change ΔT₂, additional stresses will be induced. The stress at the end of the second time interval t₂ = Δt₁ + Δt₂ is an accumulation (sum) of the first and second intervals. This procedure can be applied to any time period t_n where t_n = Δt₁ + Δt₂... Δt_n. When Δt_i approaches zero both the time and temperature change become continuous and the stress at any one time is given by the transfer function of equation 28 below. The temperature mentioned above is the externally applied temperature (the excitation temperature) that will cause stress changes in the plate. In this invention, the temperature is that of the fluid in contact with the surface of the mechanical component. If the temperature is essentially non-spatially dependent, as is the case in nuclear power fluid carrying devices, but is only time dependent, the method described above is sufficient to calculate the stress at any point on the plate. If a situation is encountered where the temperature is also spatially dependent, the spatial interaction can be taken into account by double summations of equation 28 with respect to both time and location.

0040The stress is determined in accordance with the following:<maths id="math0006" num=""><img file="EP0280489A2_D0008.tif" /></maths> where Q_ij is the stress response to the unit step function (f) and is called the kernel. Equation (28) is known as Duhamel's integral and additional discussion concerning this integral can be obtained in Hildebrand, <u style="single">Advanced Calculus for Applications</u>, Prentice-Hall, Englewood Cliff, N.J., 1976. As an alternative, Greens function, as described in Courant and Hilber, <u style="single">Methods</u><u style="single">of</u><u style="single">Mathematical</u><u style="single">Physics</u>, First Edition, Interscience Publishers, New Yor, 1937 can be used to calculate transient thermal stress. Additional information on calculations of thermal stress are: Carslow and Jaeger, <u style="single">Conduction</u><u style="single">of</u><u style="single">Heat</u><u style="single">in</u><u style="single">Solids</u>, 2d Ed. Oxford University Press, 1959; and Boley and Weiner, <u style="single">Theory</u><u style="single">of</u><u style="single">Thermal</u><u style="single">Stress</u>, John Wiley &amp; Sons, Inc., 1960.

0041The process for obtaining the stress at each location is illustrated in Fig. 13. The kernel data 350 is input from the input data base 60. The kernel data can be calculated at plant design time or later if the plant is already in service by doing a well-known finite element stress analysis on the mechanical component of interest. For example, assume that the kernel data in table 1 below is input for a boiler plate, as illustrated in Fig. 14 and the stress is desired at point A. <tables id="tabl0001" num="0001"><img file="EP0280489A2_D0009.tif" /></tables> The kernel data depicted in table 1 and illustrated as a curve in Fig. 15 is obtained at the component design and testing stage and represents the response of the component to a step temperature transient where the temperature T of the fluid at point A in contact with the plate is changed from T=0 to T=1. A finite element temperature/stress analysis is performed and the kernel function is the stress solution caused by the sudden temperature change. The kernel or transfer function takes into account not only the shape and dimensions of the component but also the properties of the material from which the component is made.

0042Next, the transient data and number of points is input 352 (Fig. 13) from the transient parameter database 64. For example, assume that the transient data, in Table 2 below, <tables id="tabl0002" num="0002"><img file="EP0280489A2_D0010.tif" /></tables> is input, where the transient is illustrated by Fig. 16. Then the time scale and data points from the transient and kernel must be correlated. That is, each point in the kernel must have a corresponding point in the transient. That is accomplished by interpolating between kernel points 354 and transient points 356 using the same time step t.<maths id="math0007" num=""><img file="EP0280489A2_D0011.tif" /></maths> where t = time, T = temperature and the subscripts 0, 1, 2 denote the data points of the input time-temperature data.

0043The output is a series of temperature data at the boundary points of each small segment along the time frame. For example, the temperature at time 5.0 seconds is 438°F, while the temperature of time 10.0 second is 375°F for the example of tables 1 and 2.

0044Next, the stress is calculated 358 for the transient using the following equation for the Duhamel integral:<maths id="math0008" num=""><img file="EP0280489A2_D0012.tif" /></maths> T equals temperature and Q equals the kernel function. Equation 30 integrates the stress response to a unit temperature input multiplied by the actual temperature input.

0045The calculations using the data of tables 1 and 2 would result in the stress set forth below in table 3. <tables id="tabl0003" num="0003"><img file="EP0280489A2_D0013.tif" /></tables> These stresses are then output to the next processing stage in which the component usage factor is determined.

0046The many features and advantages of the invention are apparent from the detailed specification and thus, it is intended by the appended claims to cover all such features and advantages of the invention which fall within the true spirit and scope thereof. Further, since numerous modifications and changes will readily occur to those skilled in the art, it is not desired to limit the invention to the exact construction and operation illustrated and described, and accordingly, all suitable modifications and equivalents may be resorted to, falling within the scope of the invention.