2.1. Flexible PV Mounting Structure Geometric Model
The constructed flexible PV support model consists of six spans, each with a span of 2 m. The spans are connected by struts, with the support cables having a height of 4.75 m, directly supporting the PV panels. The wind-resistant cables are 4 m high and are connected to the lower ends of the struts. The end support beams are 4 m high, with tie rods connected to the end support beams at a 45° angle, each measuring 5.657 m in length. There are six sets of struts, spaced 2 m apart.
In accordance with the “Steel Strand for Prestressed Concrete” (GB/T 5224-2014) [
34], we selected 1 × 7 unbonded prestressed hot-dip galvanized steel strands as the load-bearing cables, which have the following properties: tensile strength
Rm = 1860 MPa, elastic modulus
E0 = 1.95 × 10
5 MPa, linear expansion coefficient
αT = 1.17 × 10
−5/°C, and Poisson’s ratio
ν = 0.3. The diameter of the support cables is 0.0127 m, while the wind-resistant cables have a diameter of 0.0152 m. The end support beams are made of HPB300 steel, with cross-sectional dimensions of 0.2 m in length and width, and a wall thickness of 0.01 m. The columns are constructed from Q355 seamless steel pipes, having an outer diameter of 0.2 m and an inner diameter of 0.05 m, with an elastic modulus of 2.06 × 10
5 N/mm
2. The struts are solid HPB300 steel rods with a diameter of 0.05 m. To simulate the PV panels, a virtual surface was employed, applying a uniform distributed load of 0.15 kN/m
2 to represent the self-weight of the PV modules. The geometric model of the flexible PV support system and the nomenclature of its components are illustrated in
Figure 1.
2.2. Boundary Condition Configuration
For all components connected to the ground, the nodes are constrained in all six degrees of freedom (DOFs): translational in the x, y, and z directions, and rotational about the x, y, and z axes. The nodes along the upper edges on both sides of the flexible PV support structure are also fixed in all six DOFs. For the nodes on the middle support beams of the flexible PV structure, constraints are applied in the x and z directions. The sections of the columns interfacing with the ground are modeled using fixed supports to simulate full constraint conditions. Given the significant geometric nonlinearity inherent in the flexible PV support system, the analysis incorporates nonlinear approaches, specifically selecting the P-Δ effect and large displacement effects. The time step is set to 1000, with a time interval of 0.1 s. Given that the direct integration method is used for time–history analysis, Rayleigh damping is applied. Calculations yield a modal damping ratio of 0.02 for the first two modes, with the proportional damping coefficients α = 0.304 and β = 0.00132.
2.3. Simulation of Fluctuating Wind Load Time History
The fluctuating wind speed time series
v(
t) at time
t is represented as a linear combination of the time series from several previous moments, plus an independent stochastic process at time
t. The expression is as follows:
where
v(
t) is the fluctuating wind speed at time
t;
is the autoregressive coefficient;
p is the order of the AR model;
v(
t −
kΔ
t) is the fluctuating wind speed at
k previous time steps; Δ
t is the time step of the fluctuating wind speed; and
N(
t) is an independent stochastic process with a mean of 0 and variance
. Multiplying both sides of Equation (2) by
v(
t −
kΔ
t) gives:
In Equation (3),
j = 1, 2,…,
p. Taking the mathematical expectation on both sides of Equation (2)
Considering the properties of the autocorrelation function: the autocorrelation function of a stochastic process
X(
t) is defined as the mean of
. If the process is a stationary stochastic process, its mean
is independent of the absolute value of time, and only depends on the time difference, that is:
Since the mean of
N(
t) is 0 and is independent of
v(
t), the relationship between the correlation function
and the autoregressive coefficient
is:
The autoregressive coefficients
in the equation can be determined by the following system of equations:
where
is the correlation function of the fluctuating wind speed at
t = jΔ
t, which can be obtained by Fourier transform from the target power spectral density function, and can be calculated as follows:
In Equation (6), is the frequency of the fluctuating wind speed. Substituting the obtained autoregressive coefficients and the independent stochastic process N(t) at time t into Equation (1), the single-point fluctuating wind speed time series can be obtained.
In MATLAB (2021b), programming was conducted to simulate the fluctuating wind speed time series for a single node based on the Davenport wind spectrum, utilizing the auto-regressive (AR) technique. The simulation does not account for the correlation of fluctuating wind speeds between different nodes. The duration of the fluctuating wind speed time series was set to 100 s with a time step of 0.1 s, and the AR model was configured with a fourth-order process. The standard wind speed was set to 30 m/s, and the frequency integration range was chosen from 0.001 Hz to 10 Hz [
35]; the simulation results of the fluctuating wind flow are shown in
Figure 2. As depicted in
Figure 3, the simulated fluctuating wind speeds predominantly vary within the range of −10 m/s to 10 m/s, with a mean value oscillating around zero, indicating that the simulation results for the fluctuating wind are realistic and satisfactory. Furthermore, as shown in
Figure 4, the simulated spectrum closely aligns with the target spectrum, demonstrating the efficacy of the AR method in replicating the desired wind speed characteristics.
The fluctuating wind speed time series for the flexible PV support structure was simulated using the linear filtering method. This time series was then converted into a fluctuating wind load time series. The relationship between wind speed and wind pressure is given by the following equation:
In Equation (8),
ω denotes the wind pressure corresponding to the wind speed; ρ represents the air density, assumed to be 1.225 kg/m
3 in this context; and
stands for the instantaneous wind speed. Utilizing the given equation and the relationship between the mean wind speed
and the fluctuating wind speed
, the time history of the fluctuating wind load can be derived.
In Equation (9),
is the shape coefficient, taken as 1.4 [
36], and
A represents the area of the PV panel.
Figure 5 illustrates the time history of the fluctuating wind load. After this conversion, fluctuating wind can be added to SAP2000. The obtained data can be directly utilized for structural dynamic analysis in SAP2000.