Temperaturfeld in Chart darstellen

affrop · 18. Dezember 2013, 13:07

hallo leute,

ich habe in kleines proggi geschrieben, wass temperaturfelder speziell bei erdverlegten rohrleitungen im boden berechnet und dass dann graphisch in einem chart darstellen soll.
hier mal der code dazu:

Spoiler anzeigen

VB.NET-Quellcode

Private Sub Button1_Click(sender As Object, e As EventArgs) Handles Button1.Click
d = Double.Parse(TextBox1.Text)
TM = Double.Parse(TextBox2.Text)
s = Double.Parse(TextBox3.Text)
Lambda = Double.Parse(TextBox4.Text)
z = Double.Parse(TextBox5.Text)
Chart1.Series.Add("Rohr")
Chart1.Series.Add("Serie1")
Chart1.Series("Rohr").ChartType = DataVisualization.Charting.SeriesChartType.Line
Chart1.Series("Serie1").ChartType = DataVisualization.Charting.SeriesChartType.Point
LR = Math.Log((d + 2 * s) / d) / (2 * Lambda * Math.PI)
j = 0
For Me.alpha = 0 To Math.Atan(12 / z)
i = 0
z1 = z / Math.Cos(alpha)
LE = arcosh((2 * z1) / ((d + 2 * s) / 1000)) / (4 * Math.PI)
q = (TM - 10) / (LR + LE)
T0 = 10 + q * LE
Do
Tx = T0 - i * ((T0 - 10) / (z1 - (d / 2000) - (s / 1000)))
x = z * Math.Tan(alpha)
y = ((d / 2000) + (s / 1000) + i) * Math.Cos(alpha)
Chart1.Series("Serie1").Points.AddXY(x, y)
If Tx >= 45 Then
Chart1.Series("Serie1").Points(j).Color = Color.Red
Else
If Tx >= 35 And Tx < 45 Then
Chart1.Series("Serie1").Points(j).Color = Color.Violet
Else
Chart1.Series("Serie1").Points(j).Color = Color.Blue
End If
End If
i = i + 0.001
j = j + 1
Loop Until i > (z1 - (d / 2000) - (s / 1000))
Next
For x As Double = -((d + 2 * s) / 1000) To ((d + 2 * s) / 1000) Step 0.001
Chart1.Series("Rohr").Points.AddXY(x, Math.Sqrt((((d + 2 * s) / 1000) ^ 2) - (x ^ 2)))
Chart1.Series("Rohr").Points.AddXY(x, -Math.Sqrt((((d + 2 * s) / 1000) ^ 2) - (x ^ 2)))
Next
Chart1.ChartAreas("ChartArea1").AxisX.Minimum = 0
Chart1.ChartAreas("ChartArea1").AxisX.Maximum = 12
Chart1.ChartAreas("ChartArea1").AxisY.Minimum = 0
Chart1.ChartAreas("ChartArea1").AxisY.Maximum = 2
End Sub

die berechnung läuft auch super, aber die bei der punktweisen darstellung im chart gibt es probleme. die punkte werden nicht vollflächig verteilt sonder liegen nur als eine art balken in dem chart.

hier mal ein bild, wie es aussieht:

hat jemand ne idee was da schief läuft?

ErfinderDesRades · 18. Dezember 2013, 13:16

dazu müsste man wissen, wies aussehen soll.

affrop · 18. Dezember 2013, 13:22

na halt ein temperaturfeld und kein balken.
die farbverläufe sollen in dem gezeigten bereich vollflächig dargestellt werden und nicht nur als balken wie es im moment ist.

RodFromGermany · 18. Dezember 2013, 17:25

Gib der Series mal den ChartType Line.

affrop · 20. Dezember 2013, 08:20

dann hab ich das problem mit den farbverläufen weil ich ja dann für jede farbe ne eigene series machen müsste. im moment kann ich jeden einzelnen punkt definieren und einfärbern.

RodFromGermany · 20. Dezember 2013, 13:35

affrop schrieb:

weil ich ja dann für jede farbe ne eigene series machen müsste.

So isses.

affrop · 20. Dezember 2013, 18:29

so ich hab jetzt für jeden berechnungssektor ne eigene series erstellt. problem ist dies hier:

Spoiler anzeigen

das kommt so nach etwa 3,2 mio. datenpunkten, die in das chart geschrieben wurden. das ist nicht mal der komplette 1. sektor.
wenn ich für jede farbe ne neue series mit datenpunkten füllen will geht es bis 3,8 mio. dann kotzt das framework ab und macht en schönes rotes kreuz durch das chart.
hier mal en screenshot mit leistungsdaten:

hat jemand ne idee wie man das problem lösen kann?
hier mal der gesamte code mit den sektoren pro series dazu:

Spoiler anzeigen

VB.NET-Quellcode

Public Class Form1
Dim d As Double
Dim TM As Double
Dim s As Double
Dim Lambda As Double
Dim z As Double
Dim z1 As Double
Dim q As Double
Dim T0 As Double
Dim Tx As Double
Dim x As Double
Dim y As Double
Dim i As Double
Dim j As Int64
Dim h As Double
Dim LR As Double
Dim LE As Double
Dim alpha As Double
Dim Tz As Double
Dim TxListSektor1 As New ArrayList
Dim TxListSektor2 As New ArrayList
Dim TxListSektor3 As New ArrayList
Dim TxListSektor4 As New ArrayList
Dim xListSektor1 As New ArrayList
Dim xListSektor2 As New ArrayList
Dim xListSektor3 As New ArrayList
Dim xListSektor4 As New ArrayList
Dim yListSektor1 As New ArrayList
Dim yListSektor2 As New ArrayList
Dim yListSektor3 As New ArrayList
Dim yListSektor4 As New ArrayList
Dim PBValue1 As Double
Dim PBValue2 As Double
Dim PBValue3 As Double
Dim PBValue4 As Double
Dim PBValue5 As Double
Delegate Sub ProgressBar(ByVal PBValue1 As Double, ByVal PBValue2 As Double, ByVal PBValue3 As Double, ByVal PBValue4 As Double, ByVal PBValue5 As Double)
Delegate Sub Charting1(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Delegate Sub Charting2(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Delegate Sub Charting3(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Delegate Sub Charting4(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Delegate Sub Points(ByVal k As Int64)
Dim Min As Integer = 0
Dim Sek As Integer = 0
Public Shared Function arcosh(ByVal x As Double) As Double
arcosh = Math.Log(x + Math.Sqrt((x ^ 2) - 1))
End Function
Private Sub Form_Load(sender As Object, e As EventArgs) Handles MyBase.Load
Label11.Text = "Berechnungszeit: " & Min & " Min. " & Sek & " Sek."
End Sub
Private Sub Button2_Click(sender As Object, e As EventArgs) Handles Button2.Click
Application.Exit()
End Sub
Private Sub Button1_Click(sender As Object, e As EventArgs) Handles Button1.Click
d = Double.Parse(TextBox1.Text)
TM = Double.Parse(TextBox2.Text)
s = Double.Parse(TextBox3.Text)
Lambda = Double.Parse(TextBox4.Text)
z = Double.Parse(TextBox5.Text)
Chart1.Series.Add("Rohr")
Chart1.Series.Add("Sektor1")
Chart1.Series.Add("Sektor2")
Chart1.Series.Add("Sektor3")
Chart1.Series.Add("Sektor4")
Chart1.Series("Rohr").ChartType = DataVisualization.Charting.SeriesChartType.Line
Chart1.Series("Sektor1").ChartType = DataVisualization.Charting.SeriesChartType.Point
Chart1.Series("Sektor2").ChartType = DataVisualization.Charting.SeriesChartType.Point
Chart1.Series("Sektor3").ChartType = DataVisualization.Charting.SeriesChartType.Point
Chart1.Series("Sektor4").ChartType = DataVisualization.Charting.SeriesChartType.Point
Chart1.Visible = False
Timer1.Start()
BackgroundWorker1.RunWorkerAsync()
End Sub
Private Sub Backgroundworker1_DoWork(sender As Object, e As System.ComponentModel.DoWorkEventArgs) Handles BackgroundWorker1.DoWork
Dim PBGel As New ProgressBar(AddressOf Me.ProgressBarChange)
LR = Math.Log((d + 2 * s) / d) / (2 * Lambda * Math.PI)
j = 0
For Me.alpha = 0 To ((180 * Math.Atan(12 / z)) / Math.PI) Step 0.1
i = 0
z1 = z / Math.Cos(alpha * (Math.PI / 180))
LE = arcosh((2 * z1) / ((d + 2 * s) / 1000)) / (4 * Math.PI)
q = (TM - 104.1) / (LR + LE)
T0 = 14.1 + q * LE
Do
Tx = T0 - i * ((T0 - 14.1) / (z1 - (d / 2000) - (s / 1000)))
x = ((d / 2000) + (s / 1000) + i) * Math.Sin(alpha * (Math.PI / 180))
y = ((d / 2000) + (s / 1000) + i) * Math.Cos(alpha * (Math.PI / 180))
TxListSektor1.Add(Tx)
xListSektor1.Add(x)
yListSektor1.Add(y)
i = i + 0.001
j = j + 1
Loop Until i > (z1 - (d / 2000) - (s / 1000))
PBValue1 = (alpha / ((180 * Math.Atan(12 / z)) / Math.PI)) * 100
Me.Invoke(PBGel, PBValue1, 0, 0, 0, 0)
Next
PBValue1 = 100
Me.Invoke(PBGel, 100, 0, 0, 0, 0)
For Me.alpha = 0 To ((180 * Math.Atan(z / 12)) / Math.PI) Step 0.1
i = 0
z1 = 12 / Math.Cos(alpha * (Math.PI / 180))
LE = arcosh((2 * z1) / ((d + 2 * s) / 1000)) / (4 * Math.PI)
h = 12 * Math.Tan(alpha * (Math.PI / 180))
Tz = 9.78 + ((14.1 - 9.78) * Math.Exp((-h) / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))) * Math.Cos(h / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000))))) - (13.16 * Math.Exp((-h) / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))) * Math.Sin(h / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))))
q = (TM - Tz) / (LR + LE)
T0 = Tz + q * LE
Do
Tx = T0 - i * ((T0 - Tz) / (z1 - (d / 2000) - (s / 1000)))
x = ((d / 2000) + (s / 1000) + i) * Math.Cos(alpha * (Math.PI / 180))
y = ((d / 2000) + (s / 1000) + i) * Math.Sin(alpha * (Math.PI / 180))
TxListSektor2.Add(Tx)
xListSektor2.Add(x)
yListSektor2.Add(y)
i = i + 0.001
j = j + 1
Loop Until i > (z1 - (d / 2000) - (s / 1000))
PBValue2 = (alpha / ((180 * Math.Atan(z / 12)) / Math.PI)) * 100
Me.Invoke(PBGel, 100, PBValue2, 0, 0, 0)
Next
PBValue2 = 100
Me.Invoke(PBGel, 100, 100, 0, 0, 0)
For Me.alpha = 0 To ((180 * Math.Atan((12 - z) / 12)) / Math.PI) Step 0.1
i = 0
z1 = 12 / Math.Cos(alpha * (Math.PI / 180))
LE = arcosh((2 * z1) / ((d + 2 * s) / 1000)) / (4 * Math.PI)
h = 12 * Math.Tan(alpha * (Math.PI / 180)) + z
Tz = 9.78 + ((14.1 - 9.78) * Math.Exp((0 - h) / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))) * Math.Cos(h / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000))))) - (13.16 * Math.Exp((0 - h) / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))) * Math.Sin(h / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))))
q = (TM - Tz) / (LR + LE)
T0 = Tz + q * LE
Do
Tx = T0 - i * ((T0 - Tz) / (z1 - (d / 2000) - (s / 1000)))
x = ((d / 2000) + (s / 1000) + i) * Math.Cos(alpha * (Math.PI / 180))
y = 0 - (((d / 2000) + (s / 1000) + i) * Math.Sin(alpha * (Math.PI / 180)))
TxListSektor3.Add(Tx)
xListSektor3.Add(x)
yListSektor3.Add(y)
i = i + 0.001
j = j + 1
Loop Until i > (z1 - (d / 2000) - (s / 1000))
PBValue3 = (alpha / ((180 * Math.Atan((12 - z) / 12)) / Math.PI)) * 100
Me.Invoke(PBGel, 100, 100, PBValue3, 0, 0)
Next
PBValue3 = 100
Me.Invoke(PBGel, 100, 100, 100, 0, 0)
For Me.alpha = 0 To (180 * (Math.Atan(12 / (12 - z)) / Math.PI)) Step 0.1
i = 0
z1 = (12 - z) / Math.Cos(alpha * (Math.PI / 180))
LE = arcosh((2 * z1) / ((d + 2 * s) / 1000)) / (4 * Math.PI)
q = (TM - 9.78) / (LR + LE)
T0 = 9.78 + q * LE
Do
Tx = T0 - i * ((T0 - 9.78) / (z1 - (d / 2000) - (s / 1000)))
x = ((d / 2000) + (s / 1000) + i) * Math.Sin(alpha * (Math.PI / 180))
y = 0 - (((d / 2000) + (s / 1000) + i) * Math.Cos(alpha * (Math.PI / 180)))
TxListSektor4.Add(Tx)
xListSektor4.Add(x)
yListSektor4.Add(y)
i = i + 0.001
j = j + 1
Loop Until i > (z1 - (d / 2000) - (s / 1000))
PBValue4 = (alpha / (180 * (Math.Atan(12 / (12 - z)) / Math.PI))) * 100
Me.Invoke(PBGel, 100, 100, 100, PBValue4, 0)
Next
PBValue4 = 100
Me.Invoke(PBGel, 100, 100, 100, 100, 0)
End Sub
Private Sub Backgroundworker1_RunWorkerCompleted(sender As Object, e As System.ComponentModel.RunWorkerCompletedEventArgs) Handles BackgroundWorker1.RunWorkerCompleted
BackgroundWorker2.RunWorkerAsync()
End Sub
Private Sub Backgroundworker2_DoWork(sender As Object, e As System.ComponentModel.DoWorkEventArgs) Handles BackgroundWorker2.DoWork
Dim ChartDel1 As New Charting1(AddressOf ChartFill1)
Dim ChartDel2 As New Charting2(AddressOf ChartFill2)
Dim ChartDel3 As New Charting3(AddressOf ChartFill3)
Dim ChartDel4 As New Charting4(AddressOf ChartFill4)
Dim PBGel As New ProgressBar(AddressOf Me.ProgressBarChange)
Dim kDel As New Points(AddressOf Me.Point)
Dim kges As Int64 = TxListSektor1.Count + TxListSektor2.Count + TxListSektor3.Count + TxListSektor4.Count
Dim k As Int64 = 1
Dim a As Int64
a = TxListSektor1.Count - 1
j = 0
i = a
Do
Tx = Double.Parse(TxListSektor1.Item(i).ToString)
x = Double.Parse(xListSektor1.Item(i).ToString)
y = Double.Parse(yListSektor1.Item(i).ToString)
PBValue5 = (Double.Parse(k.ToString) / Double.Parse(kges.ToString)) * 100
Me.Invoke(ChartDel1, Tx, x, y, j)
Me.Invoke(PBGel, 100, 100, 100, 100, PBValue5)
Me.Invoke(kDel, k)
TxListSektor1.RemoveAt(i)
xListSektor1.RemoveAt(i)
yListSektor1.RemoveAt(i)
i = i - 1
k = k + 1
j = j + 1
Loop Until j > a
a = TxListSektor2.Count - 1
j = 0
i = a
Do
Tx = Double.Parse(TxListSektor2.Item(i).ToString)
x = Double.Parse(xListSektor2.Item(i).ToString)
y = Double.Parse(yListSektor2.Item(i).ToString)
PBValue5 = (Double.Parse(k.ToString) / Double.Parse(kges.ToString)) * 100
Me.Invoke(ChartDel2, Tx, x, y, j)
Me.Invoke(PBGel, 100, 100, 100, 100, PBValue5)
Me.Invoke(kDel, k)
TxListSektor2.RemoveAt(i)
xListSektor2.RemoveAt(i)
yListSektor2.RemoveAt(i)
i = i - 1
k = k + 1
j = j + 1
Loop Until j > a
a = TxListSektor3.Count - 1
j = 0
i = a
Do
Tx = Double.Parse(TxListSektor3.Item(i).ToString)
x = Double.Parse(xListSektor3.Item(i).ToString)
y = Double.Parse(yListSektor3.Item(i).ToString)
PBValue5 = (Double.Parse(k.ToString) / Double.Parse(kges.ToString)) * 100
Me.Invoke(ChartDel3, Tx, x, y, j)
Me.Invoke(PBGel, 100, 100, 100, 100, PBValue5)
Me.Invoke(kDel, k)
TxListSektor3.RemoveAt(i)
xListSektor3.RemoveAt(i)
yListSektor3.RemoveAt(i)
i = i - 1
k = k + 1
j = j + 1
Loop Until j > a
a = TxListSektor4.Count - 1
j = 0
i = a
Do
Tx = Double.Parse(TxListSektor4.Item(i).ToString)
x = Double.Parse(xListSektor4.Item(i).ToString)
y = Double.Parse(yListSektor4.Item(i).ToString)
PBValue5 = (Double.Parse(k.ToString) / Double.Parse(kges.ToString)) * 100
Me.Invoke(ChartDel4, Tx, x, y, j)
Me.Invoke(PBGel, 100, 100, 100, 100, PBValue5)
Me.Invoke(kDel, k)
TxListSektor4.RemoveAt(i)
xListSektor4.RemoveAt(i)
yListSektor4.RemoveAt(i)
i = i - 1
k = k + 1
j = j + 1
Loop Until j > a
Me.Invoke(PBGel, 100, 100, 100, 100, 100)
End Sub
Private Sub Backgroundworker2_RunWorkerCompleted(sender As Object, e As System.ComponentModel.RunWorkerCompletedEventArgs) Handles BackgroundWorker2.RunWorkerCompleted
For x As Double = -((d + 2 * s) / 1000) To ((d + 2 * s) / 1000) Step 0.001
Chart1.Series("Rohr").Points.AddXY(x, Math.Sqrt((((d + 2 * s) / 1000) ^ 2) - (x ^ 2)))
Chart1.Series("Rohr").Points.AddXY(x, -Math.Sqrt((((d + 2 * s) / 1000) ^ 2) - (x ^ 2)))
Next
Chart1.ChartAreas("ChartArea1").AxisX.Minimum = 0
Chart1.ChartAreas("ChartArea1").AxisX.Maximum = 12
Chart1.ChartAreas("ChartArea1").AxisY.Minimum = z - 12
Chart1.ChartAreas("ChartArea1").AxisY.Maximum = z
Chart1.Visible = True
Timer1.Stop()
End Sub
Public Sub ProgressBarChange(ByVal PBValue1 As Double, ByVal PBValue2 As Double, ByVal PBValue3 As Double, ByVal PBValue4 As Double, ByVal PBValue5 As Double)
ProgressBar1.Value = PBValue1
ProgressBar2.Value = PBValue2
ProgressBar3.Value = PBValue3
ProgressBar4.Value = PBValue4
ProgressBar5.Value = PBValue5
End Sub
Public Sub ChartFill1(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Chart1.Series("Sektor1").Points.AddXY(x, y)
If Tx >= 45 Then
Chart1.Series("Sektor1").Points(j).Color = Color.Red
Else
If Tx >= 35 And Tx < 45 Then
Chart1.Series("Sektor1").Points(j).Color = Color.Violet
Else
Chart1.Series("Sektor1").Points(j).Color = Color.Blue
End If
End If
End Sub
Public Sub ChartFill2(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Chart1.Series("Sektor2").Points.AddXY(x, y)
If Tx >= 45 Then
Chart1.Series("Sektor2").Points(j).Color = Color.Red
Else
If Tx >= 35 And Tx < 45 Then
Chart1.Series("Sektor2").Points(j).Color = Color.Violet
Else
Chart1.Series("Sektor2").Points(j).Color = Color.Blue
End If
End If
End Sub
Public Sub ChartFill3(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Chart1.Series("Sektor3").Points.AddXY(x, y)
If Tx >= 45 Then
Chart1.Series("Sektor3").Points(j).Color = Color.Red
Else
If Tx >= 35 And Tx < 45 Then
Chart1.Series("Sektor3").Points(j).Color = Color.Violet
Else
Chart1.Series("Sektor3").Points(j).Color = Color.Blue
End If
End If
End Sub
Public Sub ChartFill4(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Chart1.Series("Sektor4").Points.AddXY(x, y)
If Tx >= 45 Then
Chart1.Series("Sektor4").Points(j).Color = Color.Red
Else
If Tx >= 35 And Tx < 45 Then
Chart1.Series("Sektor4").Points(j).Color = Color.Violet
Else
Chart1.Series("Sektor4").Points(j).Color = Color.Blue
End If
End If
End Sub
Public Sub Point(ByVal k As Int64)
Label12.Text = "Anzahl Datenpunkte: " & k
End Sub
Private Sub Timer1_Tick(sender As Object, e As EventArgs) Handles Timer1.Tick
Sek = Sek + 1
If Sek = 60 Then
Min = Min + 1
Sek = 0
End If
Label11.Text = "Berechnungszeit: " & Min & " Min. " & Sek & " Sek."
End Sub
End Class

und hier die geänderte variante mit den farbweisen series:

Spoiler anzeigen

VB.NET-Quellcode

Public Class Form1
Dim d As Double
Dim TM As Double
Dim s As Double
Dim Lambda As Double
Dim z As Double
Dim z1 As Double
Dim q As Double
Dim T0 As Double
Dim Tx As Double
Dim x As Double
Dim y As Double
Dim i As Double
Dim j As Int64
Dim h As Double
Dim LR As Double
Dim LE As Double
Dim alpha As Double
Dim Tz As Double
Dim TxListSektor1 As New ArrayList
Dim TxListSektor2 As New ArrayList
Dim TxListSektor3 As New ArrayList
Dim TxListSektor4 As New ArrayList
Dim xListSektor1 As New ArrayList
Dim xListSektor2 As New ArrayList
Dim xListSektor3 As New ArrayList
Dim xListSektor4 As New ArrayList
Dim yListSektor1 As New ArrayList
Dim yListSektor2 As New ArrayList
Dim yListSektor3 As New ArrayList
Dim yListSektor4 As New ArrayList
Dim PBValue1 As Double
Dim PBValue2 As Double
Dim PBValue3 As Double
Dim PBValue4 As Double
Dim PBValue5 As Double
Delegate Sub ProgressBar(ByVal PBValue1 As Double, ByVal PBValue2 As Double, ByVal PBValue3 As Double, ByVal PBValue4 As Double, ByVal PBValue5 As Double)
Delegate Sub Charting1(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Delegate Sub Charting2(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Delegate Sub Charting3(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Delegate Sub Charting4(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
Delegate Sub Points(ByVal k As Int64)
Dim Min As Integer = 0
Dim Sek As Integer = 0
Public Shared Function arcosh(ByVal x As Double) As Double
arcosh = Math.Log(x + Math.Sqrt((x ^ 2) - 1))
End Function
Private Sub Form_Load(sender As Object, e As EventArgs) Handles MyBase.Load
Label11.Text = "Berechnungszeit: " & Min & " Min. " & Sek & " Sek."
End Sub
Private Sub Button2_Click(sender As Object, e As EventArgs) Handles Button2.Click
Application.Exit()
End Sub
Private Sub Button1_Click(sender As Object, e As EventArgs) Handles Button1.Click
d = Double.Parse(TextBox1.Text)
TM = Double.Parse(TextBox2.Text)
s = Double.Parse(TextBox3.Text)
Lambda = Double.Parse(TextBox4.Text)
z = Double.Parse(TextBox5.Text)
Chart1.Series.Add("Rohr")
Chart1.Series.Add("Red")
Chart1.Series.Add("Violet")
Chart1.Series.Add("Blue")
Chart1.Series("Rohr").ChartType = DataVisualization.Charting.SeriesChartType.Line
Chart1.Series("Red").ChartType = DataVisualization.Charting.SeriesChartType.Point
Chart1.Series("Violet").ChartType = DataVisualization.Charting.SeriesChartType.Point
Chart1.Series("Blue").ChartType = DataVisualization.Charting.SeriesChartType.Point
Chart1.Visible = False
Timer1.Start()
BackgroundWorker1.RunWorkerAsync()
End Sub
Private Sub Backgroundworker1_DoWork(sender As Object, e As System.ComponentModel.DoWorkEventArgs) Handles BackgroundWorker1.DoWork
Dim PBGel As New ProgressBar(AddressOf Me.ProgressBarChange)
LR = Math.Log((d + 2 * s) / d) / (2 * Lambda * Math.PI)
j = 0
For Me.alpha = 0 To ((180 * Math.Atan(12 / z)) / Math.PI) Step 0.1
i = 0
z1 = z / Math.Cos(alpha * (Math.PI / 180))
LE = arcosh((2 * z1) / ((d + 2 * s) / 1000)) / (4 * Math.PI)
q = (TM - 104.1) / (LR + LE)
T0 = 14.1 + q * LE
Do
Tx = T0 - i * ((T0 - 14.1) / (z1 - (d / 2000) - (s / 1000)))
x = ((d / 2000) + (s / 1000) + i) * Math.Sin(alpha * (Math.PI / 180))
y = ((d / 2000) + (s / 1000) + i) * Math.Cos(alpha * (Math.PI / 180))
TxListSektor1.Add(Tx)
xListSektor1.Add(x)
yListSektor1.Add(y)
i = i + 0.001
j = j + 1
Loop Until i > (z1 - (d / 2000) - (s / 1000))
PBValue1 = (alpha / ((180 * Math.Atan(12 / z)) / Math.PI)) * 100
Me.Invoke(PBGel, PBValue1, 0, 0, 0, 0)
Next
PBValue1 = 100
Me.Invoke(PBGel, 100, 0, 0, 0, 0)
For Me.alpha = 0 To ((180 * Math.Atan(z / 12)) / Math.PI) Step 0.1
i = 0
z1 = 12 / Math.Cos(alpha * (Math.PI / 180))
LE = arcosh((2 * z1) / ((d + 2 * s) / 1000)) / (4 * Math.PI)
h = 12 * Math.Tan(alpha * (Math.PI / 180))
Tz = 9.78 + ((14.1 - 9.78) * Math.Exp((-h) / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))) * Math.Cos(h / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000))))) - (13.16 * Math.Exp((-h) / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))) * Math.Sin(h / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))))
q = (TM - Tz) / (LR + LE)
T0 = Tz + q * LE
Do
Tx = T0 - i * ((T0 - Tz) / (z1 - (d / 2000) - (s / 1000)))
x = ((d / 2000) + (s / 1000) + i) * Math.Cos(alpha * (Math.PI / 180))
y = ((d / 2000) + (s / 1000) + i) * Math.Sin(alpha * (Math.PI / 180))
TxListSektor2.Add(Tx)
xListSektor2.Add(x)
yListSektor2.Add(y)
i = i + 0.001
j = j + 1
Loop Until i > (z1 - (d / 2000) - (s / 1000))
PBValue2 = (alpha / ((180 * Math.Atan(z / 12)) / Math.PI)) * 100
Me.Invoke(PBGel, 100, PBValue2, 0, 0, 0)
Next
PBValue2 = 100
Me.Invoke(PBGel, 100, 100, 0, 0, 0)
For Me.alpha = 0 To ((180 * Math.Atan((12 - z) / 12)) / Math.PI) Step 0.1
i = 0
z1 = 12 / Math.Cos(alpha * (Math.PI / 180))
LE = arcosh((2 * z1) / ((d + 2 * s) / 1000)) / (4 * Math.PI)
h = 12 * Math.Tan(alpha * (Math.PI / 180)) + z
Tz = 9.78 + ((14.1 - 9.78) * Math.Exp((0 - h) / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))) * Math.Cos(h / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000))))) - (13.16 * Math.Exp((0 - h) / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))) * Math.Sin(h / Math.Sqrt((2 * (2 / (1500 * 2000))) / ((2 * Math.PI) / (4 * 2592000)))))
q = (TM - Tz) / (LR + LE)
T0 = Tz + q * LE
Do
Tx = T0 - i * ((T0 - Tz) / (z1 - (d / 2000) - (s / 1000)))
x = ((d / 2000) + (s / 1000) + i) * Math.Cos(alpha * (Math.PI / 180))
y = 0 - (((d / 2000) + (s / 1000) + i) * Math.Sin(alpha * (Math.PI / 180)))
TxListSektor3.Add(Tx)
xListSektor3.Add(x)
yListSektor3.Add(y)
i = i + 0.001
j = j + 1
Loop Until i > (z1 - (d / 2000) - (s / 1000))
PBValue3 = (alpha / ((180 * Math.Atan((12 - z) / 12)) / Math.PI)) * 100
Me.Invoke(PBGel, 100, 100, PBValue3, 0, 0)
Next
PBValue3 = 100
Me.Invoke(PBGel, 100, 100, 100, 0, 0)
For Me.alpha = 0 To (180 * (Math.Atan(12 / (12 - z)) / Math.PI)) Step 0.1
i = 0
z1 = (12 - z) / Math.Cos(alpha * (Math.PI / 180))
LE = arcosh((2 * z1) / ((d + 2 * s) / 1000)) / (4 * Math.PI)
q = (TM - 9.78) / (LR + LE)
T0 = 9.78 + q * LE
Do
Tx = T0 - i * ((T0 - 9.78) / (z1 - (d / 2000) - (s / 1000)))
x = ((d / 2000) + (s / 1000) + i) * Math.Sin(alpha * (Math.PI / 180))
y = 0 - (((d / 2000) + (s / 1000) + i) * Math.Cos(alpha * (Math.PI / 180)))
TxListSektor4.Add(Tx)
xListSektor4.Add(x)
yListSektor4.Add(y)
i = i + 0.001
j = j + 1
Loop Until i > (z1 - (d / 2000) - (s / 1000))
PBValue4 = (alpha / (180 * (Math.Atan(12 / (12 - z)) / Math.PI))) * 100
Me.Invoke(PBGel, 100, 100, 100, PBValue4, 0)
Next
PBValue4 = 100
Me.Invoke(PBGel, 100, 100, 100, 100, 0)
End Sub
Private Sub Backgroundworker1_RunWorkerCompleted(sender As Object, e As System.ComponentModel.RunWorkerCompletedEventArgs) Handles BackgroundWorker1.RunWorkerCompleted
BackgroundWorker2.RunWorkerAsync()
End Sub
Private Sub Backgroundworker2_DoWork(sender As Object, e As System.ComponentModel.DoWorkEventArgs) Handles BackgroundWorker2.DoWork
Dim ChartDel1 As New Charting1(AddressOf ChartFill1)
Dim ChartDel2 As New Charting2(AddressOf ChartFill2)
Dim ChartDel3 As New Charting3(AddressOf ChartFill3)
Dim ChartDel4 As New Charting4(AddressOf ChartFill4)
Dim PBGel As New ProgressBar(AddressOf Me.ProgressBarChange)
Dim kDel As New Points(AddressOf Me.Point)
Dim kges As Int64 = TxListSektor1.Count + TxListSektor2.Count + TxListSektor3.Count + TxListSektor4.Count
Dim k As Int64 = 1
Dim a As Int64
a = TxListSektor1.Count - 1
j = 0
i = a
Do
Tx = Double.Parse(TxListSektor1.Item(i).ToString)
x = Double.Parse(xListSektor1.Item(i).ToString)
y = Double.Parse(yListSektor1.Item(i).ToString)
PBValue5 = (Double.Parse(k.ToString) / Double.Parse(kges.ToString)) * 100
Me.Invoke(ChartDel1, Tx, x, y, j)
Me.Invoke(PBGel, 100, 100, 100, 100, PBValue5)
Me.Invoke(kDel, k)
TxListSektor1.RemoveAt(i)
xListSektor1.RemoveAt(i)
yListSektor1.RemoveAt(i)
i = i - 1
k = k + 1
j = j + 1
Loop Until j > a
a = TxListSektor2.Count - 1
j = 0
i = a
Do
Tx = Double.Parse(TxListSektor2.Item(i).ToString)
x = Double.Parse(xListSektor2.Item(i).ToString)
y = Double.Parse(yListSektor2.Item(i).ToString)
PBValue5 = (Double.Parse(k.ToString) / Double.Parse(kges.ToString)) * 100
Me.Invoke(ChartDel2, Tx, x, y, j)
Me.Invoke(PBGel, 100, 100, 100, 100, PBValue5)
Me.Invoke(kDel, k)
TxListSektor2.RemoveAt(i)
xListSektor2.RemoveAt(i)
yListSektor2.RemoveAt(i)
i = i - 1
k = k + 1
j = j + 1
Loop Until j > a
a = TxListSektor3.Count - 1
j = 0
i = a
Do
Tx = Double.Parse(TxListSektor3.Item(i).ToString)
x = Double.Parse(xListSektor3.Item(i).ToString)
y = Double.Parse(yListSektor3.Item(i).ToString)
PBValue5 = (Double.Parse(k.ToString) / Double.Parse(kges.ToString)) * 100
Me.Invoke(ChartDel3, Tx, x, y, j)
Me.Invoke(PBGel, 100, 100, 100, 100, PBValue5)
Me.Invoke(kDel, k)
TxListSektor3.RemoveAt(i)
xListSektor3.RemoveAt(i)
yListSektor3.RemoveAt(i)
i = i - 1
k = k + 1
j = j + 1
Loop Until j > a
a = TxListSektor4.Count - 1
j = 0
i = a
Do
Tx = Double.Parse(TxListSektor4.Item(i).ToString)
x = Double.Parse(xListSektor4.Item(i).ToString)
y = Double.Parse(yListSektor4.Item(i).ToString)
PBValue5 = (Double.Parse(k.ToString) / Double.Parse(kges.ToString)) * 100
Me.Invoke(ChartDel4, Tx, x, y, j)
Me.Invoke(PBGel, 100, 100, 100, 100, PBValue5)
Me.Invoke(kDel, k)
TxListSektor4.RemoveAt(i)
xListSektor4.RemoveAt(i)
yListSektor4.RemoveAt(i)
i = i - 1
k = k + 1
j = j + 1
Loop Until j > a
Me.Invoke(PBGel, 100, 100, 100, 100, 100)
End Sub
Private Sub Backgroundworker2_RunWorkerCompleted(sender As Object, e As System.ComponentModel.RunWorkerCompletedEventArgs) Handles BackgroundWorker2.RunWorkerCompleted
For x As Double = -((d + 2 * s) / 1000) To ((d + 2 * s) / 1000) Step 0.001
Chart1.Series("Rohr").Points.AddXY(x, Math.Sqrt((((d + 2 * s) / 1000) ^ 2) - (x ^ 2)))
Chart1.Series("Rohr").Points.AddXY(x, -Math.Sqrt((((d + 2 * s) / 1000) ^ 2) - (x ^ 2)))
Next
Chart1.Series("Red").Color = Color.Red
Chart1.Series("Violet").Color = Color.Violet
Chart1.Series("Blue").Color = Color.Blue
Chart1.ChartAreas("ChartArea1").AxisX.Minimum = 0
Chart1.ChartAreas("ChartArea1").AxisX.Maximum = 12
Chart1.ChartAreas("ChartArea1").AxisY.Minimum = z - 12
Chart1.ChartAreas("ChartArea1").AxisY.Maximum = z
Chart1.Visible = True
Timer1.Stop()
End Sub
Public Sub ProgressBarChange(ByVal PBValue1 As Double, ByVal PBValue2 As Double, ByVal PBValue3 As Double, ByVal PBValue4 As Double, ByVal PBValue5 As Double)
ProgressBar1.Value = PBValue1
ProgressBar2.Value = PBValue2
ProgressBar3.Value = PBValue3
ProgressBar4.Value = PBValue4
ProgressBar5.Value = PBValue5
End Sub
Public Sub ChartFill1(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
If Tx >= 45 Then
Chart1.Series("Red").Points.AddXY(x, y)
Else
If Tx >= 35 And Tx < 45 Then
Chart1.Series("Violet").Points.AddXY(x, y)
Else
Chart1.Series("Blue").Points.AddXY(x, y)
End If
End If
End Sub
Public Sub ChartFill2(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
If Tx >= 45 Then
Chart1.Series("Red").Points.AddXY(x, y)
Else
If Tx >= 35 And Tx < 45 Then
Chart1.Series("Violet").Points.AddXY(x, y)
Else
Chart1.Series("Blue").Points.AddXY(x, y)
End If
End If
End Sub
Public Sub ChartFill3(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
If Tx >= 45 Then
Chart1.Series("Red").Points.AddXY(x, y)
Else
If Tx >= 35 And Tx < 45 Then
Chart1.Series("Violet").Points.AddXY(x, y)
Else
Chart1.Series("Blue").Points.AddXY(x, y)
End If
End If
End Sub
Public Sub ChartFill4(ByVal Tx As Double, ByVal x As Double, ByVal y As Double, ByVal j As Int64)
If Tx >= 45 Then
Chart1.Series("Red").Points.AddXY(x, y)
Else
If Tx >= 35 And Tx < 45 Then
Chart1.Series("Violet").Points.AddXY(x, y)
Else
Chart1.Series("Blue").Points.AddXY(x, y)
End If
End If
End Sub
Public Sub Point(ByVal k As Int64)
Label12.Text = "Anzahl Datenpunkte: " & k
End Sub
Private Sub Timer1_Tick(sender As Object, e As EventArgs) Handles Timer1.Tick
Sek = Sek + 1
If Sek = 60 Then
Min = Min + 1
Sek = 0
End If
Label11.Text = "Berechnungszeit: " & Min & " Min. " & Sek & " Sek."
End Sub
End Class

RodFromGermany · 21. Dezember 2013, 16:41

affrop schrieb:

das kommt so nach etwa 3,2 mio. datenpunkten

Das ist offensichtlich zu viel.
Um einen Algorithmus zu entwickeln / zu testen nimm 30 Datensätze.
Wenn alles läuft, versuch Dich an mehr Daten.
Überlege, wieviel Pixel Dein Bildschirm breit und hoch ist, überlege, wie breit / hoch Dein ChartControl ist und überleg Dir dann, wieviel Datensätze überhaupt sinnvoll dargestellt werden können.
Das Zauberwort heißt Datenvorverarbeitung.

affrop · 21. Dezember 2013, 21:06

das problem liegt in der mathematik dahinter. stichwort finite-elemente-rechnung.
das rechenbeispiel, was ich verwende kommt am ende etwa auf 20mio datenpunkte. das ist noch extrem wenig für ne fer.
am ende muss eine darstellung rauskommen, die vollflächig farbverläufe darstellt. algorithmisch die daten zu filtern wird nahezu unmöglich, sonst gäbe es keine über den gesamten wertebereich der berechnung gültige lösung.
das chart ist zu testzwecken erstmal fix auf 675x526px gesetzt. wenns dann mal irgendwann funktioniert soll der anwender die möglichkeit haben die form zu maximieren und was der dann für einen monitor oder grafikeinstellungen hat können nicht abgeschätzt werden. deshalb wäre es gut 20mio datenpunkte und mehr in eine chart-series prügeln zu können. es gibt ja auf der basis von cad diverse software die solche finite-elemente-rechnungen durchführt. also muss es ja irgendwie gehen diese vielzahl von datenpunkten in ein chart oder was ähnliches zu bringen.

ErfinderDesRades · 21. Dezember 2013, 21:27

also ich kenne Charts nicht so, dass man da vollflächig ausmalt - millionen Punkte in allen möglichen Farben.

Charts verstehe ich so, dass paar (< 100) Datenpunkte aufgenommen werden, und dann werden Linien dazwischen interpoliert.

Bist du sicher, dass MS-ChartControl für dich das richtige ist?

affrop · 21. Dezember 2013, 22:34

ich bin mir nicht sicher ob das chart das richtige control ist, aber auf jedenfall dass einzige was mir auf die schnelle einviel. ich bin gerne offen für alternative vorschläge.
fakt ist jedenfalls, dass durch die finite-elemente-rechnung exakte datenpunkte durch mehrdimensionale matrizen ermittelt werden, die wesentlich für die analytische lösung sind. im zweifel müsste das chart zwischen jedem gebenenen datenpunkt eine linie zeichnen und auf der interpolieren.
die anzahl der datenpunkte ist wesentlich durch die berechnung selbst definiert und nicht durch die eingabewerte. insofern kann mann daran nicht viel ändern. und da das ganze farbverläufe darstellen soll (sonst wäre es ja kein temperaturfeld) kanns unter umständen sein, dass jeder datenpunkt ne eigene farbe hat, je nach zugehöriger temperaur.

ErfinderDesRades · 21. Dezember 2013, 22:42

jetzt habichmal gewiki-t, und das sieht mir doch sehr nach olle Bitmap aus - da magst du jedes pixel berechnen, und angezeigt wirds inne Picturebox.

vb-checker · 21. Dezember 2013, 22:48

zeichne dir doch erstmal alle Punkte in ein Usercontrol mittels GDI+. Das Control sollte die Punkte zum Beispiel in einem Dictionary verwalten (Location (Typ Point) und Value (Typ Color)). Du zeichnest die Punkte auf eine Bitmap, die du nur neu erstellst, wenn ein neuer Punkt hinzugefügt wird (oder sogar seltener). Im Paint Event deines Controls zeichnest du diese Bitmap, ggf verkleinert. So hast du alle Freiheiten und bist nicht durch das Chartcontrol limitiert.

affrop · 21. Dezember 2013, 22:52

@ erfinder: die frage ist nur, wie die punkte in die picbox bringen, denn die berechneten datenpunkte sind x-y-relativ, also kartesisch. in ner picbox gibts ja kein koordinatensystem.

@vb-checker: gibts da schon vorgefertigte controls oder müsste so ein ding selbst deklariert werden?

ErfinderDesRades · 21. Dezember 2013, 22:54

affrop schrieb:

die frage ist nur, wie die punkte in die picbox bringen

genau lesen:
Die punkte sind in einer Bitmap zu berechnen und zu setzen.

gibts da schon vorgefertigte controls

PictureBox

affrop · 21. Dezember 2013, 22:58

ich glaub du hast mich falsch verstanden. die brechnung gibt x-y-koordinaten aus, die in ein kartesisches koordinatensystem gezeichnet werden müssen. das koordinatensystem liefert das chart mit. eine bitmap oder eine picbox nicht.

00yoshi · 21. Dezember 2013, 23:12

VB.NET-Quellcode

Dim bmp as new bitmap(20000,1000)
Dim col as color = color.fromargb(alpha,rot,grün,blau)
bmp.setpixel(x,y,col)
picturebox1.image = bmp

affrop · 21. Dezember 2013, 23:45

@00yoshi: welchen bezuigpunkt haben deine x-y-koordinaten?
die datenpunkte, die in meinem programm berechnet werden, beziehen sich auf ein standardisiertes kartesisches koordinatensystem dessen mittelpunkt (0;0) irgendwo liegen kann. währe dann also der mittelpuinkt des koordinatensystems in pixel in der bitmap die obere linke ecke wie in allen andren controls oder ist das variabel?

RodFromGermany · 22. Dezember 2013, 09:09

affrop schrieb:

finite-elemente-rechnung.

FEM
Solch unwichtige Information gehört in den 1. Post!
Da ist ein ChartControl völlig daneben, da brauchst Du was dynamisches, was mit Hardwarebeschleunigung. DirectX & Co oder WPF.

Higlav · 22. Dezember 2013, 10:34

Juhuu - FEM!

Aber dann greif lieber zu MatLab, wenn du das hast!

Temperaturfeld in Chart darstellen

Temperaturfeld in Chart darstellen

VB.NET-Quellcode

affrop schrieb:

VB.NET-Quellcode

VB.NET-Quellcode

affrop schrieb:

Skybird schrieb:

affrop schrieb:

VB.NET-Quellcode

affrop schrieb:

Ähnliche Themen

5 Benutzer haben hier geschrieben