---

Hallo allerseits!
Ich bin momentan mit Beulanalysen von Laminaten beschäftigt. Dazu verwende ich Shell181-Elemente mit Section-Befehlen.
Bei der Konvergenzstudie bin ich stutzig geworden, dass die Beullast (signifikant) von der Netzfeinheit abhängt. Diese HALBIERT sich bei einer sehr feinen Elementierung.

Gibt es eine Erklärung dafür??

Das gleiche tritt beispielsweise auch bei dem VM-82 Beispiel der Ansys-Hilfe auf:

FINISH
/CLEAR,NOSTART
/PREP7
SMRT,OFF
/TITLE, VM82, SIMPLY SUPPORTED LAMINATED PLATE UNDER PRESSURE
C*** USING SHELL181
/PREP7
SMRT,OFF
ANTYPE,STATIC
ET,1,SHELL181              ! 4 NODE LAYERED SHELL
KEYOPT,1,3,2                ! FULL INTEGRATION
KEYOPT,1,8,1                ! WRITE LAYER RESULTS
SECTYPE,1,SHELL
SECDATA,0.025,1,0          ! LAYER 1: 0.025 THK, THETA 0
SECDATA,0.025,1,90          ! LAYER 2: 0.025 THK, THETA 90
SECDATA,0.025,1,90          ! LAYER 3: 0.025 THK, THETA 90
SECDATA,0.025,1,0          ! LAYER 4: 0.025 THK, THETA 0
MP,EX,1,25E6                ! ORTHOTROPIC MATERIAL PROPERTIES
MP,EY,1,1E6
MP,EZ,1,1E6                ! EZ=EY ASSUMED
MP,GXY,1,5E5
MP,GYZ,1,2E5
MP,GXZ,1,5E5
MP,PRXY,1,0.25              ! MAJOR POISSONS RATIO
MP,PRYZ,1,0.01              ! MAJOR POISSONS RATIO
MP,PRXZ,1,0.25              ! MAJOR POISSONS RATIO
K,1                        ! CORNER KEYPOINTS OF QUADRANT (AREA)
K,2,5
K,3,5,5
K,4,,5
A,1,2,3,4
ESIZE,,150                  ! (Original: 6)  Sprung zwischen 150 und 200!??!
AMESH,1
NSEL,S,LOC,X,0
DSYM,SYMM,X
NSEL,S,LOC,Y,0
DSYM,SYMM,Y
NSEL,S,LOC,X,5              ! APPLY FREELY SUPPORTED B.C.
D,ALL,UZ,,,,,UY
NSEL,S,LOC,Y,5
D,ALL,UZ,,,,,UX
NSEL,ALL
SFE,ALL,2,PRES,,1          ! APPLY UNIFORM PRESSURE
OUTPR,NSOL,1
OUTPR,RSOL,1
FINISH
/SOLU   
SOLVE
FINISH
/POST1
NSEL,S,LOC,X                ! SELECT CENTER NODE
NSEL,R,LOC,Y
PRNSOL,U,Z                  ! PRINT CENTER DEFLECTION
*GET,DEF99,NODE,1,U,Z
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'DEF,m (S'
LABEL(1,2) = 'HELL181)'
*VFILL,VALUE(1,1),DATA,.0683
*VFILL,VALUE(1,2),DATA,ABS(DEF99)
*VFILL,VALUE(1,3),DATA,ABS(DEF99/.0683)
SAVE,INF4
FINISH

Wird dort eine Unterteilung von 200 vorgenommen, so halbiert sich die Durchsenkung der Platte ebenfalls (0,068m vs. 0,03m).

Kann mir jemand dafür eine schlüssige Erklärung geben?

Viele Grüße

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

CADFEM liefert die Antwort:

"man kann ein Finites Elemente Netz nicht unendlich verfeinern, da es genau zu so einem Effekt kommt den Sie hier beobachten.
Theoretisch ist es so, dass je feiner man vernetzt, umso genauer wird das Ergebnis. Leider macht uns ab einer gewissen Verfeinerung die Nummeric einen Strich durch die Rechnung. Die Integrationspunkte der Schalen liegen immer näher beieinander, bis sie sich fast in einem Punkt treffen. Die Steifigkeitsmatrix wird dann immer steifer.

Eine Elementierung von 200 ist in diesem Beispiel einfach nicht sinnvoll.

Es ist wie im Leben: Viel hilft nicht immer viel.

Mit freundlichen Grüßen

CADFEM GmbH"

Vielen Dank!!

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP