merhaba iyi gunler su sekilde bir kodum var 'i'nin her degeri için matirse bir değer kaydetmek ve en son o degerlerin toptan graifigini cizdirmek istiyorum yardımcı olur musunuz lutfen?
Sub staticforce()
Weight = 35.06 + 1 ' in kg ; weight of one landing gear including wheel, tire, oil etc. (+1 kg for oil assumed)
Z_axis_of_CM = -42.2885 'in inch; center of gravity of one landing gear
X_axis_of_CM = 110.2879 'in inch; center of gravity of one landing gear
Z_axis_of_A = -15.362 'in inch; trunion pin coordinate
X_axis_of_A = 119.971 'in inch; trunion pin coordinate
X_axis_of_F = 117.3885 'in inch; actuator & LG attachment coordinate
Z_axis_of_F = -24.1542 'in inch; actuator & LG attachment coordinate
X_axis_of_D = 149.777 'in inch; actuator & Aircraft attachment coordinate
Z_axis_of_D = -12.27 'in inch; actuator & Aircraft attachment coordinate
Dim static_force_matrix() As Double
Dim angle_matrix() As Double
Dim actuator_lenght_matrix() As Double
x = ((X_axis_of_A - X_axis_of_F) ^ 2 + (Z_axis_of_A - Z_axis_of_F) ^ 2) ^ 0.5
y = ((X_axis_of_A - X_axis_of_D) ^ 2 + (Z_axis_of_A - Z_axis_of_D) ^ 2) ^ 0.5
distance_between_AB_line_and_CM = ((X_axis_of_A - X_axis_of_CM) ^ 2 + (Z_axis_of_A - Z_axis_of_CM) ^ 2) ^ 0.5
angle1 = Application.Acos((distance_between_AB_line_and_CM ^ 2 + (X_axis_of_A - X_axis_of_F) ^ 2 + (Z_axis_of_F - Z_axis_of_A) ^ 2 - (X_axis_of_F - X_axis_of_CM) ^ 2 - (Z_axis_of_CM - Z_axis_of_F) ^ 2) / 2 / distance_between_AB_line_and_CM / (((X_axis_of_A - X_axis_of_F) ^ 2 + (Z_axis_of_F - Z_axis_of_A) ^ 2)) ^ 0.5)
angle2 = Application.Asin((Z_axis_of_D - Z_axis_of_A) / y)
fully_down_angle = Atn((Abs(Z_axis_of_CM - Z_axis_of_A)) / (Abs(X_axis_of_CM - X_axis_of_A)))
fully_up_angle = Pi * 1.1
Angle = fully_down_angle
i = 0
While Angle < fully_up_angle
actuator_lenght_DF = ((x ^ 2 + y ^ 2) - (2 * x * y) * Cos(Pi - angle1 + angle2 - Angle)) ^ 0.5
static_force = Weight * 9.81 * distance_between_AB_line_and_CM * Sin(Angle - Pi / 2) / ((x * Cos(Pi - Angle - angle1) * (y * Sin(angle2) + x * Sin(Pi - Angle - angle1)) / actuator_lenght_DF + x * Sin(Pi - Angle - angle1) * (X_axis_of_D - X_axis_of_A - x * Cos(Pi - Angle - angle1)) / (actuator_lenght_DF)))
cg = distance_between_AB_line_and_CM * Sin(Angle - Pi / 2)
aa = x * Cos(Pi - Angle - angle1)
epsilon_angle1 = Application.Asin((y * Sin(angle2) + x * Sin(Pi - Angle - angle1)) / actuator_lenght_DF) * 180 / Pi
epsilon_angle2 = Application.Acos((X_axis_of_D - X_axis_of_A - x * Cos(Pi - Angle - angle1)) / (actuator_lenght_DF)) * 180 / Pi
bb = x * Sin(Pi - Angle - angle1)
If actuator_lenght_DF < 20.8592 Then 'in inch
GoTo 1020
End If
Angle = Angle + 0.0001
i = i + 1
static_force_matrix(i) = static_force 'Newton
angle_matrix(i) = Angle * 180 / Pi 'degree
actuator_lenght_matrix(i) = actuator_lenght_DF
Wend
1020
End Sub
Sub staticforce()
Weight = 35.06 + 1 ' in kg ; weight of one landing gear including wheel, tire, oil etc. (+1 kg for oil assumed)
Z_axis_of_CM = -42.2885 'in inch; center of gravity of one landing gear
X_axis_of_CM = 110.2879 'in inch; center of gravity of one landing gear
Z_axis_of_A = -15.362 'in inch; trunion pin coordinate
X_axis_of_A = 119.971 'in inch; trunion pin coordinate
X_axis_of_F = 117.3885 'in inch; actuator & LG attachment coordinate
Z_axis_of_F = -24.1542 'in inch; actuator & LG attachment coordinate
X_axis_of_D = 149.777 'in inch; actuator & Aircraft attachment coordinate
Z_axis_of_D = -12.27 'in inch; actuator & Aircraft attachment coordinate
Dim static_force_matrix() As Double
Dim angle_matrix() As Double
Dim actuator_lenght_matrix() As Double
x = ((X_axis_of_A - X_axis_of_F) ^ 2 + (Z_axis_of_A - Z_axis_of_F) ^ 2) ^ 0.5
y = ((X_axis_of_A - X_axis_of_D) ^ 2 + (Z_axis_of_A - Z_axis_of_D) ^ 2) ^ 0.5
distance_between_AB_line_and_CM = ((X_axis_of_A - X_axis_of_CM) ^ 2 + (Z_axis_of_A - Z_axis_of_CM) ^ 2) ^ 0.5
angle1 = Application.Acos((distance_between_AB_line_and_CM ^ 2 + (X_axis_of_A - X_axis_of_F) ^ 2 + (Z_axis_of_F - Z_axis_of_A) ^ 2 - (X_axis_of_F - X_axis_of_CM) ^ 2 - (Z_axis_of_CM - Z_axis_of_F) ^ 2) / 2 / distance_between_AB_line_and_CM / (((X_axis_of_A - X_axis_of_F) ^ 2 + (Z_axis_of_F - Z_axis_of_A) ^ 2)) ^ 0.5)
angle2 = Application.Asin((Z_axis_of_D - Z_axis_of_A) / y)
fully_down_angle = Atn((Abs(Z_axis_of_CM - Z_axis_of_A)) / (Abs(X_axis_of_CM - X_axis_of_A)))
fully_up_angle = Pi * 1.1
Angle = fully_down_angle
i = 0
While Angle < fully_up_angle
actuator_lenght_DF = ((x ^ 2 + y ^ 2) - (2 * x * y) * Cos(Pi - angle1 + angle2 - Angle)) ^ 0.5
static_force = Weight * 9.81 * distance_between_AB_line_and_CM * Sin(Angle - Pi / 2) / ((x * Cos(Pi - Angle - angle1) * (y * Sin(angle2) + x * Sin(Pi - Angle - angle1)) / actuator_lenght_DF + x * Sin(Pi - Angle - angle1) * (X_axis_of_D - X_axis_of_A - x * Cos(Pi - Angle - angle1)) / (actuator_lenght_DF)))
cg = distance_between_AB_line_and_CM * Sin(Angle - Pi / 2)
aa = x * Cos(Pi - Angle - angle1)
epsilon_angle1 = Application.Asin((y * Sin(angle2) + x * Sin(Pi - Angle - angle1)) / actuator_lenght_DF) * 180 / Pi
epsilon_angle2 = Application.Acos((X_axis_of_D - X_axis_of_A - x * Cos(Pi - Angle - angle1)) / (actuator_lenght_DF)) * 180 / Pi
bb = x * Sin(Pi - Angle - angle1)
If actuator_lenght_DF < 20.8592 Then 'in inch
GoTo 1020
End If
Angle = Angle + 0.0001
i = i + 1
static_force_matrix(i) = static_force 'Newton
angle_matrix(i) = Angle * 180 / Pi 'degree
actuator_lenght_matrix(i) = actuator_lenght_DF
Wend
1020
End Sub