SOFiSTiK Forum

It is currently Sun Jul 21, 2019 9:22 pm

All times are UTC




Post new topic Reply to topic  [ 2 posts ] 
Author Message
PostPosted: Tue Oct 10, 2017 3:09 pm 
Offline

Joined: Mon Jul 03, 2017 12:43 pm
Posts: 3
Schönen guten Tag!

ich versuche für meine Masterthesis einen Stab zu optimieren, der eine Membran aufspannt.
Das Optimierungsziel ist es, die Membran möglichst homogen in Spannung zu versetzen.
Nach Überlegung müsste demnach ein hyperbolischer Paraboloid herauskommen.

Den prinzipiellen Programmablauf verstehe ich, aber ich verzweifle etwas an meiner Zielfunktion und wie ich den Stab als variabel kennzeichne.

Code:
+PROG AQUA urs:1
KOPF Material und Querschnitte
SEIT UNIE 1
MATE 602 E  2.38/0.5 E90  14.1/0.5 MUE 0.3 GAM 0.00098/0.5  BEZ 'Membran'
NMAT 602 MEMB p1=0.0  p2=0.0
STAH NR 8 ART S  FY 235 FT 360 ES 210000 GAM 78.5
MATE NR 10 TYP PVC BEZ  'PVC  (Polyvinylchlorid)'  E 3000
ROHR 11 D 10 T 0 MNR 8 BC '0' bez 'Fundament'
ROHR 1 D 20 T 0 MNR 10 BC '0' bez 'PVC Stab 20'
ROHR 10 D 5 T 0 MNR 10 BC '0' bez 'PVC Stab 5'
ENDE

+PROG SOFIMSHC  urs:2 $ system definition
HEAD system definition
PAGE UNII 0
syst SPAC gdiv -1000 gdir NEGZ
CTRL TOPO 0
CTRL TOLG -0.01[m]
CTRL MESH 1
CTRL HMIN 0.07[m]
CTRL EFAC 2.68097
CTRL FINE -0.694
CTRL PROG 1.484

sar 1002 grp 1 mno 602 nra 6 qref 'CENT' t 0.0005[m] mctl 'auto' NX   0.000000  -1.000000   0.000000
 GUID ID "63dc3b49-fe58-4faa-ab71-8fd7de11bd09"
sarb OUT
 slnb x1   0.0000000000   0.0000000000   0.0000000000   x2   1.0000000000   0.0000000000   0.0000000000
sarb OUT
 slnb x1   1.0000000000   0.0000000000   0.0000000000   x2   0.2751834044   0.0000000000   0.4167772054
sarb OUT
 slnb x1   0.2751834044   0.0000000000   0.4167772054   x2   0.0000000000   0.0000000000   0.0000000000
 sarn s   0.0000000000 degs   1
 sarn s   0.8360994774
 sarn t  -1.0000000000 degt   1
 sarn t   0.0000000000
 sarp NURB    1    1    0.0000000000    0.0000000000    0.0000000000    1.0000000000
 sarp NURB    1    2    0.2751834044    0.0000000000    0.4167772054    1.0000000000
 sarp NURB    2    1    0.0000000000    0.0000000000    0.0000000000    1.0000000000
 sarp NURB    2    2    1.0000000000    0.0000000000    0.0000000000    1.0000000000
sln 1003 grp 1 sno '1.10' styp 'NE' kr 'NEGZ' drot 0.0
 slnb x1   0.0000000000   0.0000000000   0.0000000000  x2   0.2751834044   0.0000000000   0.4167772054
sln 1002 fix 'PXPYPZMXMYMZ' grp 1 sno 11 styp 'TE' kr 'NEGZ' drot 0.0 fixa 'MTMYMZMB' fixe 'MTMYMZMB' xflg 'PA'
 slnb x1   1.0000000000   0.0000000000   0.0000000000  x2   0.0000000000   0.0000000000   0.0000000000

end

+PROG AQUA  urs:3 $ interpolation of intermediate sections
HEAD interpolation of intermediate sections
PAGE UNII 0
INTE 0
END

+prog sofiload urs:4
kopf Aufstellen
        lf 20   bez 'Aufstellen'    fakt 1.0  egz -1.0
        knot (1001) typ dyy  p1 -450
ende

+prog ase urs:5
kopf  Aufstellen
syst prob th3 nmat ja TOL -0.2
STEU ITER 3 W2 1
STEU WARN 197
         lf (20)
ende

PROG TEMPLATE urs:6
        STO#quad_dst 2 ; STO#Rohr 1 D 20,5,20
ENDE
PROG AQUA ITER  urs:7
        RCL#Rohr
ENDE
PROG ASE ITER  urs:8
kopf  Aufstellen
syst prob th3 nmat ja TOL -0.2
STEU ITER 3 W2 1 STEU WARN 197
        lf  (20)
ende
@KEY QUAD_DST 2
LET# MAX
IF ABS(#MAX) < 0.001
        EXIT_ITERATION
ENDIF


Vielen Dank für eure Hilfe!

Mit freundlichen Grüßen
Daniel


Top
 Profile  
 
PostPosted: Thu Oct 12, 2017 3:04 pm 
Offline

Joined: Fri Jul 27, 2012 1:24 pm
Posts: 10
hi,

dein beispiel läuft leider nicht ohne fehlermeldung durch, programmabbruch im modul template:
"+++++ Fehler Nr. 10132 in Programm SOF_VAR
Ungültige Zahl/Ausdruck Typ 501 (unzulässiger Operator/CDB-Item-Name) :D"

des weiteren sind soweit ich gesehen habe alle knoten deines untenliegenden fachwerkstabes gehalten -> dann könntest du auf den verzichten.
wenn er sich nicht verformen kann, kriegt er ohnehin keine last.

auch ist mir anhand deiner beschreibung nicht klar, was optimiert werden soll -> die geometrie sprich die knotenkoordinaten deines aufspannstabes, oder soll die
querschnittsform angepasst werden ?

es gibt soweit ich weiß, eine interne möglichkeit in sofistik zur formfindung: viewtopic.php?f=18&t=547&view=next, das habe ich jedoch selbst noch nie verwendet.

lg


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 2 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 3 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group