Method and system for controlling the rolling of copper alloy strips in a multi-acoustic field cooperation

Abstract

This invention provides a method and system for modulation control of copper alloy rolled material based on multi-sound field coordination. The field includes acquiring microstructure information, dividing the sound field into zones, calculating phase control parameters based on propagation time difference, establishing an inverse mapping relationship between microstructure density and sound field energy, collecting dislocation line orientation distribution, calculating the angle concentration index, and optimizing the phase difference and energy ratio between adjacent sound field excitation sources accordingly to achieve precise control of dislocation orientation in rolled material, thereby improving the rolling performance of copper alloy.

Description

Technical Field

[0001] This invention relates to metal rolling processing technology, and more particularly to a method and system for modulating and controlling copper alloy rolled materials with multiple acoustic fields. Background Technology

[0002] With the development of high-end equipment manufacturing and electronic information technology, higher requirements have been placed on the performance of copper alloy rolled materials. As a non-contact energy transfer method, sound field has shown unique advantages in the field of materials processing, providing a new technical approach for the control of copper alloy rolled materials.

[0003] Current acoustic field manipulation techniques for copper alloy rolled materials have significant shortcomings. Existing acoustic field manipulation methods mostly employ single-source excitation, failing to address differentiated manipulation based on the microstructure gradient differences along the thickness direction of the rolled material. This results in uneven acoustic energy distribution and unsatisfactory manipulation effects. Furthermore, existing technologies have not established a precise mapping relationship between dislocation orientation distribution and acoustic field parameters, making it difficult to achieve directional manipulation of dislocation orientation. This limits the performance improvement of rolled materials and fails to meet the precise control requirements of high-end applications for the anisotropic properties of copper alloy rolled materials. Summary of the Invention

[0004] This invention provides a method and system for modulating and controlling copper alloy rolled material with multiple sound fields in coordination, which can solve the problems in the prior art.

[0005] A first aspect of the present invention provides a method for modulating and controlling copper alloy rolled material with multi-field acoustic coordination, comprising:

[0006] Obtain information on the microstructure and target rolling performance of the copper alloy material to be rolled;

[0007] Based on the microstructure information, the calendered material is divided into multiple acoustic field zones along the thickness direction, and acoustic field excitation sources are arranged in each zone.

[0008] The propagation time of each sound field excitation source at different depth layers is detected, the propagation time difference is calculated, and the propagation time difference is converted into a phase modulation parameter.

[0009] Extract the microstructure density distribution curve of the calendered material, establish the inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic field excitation source, and calculate the energy distribution coefficient of each depth layer;

[0010] Based on the phase modulation parameters and energy distribution coefficient, each sound field excitation source is driven to apply a synergistic sound field effect to the calendered material. The orientation distribution of dislocation lines in the calendered material under the synergistic sound field effect is collected, the angle distribution between the dislocation lines and the calendering direction is statistically analyzed, and the angle concentration index is calculated.

[0011] The target enhancement coefficient for dislocation orientation is determined based on the target rolling performance index. The acoustic wave superposition gain is calculated based on the included angle concentration index. The acoustic wave superposition gain is compared with the target enhancement coefficient. Based on the comparison result, the phase difference and energy ratio between adjacent acoustic field excitation sources are optimized to achieve coordinated control of dislocation orientation in the rolled material.

[0012] Based on the microstructure information, the calendered material is divided into multiple acoustic field zones along the thickness direction, and acoustic field excitation sources are arranged in each zone, including:

[0013] Based on the microstructure information, extract the microstructure parameters of the calender at different thickness positions, calculate the microstructure difference value between adjacent thickness positions, determine the microstructure abrupt change position along the thickness direction of the calender based on the microstructure difference value, divide the calender into multiple acoustic field action zones with the microstructure abrupt change position as the boundary, and calculate the microstructure characteristic value of each acoustic field action zone.

[0014] The number of sound field excitation sources required in each sound field action zone is calculated based on the tissue characteristic values ​​and the size of the sound field action zone. A sound field superposition experiment is conducted on the sound field excitation sources to obtain the action radius of the sound field excitation sources. The layout spacing of the sound field excitation sources is determined according to the action radius and the number of sound field excitation sources.

[0015] Calculate the layout coordinates of the sound field excitation source and the corresponding sound field intensity distribution based on the layout spacing and the size of the sound field effect zone. Select the layout coordinates that satisfy the preset sound field coverage threshold, and place the sound field excitation source at the selected layout coordinates.

[0016] The propagation time of each sound field excitation source at different depth layers is detected, the propagation time difference is calculated, and the propagation time difference is converted into phase modulation parameters, including:

[0017] Acquire sound wave propagation signals from various sound field excitation sources at different depth layers, and perform time-frequency domain joint analysis on the sound wave propagation signals to extract the sound wave frequency drift and sound wave amplitude attenuation rate;

[0018] The dislocation density distribution of each depth layer is inverted by the acoustic frequency drift, and the grain boundary density distribution of each depth layer is inverted by the acoustic amplitude attenuation rate. The acoustic impedance correction value of each depth layer is calculated based on the dislocation density distribution and the grain boundary density distribution. The propagation time of each depth layer is corrected according to the acoustic impedance correction value to obtain the corrected propagation time.

[0019] The propagation time difference is obtained by calculating the difference between the corrected propagation times. The propagation time difference is then converted into the phase of the operating frequency of each sound field excitation source to obtain the excitation phase difference value. A spatiotemporal distribution map of the excitation phase difference value as the depth layer position is constructed.

[0020] Extract the spatial gradient vector field of the excitation phase difference from the spatiotemporal distribution map, identify the depth layer position with the largest gradient amplitude in the spatial gradient vector field, and take the excitation phase difference of the sound field excitation source corresponding to the depth layer position with the largest gradient amplitude as the reference phase difference.

[0021] The phase deviation between the excitation phase difference of other depth layer sound field excitation sources and the reference phase difference is calculated as the phase control parameter.

[0022] Extract the microstructure density distribution curve of the calendered material, establish the inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic field excitation source, and calculate the energy distribution coefficient of each depth layer, including:

[0023] Obtain the microstructure density values ​​of calendered material at different depths along the thickness direction, arrange the microstructure density values ​​according to the depth layer position, and construct a microstructure density distribution curve.

[0024] Calculate the rate of change of tissue density values ​​between adjacent depth layers, identify depth layers whose rate of change deviates from the mean, and mark them as areas of abnormal tissue density.

[0025] Read the tissue density values ​​of the abnormal tissue density area, establish the acoustic wave propagation impedance distribution, deduce the acoustic wave propagation path based on the acoustic wave propagation impedance distribution, and calculate the cumulative attenuation of acoustic wave energy.

[0026] The cumulative attenuation of acoustic energy is used as the incremental energy requirement of the acoustic field excitation source at each depth layer, and an inverse mapping relationship between the tissue density distribution curve and the output energy of the acoustic field excitation source is established.

[0027] Read the incremental output energy demand of the sound field excitation source at each depth layer, and construct the depth transmission link of the incremental output energy demand of the sound field excitation source;

[0028] Identify the anchor depth layer where abrupt changes occur in the depth transmission link, and perform bidirectional diffusion correction on the incremental output energy demand of the sound field excitation source in adjacent depth layers starting from the anchor depth layer. Normalize the corrected incremental output energy demand of the sound field excitation source to obtain the energy distribution coefficient of each depth layer.

[0029] Based on the phase modulation parameters and energy distribution coefficients, each acoustic excitation source is driven to apply a synergistic acoustic field to the calendered material. The orientation distribution of dislocation lines in the calendered material under the synergistic acoustic field is collected, and the angle distribution between the dislocation lines and the calendering direction is statistically analyzed. The angle concentration index is calculated, including:

[0030] The phase modulation parameters are converted into frequency control quantities of the sound field excitation source, and the energy distribution coefficient is converted into power control quantities of the sound field excitation source. Sound waves are generated based on the frequency control quantities and power control quantities.

[0031] The propagation time and pressure amplitude of sound waves in calendered material are collected, and the propagation time difference and pressure amplitude difference between adjacent sound field excitation sources are calculated to obtain the propagation characteristics of sound waves in the depth direction.

[0032] The beam pointing angle of the sound field excitation source is adjusted according to the propagation characteristics. The beam pointing angle is combined with the pressure amplitude difference to determine the sound wave superposition position and form a cooperative sound field. The stress state of the calendered material is collected in the area of ​​the cooperative sound field. The driving force of dislocation movement is calculated according to the stress state to determine the starting point of dislocation slip.

[0033] Record the position coordinates of the dislocation line from the starting point, construct the spatial trajectory of the dislocation line based on the position coordinates, extract the node information of the dislocation line, calculate the local direction of the dislocation line based on the node information, and project the local direction of the dislocation line onto the rolling direction to obtain the dislocation orientation distribution.

[0034] The included angle values ​​in the dislocation orientation distribution are statistically analyzed, the number distribution of included angle values ​​is calculated, and the number distribution is normalized to obtain the included angle concentration index.

[0035] The target enhancement coefficient for dislocation orientation is determined based on the target rolling performance index. The acoustic wave superposition gain is calculated based on the included angle concentration index. The acoustic wave superposition gain is compared with the target enhancement coefficient. Based on the comparison results, the phase difference and energy ratio between adjacent acoustic field excitation sources are optimized to achieve coordinated control of dislocation orientation in the rolled material. This includes:

[0036] Obtain the stress and texture requirements in the target calendering performance indicators, calculate the crystal orientation distribution and dislocation slip distribution of the calendered material, determine the critical shear force for dislocation motion based on the crystal orientation distribution and dislocation slip distribution, and convert the critical shear force into the target enhancement coefficient for dislocation orientation;

[0037] Acquire sound wave propagation signals from adjacent sound field excitation sources, calculate the sound wave propagation time difference and pressure amplitude difference, extract the angle value distribution from the angle concentration index, establish a dynamic response function between the angle value distribution and the sound wave propagation signal, adjust the output power of adjacent sound field excitation sources according to the dynamic response function, and convert the adjusted sound wave propagation signal into sound wave superposition gain.

[0038] The difference between the acoustic wave superposition gain and the target enhancement coefficient is used as the control deviation of dislocation orientation. The phase compensation value and power distribution ratio of adjacent acoustic field excitation sources are calculated based on the control deviation. The output parameters of the acoustic field excitation sources are adjusted based on the phase compensation value and power distribution ratio to achieve coordinated control of dislocation orientation in calendered material.

[0039] Establishing a dynamic response function between the included angle distribution and the sound wave propagation signal, and adjusting the output power of adjacent sound field excitation sources according to the dynamic response function includes:

[0040] Statistical analysis of the angle distribution data yields the peak position and distribution width of the angle values. The peak position is used as the dominant orientation parameter of dislocation orientation, and the distribution width is used as the dispersion parameter of dislocation orientation.

[0041] Acquire sound wave propagation signals from adjacent sound field excitation sources inside the calendered material, perform time-domain analysis on the sound wave propagation signals to extract the arrival time and peak amplitude of the sound waves, and calculate the sound wave propagation time difference and peak amplitude difference between adjacent sound field excitation sources.

[0042] The dominant direction parameter and the dispersion parameter are set as the output variables of the dynamic response function, and the sound wave propagation time difference and the sound wave peak amplitude difference are set as the input variables of the dynamic response function. The dynamic response function of the angle value distribution and the sound wave propagation signal is constructed by fitting the functional relationship between the input variables and the output variables.

[0043] Calculate the partial derivatives of the dynamic response function with respect to the sound wave propagation time difference and the sound wave peak amplitude difference. Based on the partial derivatives, convert the target values ​​of the sound wave propagation time difference and the sound wave peak amplitude difference into the output power adjustment amount of the adjacent sound field excitation source. Apply the output power adjustment amount to the corresponding sound field excitation source to complete the adjustment of the output power of the adjacent sound field excitation source.

[0044] A second aspect of the present invention provides a multi-sound-field coordinated modulation control system for copper alloy rolled materials, comprising:

[0045] The first unit is used to obtain information on the microstructure and target rolling performance of the copper alloy material to be rolled.

[0046] The second unit is used to divide the calendered material into multiple acoustic field action zones along the thickness direction according to the microstructure information, and to set up acoustic field excitation sources in each zone.

[0047] The third unit is used to detect the propagation time of each sound field excitation source at different depth layers, calculate the propagation time difference, and convert the propagation time difference into phase modulation parameters.

[0048] The fourth unit is used to extract the microstructure density distribution curve of the calendered material, establish the inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic field excitation source, and calculate the energy distribution coefficient of each depth layer.

[0049] The fifth unit is used to drive each sound field excitation source to apply a synergistic sound field effect to the calendered material according to the phase control parameters and energy distribution coefficient, collect the dislocation line orientation distribution of the calendered material under the synergistic sound field effect, statistically analyze the angle distribution between the dislocation lines and the calendering direction, and calculate the angle concentration index.

[0050] The sixth unit is used to determine the target enhancement coefficient of dislocation orientation based on the target rolling performance index, calculate the acoustic wave superposition gain according to the included angle concentration index, compare the acoustic wave superposition gain with the target enhancement coefficient, and optimize the phase difference and energy ratio between adjacent acoustic field excitation sources based on the comparison results, so as to realize the coordinated control of dislocation orientation of rolling material.

[0051] A third aspect of the present invention provides an electronic device, comprising:

[0052] processor;

[0053] Memory used to store processor-executable instructions;

[0054] The processor is configured to invoke instructions stored in the memory to execute the aforementioned method.

[0055] A fourth aspect of the present invention provides a computer-readable storage medium having stored thereon computer program instructions that, when executed by a processor, implement the aforementioned method.

[0056] In this embodiment, by performing zoned acoustic field control on the rolled material, the acoustic field can be precisely applied to different depth layers, thereby improving the utilization efficiency of acoustic field energy. The acoustic field energy of each depth layer can be accurately configured based on the inverse mapping relationship established according to the microstructure density distribution, ensuring the uniformity of the acoustic field effect. By using the concentration index of the angle between the dislocation orientation distribution and the rolling direction, combined with the dynamic optimization of the acoustic wave superposition gain, the dislocation orientation of the rolled material can be synergistically controlled, improving the microstructure uniformity and performance stability of the rolled material, and ultimately obtaining a high-performance copper alloy rolled product. Attached Figure Description

[0057] Figure 1 This is a flowchart illustrating the method for modulating and controlling copper alloy rolled material with multiple sound fields in accordance with an embodiment of the present invention.

[0058] Figure 2 This is a flowchart of the energy allocation correction algorithm according to an embodiment of the present invention. Detailed Implementation

[0059] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0060] The technical solution of the present invention will be described in detail below with reference to specific embodiments. These specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described again in some embodiments.

[0061] Figure 1 This is a flowchart illustrating the multi-sound-field coordinated modulation control method for copper alloy rolled material according to an embodiment of the present invention, as shown below. Figure 1 As shown, the method includes:

[0062] Obtain information on the microstructure and target rolling performance of the copper alloy material to be rolled;

[0063] Based on the microstructure information, the calendered material is divided into multiple acoustic field zones along the thickness direction, and acoustic field excitation sources are arranged in each zone.

[0064] The propagation time of each sound field excitation source at different depth layers is detected, the propagation time difference is calculated, and the propagation time difference is converted into a phase modulation parameter.

[0065] Extract the microstructure density distribution curve of the calendered material, establish the inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic field excitation source, and calculate the energy distribution coefficient of each depth layer;

[0066] Based on the phase modulation parameters and energy distribution coefficient, each sound field excitation source is driven to apply a synergistic sound field effect to the calendered material. The orientation distribution of dislocation lines in the calendered material under the synergistic sound field effect is collected, the angle distribution between the dislocation lines and the calendering direction is statistically analyzed, and the angle concentration index is calculated.

[0067] The target enhancement coefficient for dislocation orientation is determined based on the target rolling performance index. The acoustic wave superposition gain is calculated based on the included angle concentration index. The acoustic wave superposition gain is compared with the target enhancement coefficient. Based on the comparison result, the phase difference and energy ratio between adjacent acoustic field excitation sources are optimized to achieve coordinated control of dislocation orientation in the rolled material.

[0068] Based on the microstructure information, the calendered material is divided into multiple acoustic field zones along the thickness direction, and acoustic field excitation sources are arranged in each zone, including:

[0069] Based on the microstructure information, extract the microstructure parameters of the calender at different thickness positions, calculate the microstructure difference value between adjacent thickness positions, determine the microstructure abrupt change position along the thickness direction of the calender based on the microstructure difference value, divide the calender into multiple acoustic field action zones with the microstructure abrupt change position as the boundary, and calculate the microstructure characteristic value of each acoustic field action zone.

[0070] The number of sound field excitation sources required in each sound field action zone is calculated based on the tissue characteristic values ​​and the size of the sound field action zone. A sound field superposition experiment is conducted on the sound field excitation sources to obtain the action radius of the sound field excitation sources. The layout spacing of the sound field excitation sources is determined according to the action radius and the number of sound field excitation sources.

[0071] Calculate the layout coordinates of the sound field excitation source and the corresponding sound field intensity distribution based on the layout spacing and the size of the sound field effect zone. Select the layout coordinates that satisfy the preset sound field coverage threshold, and place the sound field excitation source at the selected layout coordinates.

[0072] Based on the microstructure information, microstructure parameters of the rolled material at different thickness locations are extracted, and the microstructure difference values ​​between adjacent thickness locations are calculated to determine the locations of microstructure abrupt changes and divide the acoustic field effect zones. Microstructure parameters such as grain size, grain boundary distribution, and second phase distribution are obtained for each thickness location. For copper alloy rolled materials, a measurement point is taken every 0.5 mm along the thickness direction, and the average grain size at that location is recorded. The difference in grain size between two adjacent measurement points is the microstructure difference value. When the microstructure difference value exceeds a preset threshold, it is determined to be a location of microstructure abrupt change. Taking a 10 mm thick brass rolled material as an example, 20 measurement points are set along the thickness direction. Measurements show that the microstructure difference values ​​at 2.5 mm and 7.5 mm are 5 μm and 6 μm, respectively, both exceeding the preset threshold of 4 μm. Therefore, these two locations are determined to be locations of microstructure abrupt changes, thus dividing the rolled material into three acoustic field effect zones: a surface zone (0-2.5 mm), an intermediate zone (2.5-7.5 mm), and a bottom zone (7.5-10 mm). Microstructural eigenvalues ​​were calculated for each partition. These eigenvalues ​​included a weighted average of parameters such as the average grain size, standard deviation of grain size, and grain boundary area ratio within the partition. The calculation results showed that the microstructural eigenvalues ​​for the three partitions were 0.68, 0.82, and 0.75, respectively.

[0073] The required number of acoustic excitation sources within each acoustic field zone is calculated based on the microstructure characteristics and the dimensions of the acoustic field zones. The number of acoustic excitation sources is related to the microstructure characteristics and the zone dimensions; a larger microstructure characteristic indicates a more non-uniform microstructure, requiring more acoustic excitation sources. A larger zone size also increases the required number of acoustic excitation sources. An acoustic field superposition experiment is conducted, placing multiple acoustic excitation sources on the surface of the rolled material, adjusting their distances, measuring the sound field intensity distribution, and determining the effective range of the acoustic field. The effective radius of the acoustic excitation source is defined as the distance at which the sound field intensity decays to 50% of its initial value. For the aforementioned brass rolled material, the effective radius of the ultrasonic transducer used is determined to be 25 mm through the acoustic field superposition experiment. The spacing of the acoustic excitation sources is determined based on the effective radius and the required number of acoustic excitation sources. The spacing should be less than twice the effective radius to ensure the continuity of the sound field coverage. For the surface, middle, and bottom zones, the required number of acoustic excitation sources is calculated to be 4, 6, and 5, respectively, based on their respective microstructure characteristics and zone dimensions. Considering the effective radius of 25 mm, the spacing between the installations is set to 40 mm to ensure that the sound fields of adjacent excitation sources have sufficient overlap.

[0074] The placement coordinates of the sound field excitation sources and the corresponding sound field intensity distribution are calculated based on the spacing between the sources and the dimensions of the sound field action zones. The placement coordinates of the sound field excitation sources are obtained through the zone dimensions and placement spacing. Taking a rectangular calendered material as an example, with a length of 500 mm and a width of 300 mm, four sound field excitation sources are placed in the surface zone with coordinates of (100, 75), (100, 225), (400, 75), and (400, 225), respectively. The coordinate unit is millimeters, and the origin is located at a corner of the calendered material. The sound field intensity distribution at each placement location is calculated, and the sound field intensity decreases with increasing distance from the excitation source. The sound field intensity distribution can be calculated using a sound field propagation model, considering parameters such as the material's acoustic impedance and absorption coefficient. Placement coordinates that satisfy a preset sound field coverage threshold are selected. The preset sound field coverage threshold is 30% of the initial sound field intensity, meaning that the sound field intensity at any location within the zone must not be less than 30% of the sound field intensity at the excitation source.

[0075] Calculations show that the four coordinates for the surface zone can meet the coverage threshold requirement for 95% of the area, the six coordinates for the middle zone can meet the requirement for 96%, and the five coordinates for the bottom zone can meet the requirement for 94%. After determining the coordinates, the sound field excitation source is placed at the selected location. For different zones, ultrasonic transducers with different frequencies and powers can be selected as the sound field excitation source: a transducer with a frequency of 20kHz and a power of 500W is selected for the surface zone, a transducer with a frequency of 25kHz and a power of 600W is selected for the middle zone, and a transducer with a frequency of 30kHz and a power of 700W is selected for the bottom zone.

[0076] This invention divides calendered material into multiple acoustic field zones based on the differences in microstructure along the thickness direction, and deploys appropriate acoustic field excitation sources for each zone's microstructure characteristics, achieving precise control of the acoustic field. Compared to traditional single acoustic field methods, this invention fully considers the inhomogeneity of the calendered material's internal microstructure. Through the synergistic effect of multiple acoustic fields, it significantly improves the utilization efficiency of acoustic field energy and enhances the modulation effect of the acoustic field on the material's microstructure.

[0077] In one optional implementation, detecting the propagation time of each sound field excitation source at different depth layers, calculating the propagation time difference, and converting the propagation time difference into a phase modulation parameter includes:

[0078] Acquire sound wave propagation signals from various sound field excitation sources at different depth layers, and perform time-frequency domain joint analysis on the sound wave propagation signals to extract the sound wave frequency drift and sound wave amplitude attenuation rate;

[0079] The dislocation density distribution of each depth layer is inverted by the acoustic frequency drift, and the grain boundary density distribution of each depth layer is inverted by the acoustic amplitude attenuation rate. The acoustic impedance correction value of each depth layer is calculated based on the dislocation density distribution and the grain boundary density distribution. The propagation time of each depth layer is corrected according to the acoustic impedance correction value to obtain the corrected propagation time.

[0080] The propagation time difference is obtained by calculating the difference between the corrected propagation times. The propagation time difference is then converted into the phase of the operating frequency of each sound field excitation source to obtain the excitation phase difference value. A spatiotemporal distribution map of the excitation phase difference value as the depth layer position is constructed.

[0081] Extract the spatial gradient vector field of the excitation phase difference from the spatiotemporal distribution map, identify the depth layer position with the largest gradient amplitude in the spatial gradient vector field, and take the excitation phase difference of the sound field excitation source corresponding to the depth layer position with the largest gradient amplitude as the reference phase difference.

[0082] The phase deviation between the excitation phase difference of other depth layer sound field excitation sources and the reference phase difference is calculated as the phase control parameter.

[0083] In this embodiment, sound wave propagation signals from various sound field excitation sources at different depth layers are collected. Time-frequency domain joint analysis is performed on the sound wave propagation signals to extract the sound wave frequency drift and amplitude attenuation rate. During the acquisition process, multiple ultrasonic transducers are arranged on the surface of the copper alloy rolled material as sound field excitation sources, and sound wave receiving probes are arranged at different depths inside the rolled material. Taking a 10 mm thick brass rolled material as an example, a receiving point is set every 1 mm in the depth direction, for a total of 10 depth layers. The pulse-echo method is used to measure the sound wave propagation signal at each depth layer. An ultrasonic flaw detector is used to collect ultrasonic echo signals, with a sampling frequency set to 100 MHz and a sampling duration of 100 microseconds. After preprocessing, the collected sound wave propagation signals are subjected to time-frequency domain joint analysis using short-time Fourier transform to extract the sound wave frequency drift and amplitude attenuation rate. In the time-frequency spectrum analysis, a Hanning window is selected as the window function, with a window length of 512 points and an overlap rate of 50%. By analyzing the changes in the peak value of the dominant frequency in the time-frequency spectrum, the frequency drift of the sound wave at each depth layer was extracted. For a sound wave with an initial frequency of 20kHz, the frequency drift measured at a depth of 5mm was -320Hz. The sound wave amplitude attenuation rate was obtained by calculating the natural logarithm of the ratio of the amplitude of the echo signal at each depth layer to the amplitude of the initial transmitted signal. Similarly, the amplitude attenuation rate measured at a depth of 5mm was 0.42dB / mm.

[0084] The dislocation density distribution at each depth layer is inverted by the acoustic wave frequency drift, and the grain boundary density distribution at each depth layer is inverted by the acoustic wave amplitude attenuation rate. Based on the dislocation density distribution and grain boundary density distribution, the acoustic impedance correction value for each depth layer is calculated. The propagation time for each depth layer is then corrected according to the acoustic impedance correction value to obtain the corrected propagation time. The acoustic wave frequency drift is affected by dislocations in the material's microstructure, and there is a corresponding relationship between the frequency drift and dislocation density. According to the physical principle of sound wave propagation in metallic materials, the frequency drift is approximately proportional to the dislocation density. The frequency drift can be converted into dislocation density using a calibration coefficient. For the measured brass rolled material, the dislocation density corresponding to a frequency drift of -320Hz is approximately 8 × 10⁻⁶. 10 / cm 2 The amplitude attenuation of acoustic waves is mainly affected by grain boundary scattering, and the amplitude attenuation rate is closely related to the grain boundary density. By establishing the correspondence between the amplitude attenuation rate and the grain boundary density, the grain boundary density distribution at each depth layer can be derived. For a measured amplitude attenuation rate of 0.42 dB / mm, the corresponding grain boundary density is approximately 4 × 10⁻⁶. 4 / mm 2Both dislocation density and grain boundary density affect the acoustic impedance of a material. By considering both parameters, a correction value for the acoustic impedance can be calculated. In the calculation, dislocation density and grain boundary density are weighted at 0.6 and 0.4 respectively, resulting in a correction value of 0.87 for the acoustic impedance at a depth of 5 mm. Multiplying the original propagation time by the correction value yields the corrected propagation time. For a depth layer with an original propagation time of 6.5 microseconds, the corrected propagation time is 5.66 microseconds.

[0085] The propagation time difference is obtained by calculating the difference between the corrected propagation times. The propagation time difference is then converted to the operating frequency of each sound field excitation source to obtain the excitation phase difference value. A spatiotemporal distribution map of the excitation phase difference value as a function of depth layer position is constructed. The difference in corrected propagation time between adjacent depth layers is the propagation time difference. For two depth layers with depths of 4 mm and 5 mm, the corrected propagation times are 4.53 microseconds and 5.66 microseconds, respectively, with a propagation time difference of 1.13 microseconds. The propagation time difference can be converted to a phase difference using the operating frequency. The formula is: phase difference = propagation time difference multiplied by operating frequency multiplied by 360 degrees. For a sound field excitation source with an operating frequency of 20 kHz, a propagation time difference of 1.13 microseconds is converted to a phase difference of 8.14 degrees. Data from all depth layers of the calendered material are processed to obtain the excitation phase difference values ​​for different depth layers, and a spatiotemporal distribution map of the excitation phase difference value as a function of depth is constructed. The spatiotemporal distribution map is represented in heatmap form, with the horizontal axis representing depth position and the vertical axis representing time.

[0086] The spatial gradient vector field of the excitation phase difference is extracted from the spatiotemporal distribution map. The depth layer position with the largest gradient amplitude in the spatial gradient vector field is identified, and the excitation phase difference value of the acoustic excitation source corresponding to the depth layer position with the largest gradient amplitude is used as the reference phase difference value. The spatial gradient vector field represents the rate of change of the excitation phase difference value in space, which can be obtained by calculating the difference of the phase difference values ​​of adjacent depth layers. For the spatiotemporal distribution map of brass rolled material, the calculated spatial gradient vector field shows that the gradient amplitude is the largest at a depth of 7 mm, reaching 4.3 degrees / mm, and the phase difference value at this point is 12.5 degrees, which is set as the reference phase difference value.

[0087] The phase deviation between the excitation phase difference of each sound field excitation source at other depth layers and the reference phase difference is calculated as a phase control parameter. The difference between the phase difference of each depth layer and the reference phase difference is the phase deviation, which serves as the phase control parameter for the sound field excitation source. For a depth of 2 mm, the phase difference is 3.6 degrees and the phase deviation is -8.9 degrees; for a depth of 5 mm, the phase difference is 8.14 degrees and the phase deviation is -4.36 degrees. These phase control parameters are used to adjust the initial phase of each sound field excitation source so that the sound field can form enhanced interference at the target depth layer, thereby achieving precise modulation of a specific depth layer. In practice, a digital signal processor controls the phase of the driving signal of each sound field excitation source to achieve precise phase control.

[0088] In this embodiment, the synergistic effect of multiple sound fields within copper alloy rolled material can be effectively improved. By accurately detecting the propagation characteristics of sound waves at different depths and combining this with material microstructure parameters for propagation time correction, precise focusing of sound field energy at the target depth layer is achieved. This phase-modulated multi-sound field synergistic method overcomes the limitation of insufficient depth-direction precision in traditional sound field control methods. It can achieve directional energy transfer to specific depth layers based on the distribution characteristics of the rolled material's internal structure, thereby precisely controlling the material's microstructure evolution process at the microscale.

[0089] like Figure 2 The diagram illustrates the energy allocation correction algorithm flow of this embodiment.

[0090] In one optional implementation, the microstructure density distribution curve of the calendered material is extracted, an inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic field excitation source is established, and the energy distribution coefficient of each depth layer is calculated, including:

[0091] Obtain the microstructure density values ​​of calendered material at different depths along the thickness direction, arrange the microstructure density values ​​according to the depth layer position, and construct a microstructure density distribution curve.

[0092] Calculate the rate of change of tissue density values ​​between adjacent depth layers, identify depth layers whose rate of change deviates from the mean, and mark them as areas of abnormal tissue density.

[0093] Read the tissue density values ​​of the abnormal tissue density area, establish the acoustic wave propagation impedance distribution, deduce the acoustic wave propagation path based on the acoustic wave propagation impedance distribution, and calculate the cumulative attenuation of acoustic wave energy.

[0094] The cumulative attenuation of acoustic energy is used as the incremental energy requirement of the acoustic field excitation source at each depth layer, and an inverse mapping relationship between the tissue density distribution curve and the output energy of the acoustic field excitation source is established.

[0095] Read the incremental output energy demand of the sound field excitation source at each depth layer, and construct the depth transmission link of the incremental output energy demand of the sound field excitation source;

[0096] Identify the anchor depth layer where abrupt changes occur in the depth transmission link, and perform bidirectional diffusion correction on the incremental output energy demand of the sound field excitation source in adjacent depth layers starting from the anchor depth layer. Normalize the corrected incremental output energy demand of the sound field excitation source to obtain the energy distribution coefficient of each depth layer.

[0097] In this embodiment, the microstructure density values ​​of the rolled material at different depths along the thickness direction are obtained. These values ​​are then arranged according to depth to construct a microstructure density distribution curve. Microstructure density characterizes the distribution of discontinuities such as pores and defects within the material and is an important indicator for evaluating material microstructure quality. The microstructure density values ​​are obtained through a combination of metallographic microscopy observation and image analysis. For example, metallographic samples are taken from different depths of a 12 mm thick H62 brass rolled material. After grinding, polishing, and etching, microstructure photographs are taken under a metallographic microscope with a resolution of 2048 × 1536 pixels and a magnification of 200x. The photographs are binarized using image analysis software, and the ratio of the black area (pores, defects, etc.) to the entire field of view is calculated. Subtracting this ratio from 1 gives the microstructure density. Samples are taken from the rolled material every 1 mm from the surface to the bottom, yielding microstructure density values ​​for 12 depths. The tissue density was found to be 0.985 for the surface layer (0-1 mm), 0.962 for the intermediate layer (5-6 mm), and 0.978 for the bottom layer (11-12 mm). These data points were arranged by depth, and a continuous tissue density distribution curve was generated using cubic spline interpolation. The x-axis of the curve represents depth, and the y-axis represents tissue density.

[0098] The rate of change of microstructure density values ​​between adjacent depth layers was calculated, and depth layers with rates of change deviating from the mean were identified and marked as areas of abnormal microstructure density. The rate of change of microstructure density between adjacent depth layers was calculated by dividing the difference in microstructure density between two adjacent points by the depth interval. For H62 brass rolled material, 11 rate of change values ​​were calculated, with an average value of 0.002 / mm. A rate of change deviating from the mean by more than twice the standard deviation was defined as an outlier; the calculated standard deviation was 0.003 / mm. Therefore, depth layers with an absolute rate of change greater than 0.008 / mm were marked as areas of abnormal microstructure density. Analysis showed that the rate of change was -0.009 / mm at a depth of 3-4 mm and 0.011 / mm at a depth of 8-9 mm; these two areas were marked as areas of abnormal microstructure density.

[0099] The tissue density values ​​of areas with abnormal tissue density are read to establish the acoustic wave propagation impedance distribution. The acoustic wave propagation path is deduced based on this distribution, and the cumulative attenuation of acoustic wave energy is calculated. Acoustic wave propagation impedance is directly related to the tissue density of the material, and the two are approximately linearly related. The formula for calculating acoustic wave propagation impedance is: acoustic wave propagation impedance equals the material density multiplied by the speed of sound propagation in the material. For H62 brass, the basic acoustic wave propagation impedance is approximately 30 × 10⁻⁶. 6 kg / m 2 Based on this, adjustments are made according to changes in tissue density. For every 0.01 decrease in tissue density, the acoustic wave propagation impedance increases by approximately 0.5 × 10⁻⁶. 6 kg / m 2 Based on this relationship, the acoustic wave propagation impedance at a depth of 3-4 mm is calculated to be 30.9 × 10⁻⁶. 6 kg / m 2 The acoustic wave propagation impedance at a depth of 8-9 mm is 29.7 × 10⁻⁶. 6 kg / m 2 During sound wave propagation, reflection and transmission occur due to changes in acoustic impedance, resulting in energy attenuation. The amount of sound wave energy attenuation is proportional to the square of the difference in acoustic impedance. Calculations using a sound propagation model show that the cumulative attenuation of sound wave energy is approximately 2.1 dB from the surface to a depth of 3-4 mm; and approximately 4.3 dB from the surface to a depth of 8-9 mm.

[0100] The cumulative attenuation of acoustic energy is used as the incremental energy requirement of the acoustic excitation source at each depth layer, establishing an inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic excitation source. The incremental energy requirement of the acoustic excitation source is proportional to the cumulative attenuation of acoustic energy. To compensate for the attenuation of acoustic waves during propagation, the output energy of the acoustic excitation source needs to be adjusted accordingly. The inverse mapping relationship is established through a fitting function. For H62 brass rolled material, an exponential function is used: the incremental energy requirement of the acoustic excitation source is equal to the base energy multiplied by an exponential function of the cumulative attenuation of acoustic energy. The base energy is set to 100W, and the exponential coefficient is 0.15. Based on this relationship, the incremental energy requirement of the acoustic excitation source for the 3-4 mm depth layer is calculated to be 36.7W, and the incremental energy requirement for the 8-9 mm depth layer is 90.2W.

[0101] The energy increment demand of the sound field excitation source at each depth layer is read, and a depth transmission link for the energy increment demand of the sound field excitation source is constructed. A depth transmission link is a data structure representing the transmission relationship of energy demand between different depth layers. The construction method is as follows: the energy increment demand of each depth layer is arranged in depth order to form a one-dimensional array, and then the transmission coefficient is calculated through the difference between adjacent elements. For the 12 depth layers of H62 brass rolled material, the constructed depth transmission link contains 12 nodes, each node recording the energy increment demand of that depth layer and the transmission coefficient with adjacent depth layers.

[0102] Anchor depth layers with abrupt changes in the depth transmission link are identified. Starting from these anchor depth layers, a bidirectional diffusion correction is performed on the incremental output energy demand of the sound field excitation source in adjacent depth layers. The corrected incremental output energy demand is then normalized to obtain the energy distribution coefficient for each depth layer. Nodes in the depth transmission link where the transmission coefficient changes by more than 50% between adjacent nodes are identified as abrupt change points, and the corresponding depth layers are designated as anchor depth layers. For H62 brass rolled material, two anchor depth layers, 3-4 mm and 8-9 mm, are identified. Starting from these anchor depth layers, a bidirectional diffusion correction is performed on the incremental energy demand of adjacent depth layers. The diffusion correction uses an exponential decay model, with the correction magnitude decreasing as distance increases. After correction, the incremental energy demand of the 12 depth layers is normalized so that their sum equals 1, yielding the energy distribution coefficient for each depth layer. After normalization, the energy distribution coefficient of the 3-4 mm depth layer is 0.12, the energy distribution coefficient of the 8-9 mm depth layer is 0.19, and the energy distribution coefficient of the other depth layers is between 0.05 and 0.10.

[0103] In this embodiment, the modulation control process of copper alloy rolled material under the synergistic effect of multiple sound fields can be effectively guided. The non-uniformity of the internal structure of the rolled material is fully considered, and precise energy compensation is performed for areas with abnormal density, overcoming the inherent defect of sound wave energy attenuation with depth in traditional sound field processing methods. By precisely distributing the sound field energy along the depth direction, a balanced effect of sound field energy is achieved in each depth layer of the rolled material, significantly improving the uniformity and effectiveness of the sound field processing.

[0104] In one optional implementation, each acoustic excitation source is driven to apply a synergistic acoustic field to the calendered material according to the phase modulation parameters and energy distribution coefficient. The dislocation line orientation distribution of the calendered material under the synergistic acoustic field is collected, the angle distribution between the dislocation lines and the calendering direction is statistically analyzed, and the angle concentration index is calculated, including:

[0105] The phase modulation parameters are converted into frequency control quantities of the sound field excitation source, and the energy distribution coefficient is converted into power control quantities of the sound field excitation source. Sound waves are generated based on the frequency control quantities and power control quantities.

[0106] The propagation time and pressure amplitude of sound waves in calendered material are collected, and the propagation time difference and pressure amplitude difference between adjacent sound field excitation sources are calculated to obtain the propagation characteristics of sound waves in the depth direction.

[0107] The beam pointing angle of the sound field excitation source is adjusted according to the propagation characteristics. The beam pointing angle is combined with the pressure amplitude difference to determine the sound wave superposition position and form a cooperative sound field. The stress state of the calendered material is collected in the area of ​​the cooperative sound field. The driving force of dislocation movement is calculated according to the stress state to determine the starting point of dislocation slip.

[0108] Record the position coordinates of the dislocation line from the starting point, construct the spatial trajectory of the dislocation line based on the position coordinates, extract the node information of the dislocation line, calculate the local direction of the dislocation line based on the node information, and project the local direction of the dislocation line onto the rolling direction to obtain the dislocation orientation distribution.

[0109] The included angle values ​​in the dislocation orientation distribution are statistically analyzed, the number distribution of included angle values ​​is calculated, and the number distribution is normalized to obtain the included angle concentration index.

[0110] In this embodiment, the phase modulation parameter is converted into a frequency control quantity for the sound field excitation source, and the energy distribution coefficient is converted into a power control quantity for the sound field excitation source. Sound waves are generated based on the frequency and power control quantities. The phase modulation parameter is converted into a frequency control quantity for the sound field excitation source through frequency modulation. For a depth layer with a phase deviation of -8.9 degrees, the frequency control quantity is calculated as the base frequency plus the frequency offset. The frequency offset is equal to the base frequency multiplied by the phase deviation divided by 360 degrees. Taking a base frequency of 20kHz as an example, the frequency offset corresponding to a phase deviation of -8.9 degrees is -0.49kHz, therefore the frequency control quantity for the sound field excitation source at this depth layer is 19.51kHz. The energy distribution coefficient is directly converted into a power control quantity for the sound field excitation source. For a depth layer with an energy distribution coefficient of 0.12, its power control quantity is equal to the total system power multiplied by this coefficient. Taking a total system power of 1000W as an example, the power control quantity for the sound field excitation source at this depth layer is 120W. Based on the calculated frequency and power control values, a control signal is generated by a digital signal processor to drive the ultrasonic power amplifier to output the corresponding excitation signal, which is ultimately converted into sound waves by the ultrasonic transducer. For H62 brass rolled material, eight acoustic field excitation sources are set up, each corresponding to a different depth layer, and each is independently controlled according to its own frequency and power control values.

[0111] The propagation time and pressure amplitude of sound waves in the rolled material were collected, and the propagation time difference and pressure amplitude difference between adjacent sound field excitation sources were calculated to obtain the propagation characteristics of sound waves in the depth direction. An ultrasonic receiving probe array was arranged on the surface of the rolled material to receive sound wave signals from various depth layers. The probe array consisted of 16 piezoelectric sensors, evenly distributed within a 100mm × 100mm area. The collected sound wave signals included two key parameters: propagation time and pressure amplitude. Propagation time was obtained by measuring the time interval from sound wave emission to reception, and pressure amplitude was obtained by measuring the peak voltage of the received signal. For H62 brass rolled material, the sound wave propagation time at a depth of 3mm was 2.16 microseconds, and the pressure amplitude was 0.85MPa; the sound wave propagation time at a depth of 4mm was 2.87 microseconds, and the pressure amplitude was 0.72MPa. The propagation time difference and pressure amplitude difference between adjacent depth layers were calculated; the propagation time difference between the 3-4mm depth layers was 0.71 microseconds, and the pressure amplitude difference was 0.13MPa. By analyzing data from all depth layers, the propagation characteristic curves of sound waves in the depth direction are obtained, characterizing the variation of sound wave velocity and attenuation with depth.

[0112] The beam pointing angle of the sound field excitation source is adjusted according to the propagation characteristics. The beam pointing angle is combined with the pressure amplitude difference to determine the sound wave superposition position, forming a cooperative sound field. The stress state of the rolled material is collected in the cooperative sound field area, and the driving force for dislocation motion is calculated based on the stress state to determine the starting point of dislocation slip. The beam pointing angle of the sound field excitation source is adjusted using phased array technology. Based on the propagation characteristics of sound waves in the depth direction, the required beam pointing angle for each depth layer is calculated. For depth layers requiring sound field enhancement, the beam pointing angle is adjusted to the direction that maximizes sound wave energy. In H62 brass rolled material, the beam pointing angle is adjusted to 5.3 degrees for a 3-4mm depth layer and 7.8 degrees for an 8-9mm depth layer. The optimal sound wave superposition position is determined by combining the beam pointing angle with the pressure amplitude difference. The sound wave superposition position is selected in the region with the largest pressure amplitude difference and a suitable beam pointing angle to achieve maximum focusing of sound field energy. Within the defined region of synergistic acoustic field interaction, in-situ transmission electron microscopy was used to observe the microstructural changes of the rolled material and collect stress state data. The stress state was obtained by measuring the degree of lattice distortion; a higher degree of lattice distortion indicates higher local stress. Based on the stress state data, the driving force for dislocation motion was calculated. The driving force for dislocation motion is equal to the applied shear stress minus the critical shear stress. When the driving force is greater than zero, dislocation slip begins. For H62 brass rolled material, the shear stress measured at a depth of 4 mm was 58 MPa, the critical shear stress was 45 MPa, and the calculated driving force for dislocation motion was 13 MPa. This location was determined as the starting point of dislocation slip.

[0113] The system records the position coordinates of dislocation lines from the starting point, constructs their spatial trajectory based on these coordinates, extracts nodal information, calculates their local orientation, and projects this local orientation onto the rolling direction to obtain the dislocation orientation distribution. Using a high-resolution transmission electron microscope, position coordinates are recorded every 50 nm along the dislocation line's extension direction, starting from the dislocation slip initiation point, forming a spatial trajectory dataset. For a typical dislocation line in H62 brass rolled material, 20 position coordinate points are recorded, forming a complete spatial trajectory. A three-dimensional spatial curve is constructed based on these coordinates, and nodal information is extracted. This nodal information includes the nodal position, curvature, and tangent direction. For the recorded dislocation line trajectory, the plane determined by three adjacent coordinate points is calculated; the normal vector of this plane is the local orientation of the dislocation line at that node. The local orientation of the dislocation line is projected onto the rolling direction, and the angle between the two is calculated to obtain the dislocation orientation distribution. For the observed dislocation lines, the angle between the local direction and the rolling direction is distributed between 15 degrees and 75 degrees, with an average angle of 42 degrees.

[0114] The included angle values ​​in the dislocation orientation distribution were statistically analyzed, and the number distribution of these included angle values ​​was calculated. The distribution was then normalized to obtain the included angle concentration index. The angle between the dislocation line and the rolling direction was divided into nine intervals: 0-10 degrees, 10-20 degrees, ..., 80-90 degrees. The number of dislocation line nodes within each interval was counted, forming a histogram of the included angle value distribution. For H62 brass rolled material, the included angle distribution exhibited a bimodal characteristic, mainly concentrated in the 30-40 degree and 60-70 degree intervals, accounting for 28% and 23% of the total number of nodes, respectively. The distribution was normalized, and the included angle concentration index was calculated. The included angle concentration index is defined as the proportion of nodes within the dominant included angle interval to the total number of nodes. The 30-40 degree interval was determined as the dominant included angle interval, with an included angle concentration index of 0.28. A higher angle concentration index indicates a more concentrated dislocation line orientation and stronger anisotropy in the material microstructure; a lower index indicates a more uniform distribution of dislocation line orientation and better isotropy in the material microstructure. The change in the angle concentration index before and after the synergistic acoustic field treatment can be used to evaluate the effectiveness of acoustic field treatment in controlling the microstructure of calendered materials.

[0115] In this embodiment, precise control over the microstructure of calendered materials can be achieved, particularly the effective adjustment of the orientation distribution of dislocation lines. Through precise control of phase and energy, a synergistic sound field is formed along the depth direction, effectively overcoming the defect of uneven energy attenuation in the depth direction of a single sound field. Under the action of the sound field, the slip and rearrangement of dislocation lines are directionally controlled, causing the dislocation line orientation to tend towards a reasonable distribution, effectively reducing stress concentration inside the calendered material, and improving the material's microstructure uniformity and performance stability.

[0116] In one optional implementation, the target enhancement coefficient for dislocation orientation is determined based on the target rolling performance index, the acoustic wave superposition gain is calculated based on the included angle concentration index, the acoustic wave superposition gain is compared with the target enhancement coefficient, and the phase difference and energy ratio between adjacent acoustic field excitation sources are optimized based on the comparison result to achieve coordinated control of dislocation orientation in the rolled material, including:

[0117] Obtain the stress and texture requirements in the target calendering performance indicators, calculate the crystal orientation distribution and dislocation slip distribution of the calendered material, determine the critical shear force for dislocation motion based on the crystal orientation distribution and dislocation slip distribution, and convert the critical shear force into the target enhancement coefficient for dislocation orientation;

[0118] Acquire sound wave propagation signals from adjacent sound field excitation sources, calculate the sound wave propagation time difference and pressure amplitude difference, extract the angle value distribution from the angle concentration index, establish a dynamic response function between the angle value distribution and the sound wave propagation signal, adjust the output power of adjacent sound field excitation sources according to the dynamic response function, and convert the adjusted sound wave propagation signal into sound wave superposition gain.

[0119] The difference between the acoustic wave superposition gain and the target enhancement coefficient is used as the control deviation of dislocation orientation. The phase compensation value and power distribution ratio of adjacent acoustic field excitation sources are calculated based on the control deviation. The output parameters of the acoustic field excitation sources are adjusted based on the phase compensation value and power distribution ratio to achieve coordinated control of dislocation orientation in calendered material.

[0120] During the rolling process, the stress and texture requirements in the target rolling performance indicators are obtained. The crystal orientation distribution and dislocation slip distribution of the rolled material are calculated. Based on the crystal orientation distribution and dislocation slip distribution, the critical shear force for dislocation motion is determined, and the critical shear force is converted into the target enhancement coefficient for dislocation orientation. The target rolling performance indicators are usually determined by the product application requirements, including mechanical properties such as tensile strength, yield strength, and elongation, as well as characteristic indicators such as conductivity and texture. For electronic C2680 brass rolled strip, the target tensile strength is 430 MPa, the yield strength is 340 MPa, the elongation is not less than 22%, the conductivity is not less than 28% IACS, and the texture requirement is that the cubic texture accounts for not less than 65%. The crystal orientation of the rolled material is analyzed by electron backscatter diffraction technology to obtain the crystal orientation distribution function. For C2680 brass rolled strip, the crystal orientation distribution is mainly concentrated in {110}. <112> and {112} <111> Two texture components, accounting for 45% and 25% respectively, were used. Possible slip systems were determined based on the crystal orientation distribution, and the dislocation slip distribution was calculated. The dislocation slip distribution refers to the distribution of dislocation density in different slip systems. Using a dislocation statistical model, the dislocation density in the primary slip system was found to be 8 × 10^10 / m², and in the secondary slip system, it was 3 × 10^10 / m². Based on the crystal orientation distribution and dislocation slip distribution, and combined with the Schmidt factor rule, the critical shear force for dislocation motion was calculated. The critical shear force represents the minimum stress required for a dislocation to begin moving. For C2680 brass rolled strip, the critical shear force of the primary slip system was 68 MPa. The critical shear force was converted into a target enhancement factor for dislocation orientation, and the conversion relationship was established based on the dislocation kinematic model. The target enhancement factor for dislocation orientation is defined as the increase in dislocation density in the primary slip system after the acoustic field is applied, relative to the original state. For a critical shear force of 68 MPa, the target enhancement factor for dislocation orientation is calculated to be 1.35, indicating that the dislocation density on the main slip system needs to be increased by 35%.

[0121] Acoustic wave propagation signals from adjacent acoustic field excitation sources were collected. The propagation time difference and pressure amplitude difference were calculated. The angle distribution in the angle concentration index was extracted, and a dynamic response function between the angle distribution and the acoustic wave propagation signal was established. Based on the dynamic response function, the output power of the adjacent acoustic field excitation sources was adjusted, and the adjusted acoustic wave propagation signal was converted into acoustic wave superposition gain. A piezoelectric sensor array was used to collect the acoustic wave propagation signals from adjacent acoustic field excitation sources. The sampling frequency was set to 100MHz, and the sampling duration was 200μs. For C2680 brass rolled strip, adjacent acoustic field excitation sources arranged at depths of 4mm and 5mm had acquired propagation times of 2.63μs and 3.28μs, respectively, with a calculated propagation time difference of 0.65μs. The acquired pressure amplitudes were 0.78MPa and 0.65MPa, respectively, with a calculated pressure amplitude difference of 0.13MPa. The orientation of dislocation lines was observed using a transmission electron microscope, and the angle distribution between the dislocation lines and the rolling direction was statistically analyzed. For C2680 brass rolled strip, the included angle values ​​are mainly concentrated in two ranges: 25°-35° and 55°-65°, with an included angle concentration index of 0.32, indicating that dislocations within the dominant included angle range account for 32% of the total number of dislocations. A dynamic response function is established for the included angle distribution and the sound wave propagation signal, describing the influence of the sound field on the dislocation orientation distribution. The dynamic response function adopts a nonlinear regression model, with the input variables being the sound wave propagation time difference and the pressure amplitude difference, and the output variable being the rate of change of the included angle concentration index. Based on the established dynamic response function, the output power of adjacent sound field excitation sources is adjusted. When the included angle concentration index is lower than the target value, the output power of the sound field excitation source in the depth layer corresponding to the dominant included angle range is increased; when the index is higher than the target value, the output power of the sound field excitation source in the depth layer corresponding to the non-dominant range is increased. For C2680 brass rolled strip, the output power of the acoustic excitation source for the 4mm depth layer was adjusted from the original 180W to 210W, and the output power for the 5mm depth layer was adjusted from 150W to 135W. The adjusted acoustic wave propagation signal was converted into acoustic wave superposition gain using an energy superposition model. The acoustic wave superposition gain is defined as the ratio of the adjusted acoustic field energy to the original acoustic field energy, and the calculated result is 1.28, indicating that the acoustic field energy was enhanced by 28%.

[0122] The difference between the acoustic wave superposition gain and the target enhancement coefficient is used as the control deviation for dislocation orientation. Based on this control deviation, the phase compensation value and power distribution ratio of adjacent acoustic field excitation sources are calculated. The output parameters of the acoustic field excitation sources are adjusted based on the phase compensation value and power distribution ratio to achieve coordinated control of dislocation orientation in the rolled material. The difference between the acoustic wave superposition gain (1.28) and the target enhancement coefficient (1.35) is 0.07; this difference represents the control deviation for dislocation orientation. A positive control deviation indicates insufficient current acoustic field intensity, requiring further enhancement; a negative deviation indicates excessive intensity, requiring appropriate reduction. The phase compensation value of adjacent acoustic field excitation sources is calculated based on the control deviation. The phase compensation value is proportional to the control deviation, and the proportionality coefficient is determined experimentally. For a control deviation of 0.07, the calculated phase compensation value is 12.6°. The phase compensation value is used to adjust the phase difference between adjacent acoustic field excitation sources to optimize the acoustic wave interference effect. The power distribution ratio is calculated based on the control deviation; the power distribution ratio is defined as the ratio of the output power of adjacent acoustic field excitation sources. For a control deviation of 0.07, the calculated power distribution ratio is 1.65, indicating that the output power of the acoustic excitation source for a 4mm depth layer should be 1.65 times that of a 5mm depth layer. Based on the phase compensation value and the power distribution ratio, the output parameters of the acoustic excitation source are adjusted by an ultrasonic control system. The phase compensation value of 12.6° is converted into a phase modulation control signal, which is used to adjust the phase through a digital phase shifter; the power distribution ratio of 1.65 is converted into a power modulation control signal, which is used to adjust the power through a programmable power amplifier. The acoustic excitation source is driven to work according to the adjusted parameters to form an optimized synergistic acoustic field. The changes in dislocation orientation distribution are monitored in real time using an electron microscope. When the angle concentration index reaches a set threshold, the acoustic field processing is completed. For C2680 brass rolled strip, after 20 minutes of synergistic acoustic field processing, the angle concentration index increased from the original 0.32 to 0.44. Dislocations are mainly concentrated in the range of 30°±5° with the rolling direction. The dislocation arrangement is more orderly, and the crystal orientation is more consistent, achieving synergistic control of dislocation orientation.

[0123] In this embodiment, by establishing a correspondence between rolling performance indicators and dislocation orientation distribution, macroscopic performance requirements are transformed into microstructure control targets, and the target enhancement coefficient is precisely matched through precise control of acoustic field parameters. Compared with traditional heat treatment or machining methods, this approach enables directional adjustment of the internal microstructure of the material without altering its chemical composition and surface state. Through the synergistic optimization of phase difference and energy ratio, the limitations of single acoustic field control are effectively overcome, significantly improving the utilization efficiency and control effect of acoustic field energy.

[0124] In one optional implementation, establishing a dynamic response function between the included angle distribution and the sound wave propagation signal, and adjusting the output power of adjacent sound field excitation sources according to the dynamic response function includes:

[0125] Statistical analysis of the angle distribution data yields the peak position and distribution width of the angle values. The peak position is used as the dominant orientation parameter of dislocation orientation, and the distribution width is used as the dispersion parameter of dislocation orientation.

[0126] Acquire sound wave propagation signals from adjacent sound field excitation sources inside the calendered material, perform time-domain analysis on the sound wave propagation signals to extract the arrival time and peak amplitude of the sound waves, and calculate the sound wave propagation time difference and peak amplitude difference between adjacent sound field excitation sources.

[0127] The dominant direction parameter and the dispersion parameter are set as the output variables of the dynamic response function, and the sound wave propagation time difference and the sound wave peak amplitude difference are set as the input variables of the dynamic response function. The dynamic response function of the angle value distribution and the sound wave propagation signal is constructed by fitting the functional relationship between the input variables and the output variables.

[0128] Calculate the partial derivatives of the dynamic response function with respect to the sound wave propagation time difference and the sound wave peak amplitude difference. Based on the partial derivatives, convert the target values ​​of the sound wave propagation time difference and the sound wave peak amplitude difference into the output power adjustment amount of the adjacent sound field excitation source. Apply the output power adjustment amount to the corresponding sound field excitation source to complete the adjustment of the output power of the adjacent sound field excitation source.

[0129] Statistical analysis of the angle distribution data yielded the peak position and distribution width of the angle values. The peak position was used as the dominant orientation parameter for dislocation orientation, and the distribution width as the dispersion parameter for dislocation orientation. Orientation data for 300 dislocation lines were collected from different regions of the calendered material, with 5 points measured for each dislocation line, resulting in a total of 1500 angle data points. The angle values ​​within the range of 0° to 90° were divided into 18 intervals, and the frequency of data within each interval was statistically analyzed and a distribution histogram was plotted. Kernel density estimation was used to smooth the distribution, and the peak position and distribution width of the angle distribution were extracted. The angle distribution exhibited a bimodal characteristic with a main peak at 32° and a secondary peak at 58°. The main peak distribution width was 7.5°, and the secondary peak distribution width was 10.2°, with the main peak accounting for 42% of the total data. The main peak position of 32° was determined as the dominant orientation parameter for dislocation orientation, and its distribution width of 7.5° was determined as the dispersion parameter for dislocation orientation. These two parameters comprehensively describe the orientation characteristics of dislocation lines in the calendered material.

[0130] Acoustic wave propagation signals from adjacent acoustic excitation sources within the calendered material were collected. Time-domain analysis was performed to extract the arrival time and peak amplitude of the sound waves, and the propagation time difference and peak amplitude difference between adjacent excitation sources were calculated. Eight piezoelectric sensors were arranged in a 5mm × 5mm grid on the surface of the calendered material to receive acoustic wave signals emitted from adjacent excitation sources at depths of 3mm and 6mm. The collected acoustic wave signals were filtered to remove environmental noise and high-frequency interference, and the signal envelope was extracted using Hilbert transform. The arrival time was defined as the moment when the acoustic wave amplitude exceeded three times the background noise, and the maximum value of the signal envelope was the peak amplitude. For a 3mm deep excitation source, the arrival time was 1.86μs, and the peak amplitude was 0.87MPa; for a 6mm deep excitation source, the arrival time was 3.72μs, and the peak amplitude was 0.63MPa. The calculated sound wave propagation time difference between adjacent sound field excitation sources is 1.86 μs, and the sound wave peak amplitude difference is 0.24 MPa. These parameters reflect the spatiotemporal characteristics of sound wave propagation inside the calendered material.

[0131] The dominant orientation parameter and the dispersion parameter were set as the output variables of the dynamic response function, while the sound wave propagation time difference and the sound wave peak amplitude difference were set as the input variables. The dynamic response function of the angle distribution and the sound wave propagation signal was constructed by fitting the functional relationship between the input and output variables. To establish the relationship between the input and output variables, 20 sound field processing experiments with different parameter combinations were designed. While keeping other conditions consistent, the frequency difference and power ratio of adjacent sound field excitation sources were varied to generate different sound wave propagation time differences (1.2 μs to 2.4 μs) and sound wave peak amplitude differences (0.15 MPa to 0.35 MPa). The dislocation orientation distribution of the samples after each experimental treatment was measured, and the changes in the dominant orientation parameter and the dispersion parameter were recorded. The dynamic response function was established using multiple regression analysis. Considering the possible nonlinear relationship between the variables, a quadratic polynomial model was used for fitting. The fitted dynamic response function shows that as the sound wave propagation time difference increases, the dominant direction parameter moves closer to the rolling direction (its value decreases), and the dispersion parameter decreases, indicating that dislocation orientations are more concentrated. Conversely, as the sound wave peak amplitude difference increases, the dominant direction parameter moves further away from the rolling direction (its value increases), and the dispersion parameter increases, indicating that the dislocation orientation distribution is more dispersed. The model fit goodness of fit reaches 0.89, indicating that the established dynamic response function can well describe the relationship between sound field parameters and dislocation orientation distribution.

[0132] The partial derivatives of the dynamic response function with respect to the sound wave propagation time difference and the sound wave peak amplitude difference are calculated. Based on these partial derivatives, the target values ​​of the sound wave propagation time difference and the sound wave peak amplitude difference are converted into output power adjustment values ​​for adjacent sound field excitation sources. These output power adjustment values ​​are then applied to the corresponding sound field excitation sources to adjust their output power. Based on the established dynamic response function, the partial derivatives of the function with respect to the two input variables at the current operating point (1.86 μs, 0.24 MPa) are calculated. The calculation results show that the partial derivative of the dominant orientation parameter with respect to the sound wave propagation time difference is -8.7° / μs, and the partial derivative with respect to the sound wave peak amplitude difference is 12.3° / MPa; the partial derivative of the dispersion parameter with respect to the sound wave propagation time difference is -3.2° / μs, and the partial derivative with respect to the sound wave peak amplitude difference is 5.6° / MPa. Based on the product performance requirements, the target parameters for dislocation orientation are determined as follows: the dominant orientation parameter is 25°, and the dispersion parameter is 5°. There are discrepancies between the current measured values ​​and the target values: the dominant direction parameter difference is 7°, and the dispersion parameter difference is 2.5°. By solving the system of equations, these parameter differences are converted into the required adjustments for the sound wave propagation time difference and the peak amplitude difference. Calculations show that the sound wave propagation time difference needs to be increased by 0.42 μs, and the peak amplitude difference needs to be decreased by 0.13 MPa. There is a corresponding relationship between the sound wave propagation parameters and the output power of the sound field excitation source: the propagation time difference is mainly controlled by the frequency difference, and the peak amplitude difference is mainly controlled by the power ratio. Calculations using the sound wave propagation model show that the output power of the 3mm depth sound field excitation source needs to be increased from the original 220W to 265W, and the output power of the 6mm depth sound field excitation source needs to be decreased from the original 180W to 160W. The calculated power adjustment is applied to the corresponding sound field excitation source using an ultrasonic control system, and the changes in the sound wave propagation parameters are monitored in real time until the target values ​​are reached, thus completing the precise adjustment of the output power of adjacent sound field excitation sources.

[0133] In this embodiment, a method for adjusting the output power of the acoustic field excitation source is adopted by establishing a dynamic response function, providing a scientific basis and implementation path for the modulation control of copper alloy rolled materials with multiple acoustic fields. Through detailed analysis of sound wave propagation characteristics and scientific construction of the dynamic response function, the influence of acoustic field parameter changes on dislocation orientation can be accurately predicted, and the output power configuration of the acoustic field excitation source can be optimized accordingly. The optimized acoustic field effect can effectively adjust the slip direction and arrangement of dislocation lines, making the dislocation orientation more reasonable and the internal stress distribution more uniform.

[0134] A second aspect of the present invention provides a multi-sound-field coordinated modulation control system for copper alloy rolled materials, the system comprising:

[0135] The first unit is used to obtain information on the microstructure and target rolling performance of the copper alloy material to be rolled.

[0136] The second unit is used to divide the calendered material into multiple acoustic field action zones along the thickness direction according to the microstructure information, and to set up acoustic field excitation sources in each zone.

[0137] The third unit is used to detect the propagation time of each sound field excitation source at different depth layers, calculate the propagation time difference, and convert the propagation time difference into phase modulation parameters.

[0138] The fourth unit is used to extract the microstructure density distribution curve of the calendered material, establish the inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic field excitation source, and calculate the energy distribution coefficient of each depth layer.

[0139] The fifth unit is used to drive each sound field excitation source to apply a synergistic sound field effect to the calendered material according to the phase control parameters and energy distribution coefficient, collect the dislocation line orientation distribution of the calendered material under the synergistic sound field effect, statistically analyze the angle distribution between the dislocation lines and the calendering direction, and calculate the angle concentration index.

[0140] The sixth unit is used to determine the target enhancement coefficient of dislocation orientation based on the target rolling performance index, calculate the acoustic wave superposition gain according to the included angle concentration index, compare the acoustic wave superposition gain with the target enhancement coefficient, and optimize the phase difference and energy ratio between adjacent acoustic field excitation sources based on the comparison results, so as to realize the coordinated control of dislocation orientation of rolling material.

[0141] A third aspect of the present invention provides an electronic device, comprising:

[0142] processor;

[0143] Memory used to store processor-executable instructions;

[0144] The processor is configured to invoke instructions stored in the memory to execute the aforementioned method.

[0145] A fourth aspect of the present invention provides a computer-readable storage medium having stored thereon computer program instructions that, when executed by a processor, implement the aforementioned method.

[0146] This invention can be a method, apparatus, system, and / or computer program product. The computer program product may include a computer-readable storage medium having computer-readable program instructions loaded thereon for performing various aspects of the invention.

[0147] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for modulating and controlling copper alloy rolled material using multi-field acoustic coordination, characterized in that, include: Obtain information on the microstructure and target rolling performance of the copper alloy material to be rolled; Based on the microstructure information, the calendered material is divided into multiple acoustic field zones along the thickness direction, and acoustic field excitation sources are arranged in each zone. The propagation time of each sound field excitation source at different depth layers is detected, the propagation time difference is calculated, and the propagation time difference is converted into a phase modulation parameter. Extract the microstructure density distribution curve of the calendered material, establish the inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic field excitation source, and calculate the energy distribution coefficient of each depth layer; Based on the phase modulation parameters and energy distribution coefficient, each sound field excitation source is driven to apply a synergistic sound field effect to the calendered material. The orientation distribution of dislocation lines in the calendered material under the synergistic sound field effect is collected, the angle distribution between the dislocation lines and the calendering direction is statistically analyzed, and the angle concentration index is calculated. The target enhancement coefficient for dislocation orientation is determined based on the target rolling performance index. The acoustic wave superposition gain is calculated based on the included angle concentration index. The acoustic wave superposition gain is compared with the target enhancement coefficient. Based on the comparison result, the phase difference and energy ratio between adjacent acoustic field excitation sources are optimized to achieve coordinated control of dislocation orientation in the rolled material.

2. The method according to claim 1, characterized in that, Based on the microstructure information, the calendered material is divided into multiple acoustic field zones along the thickness direction, and acoustic field excitation sources are arranged in each zone, including: Based on the microstructure information, extract the microstructure parameters of the calender at different thickness positions, calculate the microstructure difference value between adjacent thickness positions, determine the microstructure abrupt change position along the thickness direction of the calender based on the microstructure difference value, divide the calender into multiple acoustic field action zones with the microstructure abrupt change position as the boundary, and calculate the microstructure characteristic value of each acoustic field action zone. The number of sound field excitation sources required in each sound field action zone is calculated based on the tissue characteristic values ​​and the size of the sound field action zone. A sound field superposition experiment is conducted on the sound field excitation sources to obtain the action radius of the sound field excitation sources. The layout spacing of the sound field excitation sources is determined according to the action radius and the number of sound field excitation sources. Calculate the layout coordinates of the sound field excitation source and the corresponding sound field intensity distribution based on the layout spacing and the size of the sound field effect zone. Select the layout coordinates that satisfy the preset sound field coverage threshold, and place the sound field excitation source at the selected layout coordinates.

3. The method according to claim 1, characterized in that, The propagation time of each sound field excitation source at different depth layers is detected, the propagation time difference is calculated, and the propagation time difference is converted into phase modulation parameters, including: Acquire sound wave propagation signals from various sound field excitation sources at different depth layers, and perform time-frequency domain joint analysis on the sound wave propagation signals to extract the sound wave frequency drift and sound wave amplitude attenuation rate; The dislocation density distribution of each depth layer is inverted by the acoustic frequency drift, and the grain boundary density distribution of each depth layer is inverted by the acoustic amplitude attenuation rate. The acoustic impedance correction value of each depth layer is calculated based on the dislocation density distribution and the grain boundary density distribution. The propagation time of each depth layer is corrected according to the acoustic impedance correction value to obtain the corrected propagation time. The propagation time difference is obtained by calculating the difference between the corrected propagation times. The propagation time difference is then converted into the phase of the operating frequency of each sound field excitation source to obtain the excitation phase difference value. A spatiotemporal distribution map of the excitation phase difference value as the depth layer position changes is then constructed. Extract the spatial gradient vector field of the excitation phase difference from the spatiotemporal distribution map, identify the depth layer position with the largest gradient amplitude in the spatial gradient vector field, and take the excitation phase difference of the sound field excitation source corresponding to the depth layer position with the largest gradient amplitude as the reference phase difference. The phase deviation between the excitation phase difference of other depth layer sound field excitation sources and the reference phase difference is calculated as the phase control parameter.

4. The method according to claim 1, characterized in that, Extract the microstructure density distribution curve of the calendered material, establish the inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic field excitation source, and calculate the energy distribution coefficient of each depth layer, including: Obtain the microstructure density values ​​of calendered material at different depths along the thickness direction, arrange the microstructure density values ​​according to the depth layer position, and construct a microstructure density distribution curve. Calculate the rate of change of tissue density values ​​between adjacent depth layers, identify depth layers whose rate of change deviates from the mean, and mark them as areas of abnormal tissue density. Read the tissue density values ​​of the abnormal tissue density area, establish the acoustic wave propagation impedance distribution, deduce the acoustic wave propagation path based on the acoustic wave propagation impedance distribution, and calculate the cumulative attenuation of acoustic wave energy. The cumulative attenuation of acoustic energy is used as the incremental energy requirement of the acoustic field excitation source at each depth layer, and an inverse mapping relationship between the tissue density distribution curve and the output energy of the acoustic field excitation source is established. Read the incremental output energy demand of the sound field excitation source at each depth layer, and construct the depth transmission link of the incremental output energy demand of the sound field excitation source; Identify the anchor depth layer where abrupt changes occur in the depth transmission link, and perform bidirectional diffusion correction on the incremental output energy demand of the sound field excitation source in adjacent depth layers starting from the anchor depth layer. Normalize the corrected incremental output energy demand of the sound field excitation source to obtain the energy distribution coefficient of each depth layer.

5. The method according to claim 1, characterized in that, Based on the phase modulation parameters and energy distribution coefficients, each acoustic excitation source is driven to apply a synergistic acoustic field to the calendered material. The orientation distribution of dislocation lines in the calendered material under the synergistic acoustic field is collected, and the angle distribution between the dislocation lines and the calendering direction is statistically analyzed. The angle concentration index is calculated, including: The phase modulation parameters are converted into frequency control quantities of the sound field excitation source, and the energy distribution coefficient is converted into power control quantities of the sound field excitation source. Sound waves are generated based on the frequency control quantities and power control quantities. The propagation time and pressure amplitude of sound waves in calendered material are collected, and the propagation time difference and pressure amplitude difference between adjacent sound field excitation sources are calculated to obtain the propagation characteristics of sound waves in the depth direction. The beam pointing angle of the sound field excitation source is adjusted according to the propagation characteristics. The beam pointing angle is combined with the pressure amplitude difference to determine the sound wave superposition position and form a cooperative sound field. The stress state of the calendered material is collected in the area of ​​the cooperative sound field. The driving force of dislocation movement is calculated according to the stress state to determine the starting point of dislocation slip. Record the position coordinates of the dislocation line from the starting point, construct the spatial trajectory of the dislocation line based on the position coordinates, extract the node information of the dislocation line, calculate the local direction of the dislocation line based on the node information, and project the local direction of the dislocation line onto the rolling direction to obtain the dislocation orientation distribution. The included angle values ​​in the dislocation orientation distribution are statistically analyzed, the number distribution of included angle values ​​is calculated, and the number distribution is normalized to obtain the included angle concentration index.

6. The method according to claim 1, characterized in that, The target enhancement coefficient for dislocation orientation is determined based on the target rolling performance index. The acoustic wave superposition gain is calculated based on the included angle concentration index. The acoustic wave superposition gain is compared with the target enhancement coefficient. Based on the comparison results, the phase difference and energy ratio between adjacent acoustic field excitation sources are optimized to achieve coordinated control of dislocation orientation in the rolled material. This includes: Obtain the stress and texture requirements in the target calendering performance indicators, calculate the crystal orientation distribution and dislocation slip distribution of the calendered material, determine the critical shear force for dislocation motion based on the crystal orientation distribution and dislocation slip distribution, and convert the critical shear force into the target enhancement coefficient for dislocation orientation; Acquire sound wave propagation signals from adjacent sound field excitation sources, calculate the sound wave propagation time difference and pressure amplitude difference, extract the angle value distribution from the angle concentration index, establish a dynamic response function between the angle value distribution and the sound wave propagation signal, adjust the output power of adjacent sound field excitation sources according to the dynamic response function, and convert the adjusted sound wave propagation signal into sound wave superposition gain. The difference between the acoustic wave superposition gain and the target enhancement coefficient is used as the control deviation of dislocation orientation. The phase compensation value and power distribution ratio of adjacent acoustic field excitation sources are calculated based on the control deviation. The output parameters of the acoustic field excitation sources are adjusted based on the phase compensation value and power distribution ratio to achieve coordinated control of dislocation orientation in calendered material.

7. The method according to claim 6, characterized in that, Establishing a dynamic response function between the included angle distribution and the sound wave propagation signal, and adjusting the output power of adjacent sound field excitation sources according to the dynamic response function includes: Statistical analysis of the angle distribution data yields the peak position and distribution width of the angle values. The peak position is used as the dominant orientation parameter of dislocation orientation, and the distribution width is used as the dispersion parameter of dislocation orientation. Acquire sound wave propagation signals from adjacent sound field excitation sources inside the calendered material, perform time-domain analysis on the sound wave propagation signals to extract the arrival time and peak amplitude of the sound waves, and calculate the sound wave propagation time difference and peak amplitude difference between adjacent sound field excitation sources. The dominant direction parameter and the dispersion parameter are set as the output variables of the dynamic response function, and the sound wave propagation time difference and the sound wave peak amplitude difference are set as the input variables of the dynamic response function. The dynamic response function of the angle value distribution and the sound wave propagation signal is constructed by fitting the functional relationship between the input variables and the output variables. Calculate the partial derivatives of the dynamic response function with respect to the sound wave propagation time difference and the sound wave peak amplitude difference. Based on the partial derivatives, convert the target values ​​of the sound wave propagation time difference and the sound wave peak amplitude difference into the output power adjustment amount of the adjacent sound field excitation source. Apply the output power adjustment amount to the corresponding sound field excitation source to complete the adjustment of the output power of the adjacent sound field excitation source.

8. A multi-sound-field coordinated modulation control system for copper alloy rolled materials, used to implement the method of any one of claims 1-7, characterized in that, include: The first unit is used to obtain information on the microstructure and target rolling performance of the copper alloy material to be rolled. The second unit is used to divide the calendered material into multiple acoustic field action zones along the thickness direction according to the microstructure information, and to set up acoustic field excitation sources in each zone. The third unit is used to detect the propagation time of each sound field excitation source at different depth layers, calculate the propagation time difference, and convert the propagation time difference into phase modulation parameters. The fourth unit is used to extract the microstructure density distribution curve of the calendered material, establish the inverse mapping relationship between the microstructure density distribution curve and the output energy of the acoustic field excitation source, and calculate the energy distribution coefficient of each depth layer. The fifth unit is used to drive each sound field excitation source to apply a synergistic sound field effect to the calendered material according to the phase control parameters and energy distribution coefficient, collect the dislocation line orientation distribution of the calendered material under the synergistic sound field effect, statistically analyze the angle distribution between the dislocation lines and the calendering direction, and calculate the angle concentration index. The sixth unit is used to determine the target enhancement coefficient of dislocation orientation based on the target rolling performance index, calculate the acoustic wave superposition gain according to the included angle concentration index, compare the acoustic wave superposition gain with the target enhancement coefficient, and optimize the phase difference and energy ratio between adjacent acoustic field excitation sources based on the comparison results, so as to realize the coordinated control of dislocation orientation of rolling material.

9. An electronic device, characterized in that, include: processor; Memory used to store processor-executable instructions; The processor is configured to invoke instructions stored in the memory to execute the method according to any one of claims 1 to 7.

10. A computer-readable storage medium having computer program instructions stored thereon, characterized in that, When the computer program instructions are executed by the processor, they implement the method described in any one of claims 1 to 7.

Images