请教关于膜结构风荷载的建立

您好：

如何建立施加在不规则曲面上的风荷载呢?如下图：对于一个圆形的界面，横向方向的风，会在曲面的不同位置产生不同的荷载。迎风面为正风压，背风面为负风压，也就是部分的真空。不知道sofistik会自动进行风载的计算，还是需要自行手动计算然后分区域施加呢？

我的代码如下，在代码最后一行有施加的风荷载的代码，总是报错。大神可以帮忙看一下吗？

+PROG AQUA urs:5.1

KOPF Material $Benenung des Moduls

ECHO MAT,QUER VOLL $ liefert Bilder in Ursula

mate 1 e 1200 mue 0.4 g 30 gam 0.005 e90 900 bez 'Textilmembrane'

nmat 1 art memb $Nichtlineares Material Nr. 1 ist eine Membran

ENDE

+PROG SOFIMSHA urs:5.2 $ Text Eingabe zur Systemgenerierung

KOPF Text Eingabe Systemgenerierung

SYST init

syst raum gdiv 10000

$Variablen definieren

let#D_Oeff 1.0 $Durchmesser der ?ffnung am Hochpunkt in [m]

let#l 10 $L?nge der ?u?eren Kanten in [m]

let#teil 10 $Teilung der Kanten

$Knoten f?r den inneren Kreis festlegen

knot 1 0 0 0 FIX F $Mittelpunktsknoten, komplett festgehalten

knot 2 #D_Oeff/2 0 0 phiz 0  fix kf nr1 1 $gekoppelt an Festhaltungen Knoten 1

knot 3 #D_Oeff/2 0 0 phiz 90 fix kf nr1 1

knot 4 #D_Oeff/2 0 0 phiz 180  fix kf nr1 1

knot 5 #D_Oeff/2 0 0 phiz 270  fix kf nr1 1

knot 23 #D_Oeff/2 0 0 phiz 45  fix kf nr1 1

knot 34 #D_Oeff/2 0 0 phiz 135 fix kf nr1 1

knot 45 #D_Oeff/2 0 0 phiz 225 fix kf nr1 1

knot 52 #D_Oeff/2 0 0 phiz 315 fix kf nr1 1

$R?nder des inneren Kreises

rand 1 bez '1.Viertel i'; rand ka 23 ke 34 km -3 teil #teil $typ fmb

rand 2 bez '2.Viertel i'; rand ka 34 ke 45 km -4 teil #teil $typ fmb

rand 3 bez '3.Viertel i'; rand ka 45 ke 52 km -5 teil #teil $typ fmb

rand 4 bez '4.Viertel i'; rand ka 52 ke 23 km -2 teil #teil $typ fmb

$Knoten des ?u?eren Quadrats

knot 100 #l/2 #l/2

knot 101 -#l/2 #l/2

knot 102 -#l/2 -#l/2

knot 103 #l/2 -#l/2

$R?nder des ?u?eren Quadrats

rand 100 bez 'Au?enrand 1'; rand 100 101 teil #teil typ fmb

rand 101 bez 'Au?enrand 2'; rand 101 102 teil #teil typ fmb

rand 102 bez 'Au?enrand 3'; rand 102 103 teil #teil typ fmb

rand 103 bez 'Au?enrand 4'; rand 103 100 teil #teil typ fmb

$Rand radial zur Elementeinteilung

rand 104 bez 'Rand radial'; rand ka 23 100 teil #teil*2  prog 1.15

rand 105 bez 'Rand radial'; rand ka 34 101 teil #teil*2  prog 1.15

rand 106 bez 'Rand radial'; rand ka 45 102 teil #teil*2  prog 1.15

rand 107 bez 'Rand radial'; rand ka 52 103 teil #teil*2  prog 1.15

$Elemente verlegen + automatisch vernetzen (fit)

grup 0 $Alle folgenden Elemente werden Gruppe 0 zugeordnet

quad fit k1 100 101 34 23 m34 3 mnr 1 t 1[mm] kr 1 $?ber kr 1 werden die lokalen x-Achsen der Elemente auf den Mittelknoten 1 ausgerichtet

quad fit k1 101 102 45 34 m34 4 mnr 1 t 1[mm] kr 1

quad fit k1 102 103 52 45 m34 5 mnr 1 t 1[mm] kr 1

quad fit k1 103 100 23 52 m34 2 mnr 1 t 1[mm] kr 1

let#KnotNr 1001  $Funktioniert nur bei #teil=10, sonst passt Knotennummerierung nicht

loop 32

knot #KnotNr fix kf nr1 1 $Kopplung der Festhaltung in Knoten 1

let#KnotNr #KnotNr+1

endloop

ENDE 

+PROG SOFILOAD urs:2

KOPF

SEIT UNIE 5 $ Einheiten: Querschnitte in mm, Geometrie+Lasten in m

LF 901 TYP NONE BEZ 'Auflagerverschiebung 1mm'

 KNOT 1 WZZ 1[mm]

ENDE

+prog ase urs:10.1

kopf formfindung

syst prob th3 iter 400 nmat nein tol 0.0001

grup - $Alle Gruppen werden aktiviert

grup 0 faks 1e-10 $Membranelemente ohne Steifigkeit in die Berechnung aufnehmen

hoch 0 0 vr1 7 vtvr 0.20 $Radiale Vorspannung 7 kN/m,tangentiale Vorsp. 20% der radialen Vorsp.

lf 1 egz 1.0 bez 'Formfindung' $Eigengewicht in z-Richtung 1.0-fach ber?cksichtigen

   LC 901 -4000  $ 4 m Hebung des zentralen Knotens 1

ende

+prog ase  urs:11.2

kopf Ausgleich der Knotenkr?fte

syst prob th3 iter 400 nmat ja tol 0.0001   plf 1

grup  faks 1

lf 2 bez 'Ausgleich der Knotenkr?fte'

ende

+PROG ASE URS:6

KOPF System updaten damit neue Verschiebungen vom

KOPF Formfindungungszustand LF 2 aus gerechnet werden

SYST PLF 2 spei Ja $Geometrie Update, Prim?rlastfall 2

ENDE

+prog sofiload urs:6.3

kopf windlast

LF 902 BEZ 'wind'

area ref qgrp nr 0 typ wind p1 1.0

LF 903

wind 1055 1 b dx +1.0

copy  902 1.0 wind

ende

+Prog ASE urs:12.1

kopf windlast

syst prob thII iter 800 nmat ja tol -.1   plf 2

grup faks 1

lf 3 bez 'windlasten'

LC 903

ende