Liebe Forumsbesucher,
ich versuche eine Transient Analyse eines Balkens zu machen, der an einem Tisch und einer Seitenleiste gelagert ist. Dabei wird der Balken mit einer zeitveränderlichen Kraft in positiver z-Richtung angeregt. Die Analyse wird für 0.5 Sekunden gemacht was lang genug sein sollte um den eingeschwungenen Spannungszustand zu erreichen.
Ich bekomme im während der Analyse sehr viele Warnungen (fehler), z.B.:
*** WARNING *** CP = 3.828 TIME= 04:31:01
Large deflection effects not included for element CONTA173.
*** WARNING *** CP = 4.234 TIME= 04:31:05
Coefficient ratio exceeds 1.0e8 - Check results.
*** WARNING *** CP = 4.406 TIME= 04:31:09
The calculated reference HT FLOW CONVERGENCE VALUE= 0 is less than the
default minimum of 1E-6 (MINREF field of CNVTOL command). Check
results carefully.
*** ERROR *** SUPPRESSED MESSAGE CP = 8376.016 TIME= 08:19:37
Solution not converged at time 1.18967856E-04 (load step 1 substep 1).
Run terminated.
*** WARNING *** SUPPRESSED MESSAGE CP = 8376.016 TIME= 08:19:37
The unconverged solution (identified as time 0.499902931 substep
999999) is output for analysis debug purposes. Results should not be
used for any other purpose.
/COM,ANSYS RELEASE 10.0 UP20050718 13:08:52 12/12/2010
*** WARNING *** SUPPRESSED MESSAGE CP = 13.609 TIME= 13:11:36
Cannot change number of variables with data stored.
The NUMV command is ignored.
Es ist offensichtlich, dass die Kontakt Parameter und der Code der Transient Analysis Fehler enthalten. Ich habe darin mittlerweile so oft Änderungen vorgenommen, dass ich die Fehler selbst nicht mehr finde. Darum bin ich hier auf der Suche nach Experten, die sich mal meinen Code anschauen können um die Fehler zu indentifizieren bevor ich die Netzdichte erhöhe und die Analyse weiter verfeinere.
Hier ist der Code:
/CLEAR
/FILNAME,TransAnalysis,1
/ANGLE,1,-90,XS,1 ! Rotates the display
/ANGLE,1,-150,YS,1
/REPLOT,FAST ! Fast replot
!*** DIMENSIONS ***
!* Desk, Steel
Swidth = 0.2
Sheight = 0.020
Sdepth = 1.195
!* Stop Rail, Steel
SRwidth = 0.020
SRheight = 0.020
SRdepth = 2.940
!* Workpiece, Wood
Wwidth = 0.1
Wheight = 0.1
Wdepth = 2
GGap = 0 ! Geometric gap, GGap=0 => contacts surfaces are closed
RGap = 1e-4 ! Gap between rigid node and rigid element
meshSize = 0.1
oH = 0.000 ! Overhang (Ueberstand) 0mm (for meshing purposes)
!* Position of rigid Node for Desk
rNodeDX = Sdepth/2
rNodeDY = Swidth/2
rNodeDZ = -Sheight - GGap - RGap
!* Position of rigid Node for Stop Rail
rNodeSRX = SRdepth/2
rNodeSRY = -SRwidth/2
rNodeSRZ = -RGap
!* Position of first driving roller
dR1x = 1.140
dR1z = Wheight
!* Position of second driving roller
dR2x = dR1x - 0.160
dR2z = Wheight
!* Position of side pressure roller
sPRx = 0.680
sPRy = Wwidth
sPRz = meshSize
FINISH
/PREP7
!*** GEOMETRY ***
!* Desk, Steel
BLOCK,0,Sdepth,0,Swidth,-GGap,-GGap-Sheight
!* Stop Rail, Steel
BLOCK,0,SRdepth,-GGap,-GGap-SRwidth,0,SRheight
!* Workpiece, Wood
BLOCK,Sdepth+oH,Sdepth+oH-Wdepth,0,Wwidth,0,Wheight
FINISH
!*** MATERIAL PROPERTIES AND REAL CONSTANTS ***
/PREP7
!* 1 * Desk, Steel
MP,EX,1,200e9
MP,NUXY,1,0.28
MP,DENS,1,0
R,1
!* 2 * Workpiece, Wood
MP,EX,2,1.3e10
MP,EY,2,5.0e8
MP,EZ,2,7.0e8
MP,GXY,2,6.32e8
MP,GYZ,2,4.70e8
MP,GXZ,2,2.22e8
MP,PRXY,2,0.5
MP,PRYZ,2,0.0133
MP,PRXZ,2,0.38
MP,DENS,2,0.47e3
R,2
!* 3 * Rigid Link Element
MP,EX,3,12e12
MP,PRXY,3,0.3
MP,DENS,3,0
R,3,1e5
!*** Element Types ***
!* 1 * Desk, Steel
ET,1,SOLID185
!* 2 * Workpiece, Wood
ET,2,SOLID185
!* 4 * Rigid Link Element
ET,3,LINK8
FINISH
! *** MESHING
/PREP7
SMRT,OFF
esize,meshSize
MSHAPE,0 ! 0: hexahedral-shaped elements; 1: triangle-shaped elements
MSHKEY,1 ! 0: free meshing; 1: mapped meshing
TYPE,1
REAL,1
MAT,1
VMESH,1
VMESH,2
TYPE,2
REAL,2
MAT,2
VMESH,3
FINISH
! *** RIGID REGIONS
/PREP7
*GET,nNr1,NODE,0,NUM,MAX
N,nNr1+1,rNodeDX,rNodeDY,rNodeDZ
N,nNr1+2,rNodeSRX,rNodeSRY,rNodeSRZ
D,nNr1+1,ALL
D,nNr1+2,ALL
!* Between Desk and Workpiece
TYPE,3
REAL,3
MAT,3
NSEL,S,LOC,X,0,Sdepth ! Select all the nodes of the desk
NSEL,R,LOC,Y,0,Swidth
NSEL,R,LOC,Z,-GGap,-GGap-Sheight
CM,nodesDesk,NODE
CMSEL,S,nodesDesk,NODE
*GET,NumDN,NODE,0,COUNT ! Retrieve the number of selected nodes, number of desk nodes
nTemp1 = 0
*DO,i,1,NumDN,1
CMSEL,S,nodesDesk,NODE
nTemp1 = NDNEXT(nTemp1) ! Get the Next Node Number of selected amount of nodes
E,nNr1+1,nTemp1
*ENDDO
ALLSEL
!* Between Stop Rail and Workpiece
TYPE,3
REAL,3
MAT,3
NSEL,S,LOC,X,0,SRdepth ! Select all the nodes of the stop rail
NSEL,R,LOC,Y,-GGap,-GGap-SRwidth
NSEL,R,LOC,Z,0,SRheight
CM,nodesSR,NODE
CMSEL,S,nodesSR,NODE
*GET,NumSRN,NODE,0,COUNT ! Retrieve the number of selected nodes, number of stop rail nodes
nTemp2 = 0
*DO,i,1,NumSRN,1
CMSEL,S,nodesSR,NODE
nTemp2 = NDNEXT(nTemp2) ! Get the Next Node Number of selected amount of nodes
E,nNr1+2,nTemp2
*ENDDO
ALLSEL
FINISH
! *** CONTACT AND TARGET ELEMENTS
/PREP7
!* Between Workpiece and Desk
NSEL,S,LOC,X,0,Sdepth+oH ! Select the nodes on the bottom of the workpiece, contact surface nodes (softer)
NSEL,R,LOC,Y,0,Wwidth
NSEL,R,LOC,Z,0
CM,nWPBottom,NODE
ALLSEL,ALL
NSEL,S,LOC,X,0,Sdepth ! Select the nodes on the surface of the desk, target surface nodes (stiffer)
NSEL,R,LOC,Y,0,Wwidth
NSEL,R,LOC,Z,-GGap
CM,nDeskSurface,NODE
ET,4,TARGE170
MP,MU,4,0.5
KEYOPT,4,1,0
KEYOPT,4,2,0
R,4,,,5,,1E-10
ET,5,CONTA173
KEYOPT,5,1,2
KEYOPT,5,2,1
KEYOPT,5,4,0
KEYOPT,5,5,3
KEYOPT,5,7,2
KEYOPT,5,9,1
KEYOPT,5,12,0
R,5,0,0,0.1,0.01,0,0.1
RMORE,0,0,1E20,0,1,1
RMORE,0,0,0,0,0,1
RMORE,0,0,1,1
TYPE,4 ! Target surface definition (Desk)
MAT,1
REAL,5
CMSEL,S,nDeskSurface,NODE
ESLN,S
ESURFACE
TYPE,5 ! Contact surface definition
MAT,2
REAL,5
CMSEL,S,nWPBottom,NODE
ESLN,S
ESURFACE
ALLSEL
FINISH
/PREP7
!* Between Workpiece and Stop Rail
NSEL,S,LOC,X,0,Sdepth+oH ! Select the nodes on the side of the workpiece, contact surface nodes (softer)
NSEL,R,LOC,Y,0
NSEL,R,LOC,Z,0,SRheight
CM,nWPSide,NODE
NSEL,S,LOC,X,0,Sdepth+oH ! Select the nodes on the side of the stop rail coincident to the workpiece, target surface nodes (stiffer)
NSEL,R,LOC,Y,-GGap
NSEL,R,LOC,Z,0,SRheight
CM,nSRSide,NODE
ET,6,TARGE170
MP,MU,6,0.5
KEYOPT,6,1,0
KEYOPT,6,2,0
R,6,,,5,,1E-10
ET,7,CONTA173
KEYOPT,7,1,2
KEYOPT,7,2,1
KEYOPT,7,4,0
KEYOPT,7,5,3
KEYOPT,7,7,2
KEYOPT,7,9,1
KEYOPT,7,12,0
R,7,0,0,0.1,0.01,0,0.1
RMORE,0,0,1E20,0,1,1
RMORE,0,0,0,0,0,1
RMORE,0,0,1,1
TYPE,6 ! Target surface definition (Stop Rail)
MAT,1
REAL,7
CMSEL,S,nSRSide,NODE
ESLN,S
ESURFACE
TYPE,7 ! Contact surface definition
MAT,2
REAL,7
CMSEL,S,nWPSide,NODE
ESLN,S
ESURFACE
ALLSEL
FINISH
! *** SELECTION OF NODES RELATED TO DRIVING AND PRESSURE ROLLERS
/PREP7
NSEL,S,LOC,X,1.095
NSEL,R,LOC,Y,0,0.1
NSEL,R,LOC,Z,0.1
CM,nDR1,NODE
NSEL,S,LOC,X,0.995
NSEL,R,LOC,Y,0,0.1
NSEL,R,LOC,Z,0.1
CM,nDR2,NODE
NSEL,S,LOC,X,0.695
NSEL,R,LOC,Y,0.1
NSEL,R,LOC,Z,0,0.1
CM,nSPR,NODE
ALLSEL
! **********************************************************
/CONFIG,NRES,1000000
/nerr,,99999999
/SOLUTION
/TITLE,Transient Analysis of a clamped wooden beam in a planing machine
ANTYPE,TRANSIENT,NEW
TRNOPT,FULL
!* DOF Constraints through Driving and Pressure Rollers
CMSEL,S,nDR1,NODE
D,ALL,UZ,,,,,UX ! UX included to ensure good convergence
CMSEL,S,nDR2,NODE
D,ALL,UZ,,,,,UX
CMSEL,S,nSPR,NODE
D,ALL,UY,,,,,UX
!* Weight
ACEL, , ,-9810
NLGEOM,ON ! Geometric Nonlinearities
OUTRES,ALL,ALL
f = 420.29 ! Frequency of trigger function (estimatation of eigenfrequency)
pi = 3,14159265
T = (2*pi)/f ! Period
ITS = 1/(20*f)
DELTIM,ITS
AUTOTS,ON ! Automatic Time Stepping is on
CycLeng = 0.5 ! Cycle length: Simulates exactly for 0.5 seconds
nFAP = CycLeng/ITS ! Number of total force application points
FZmax = 1000 ! 1000 N
F_Ax = -1000 ! -1000 N (Axial force)
IniT = 1e-5 ! Initial time when trigger force starts acting (convergence criterion: must be > 0)
TM_START = IniT
TM_END = CycLeng + IniT ! Ending time of the transient
TM_INCR = ITS ! Time increment
*DIM,TFORCE,TABLE,nFAP+1 ! Defining a (nFAP+1)x1x1 table
*DO,i,1,nFAP+1,1
t = ((i-1)*T/20)+IniT
TFORCE(i,1) = FZmax*sin(f*t-IniT) ! Force values in column 1
TFORCE(i,0) = t ! Corresponding time values in column 0
*ENDDO
j = 1 ! Temporary Variable for second *DO loop
NSEL,S,LOC,X,Sdepth+oH
NSEL,R,LOC,Y,Wwidth/2
NSEL,R,LOC,Z,Wheight/2
CM,fAppNodes,NODE
*DO,TM,TM_START,TM_END,TM_INCR
TIME,TM ! Time value
CMSEL,S,fAppNodes,NODE ! Select last part of beam for load application
F,ALL,FZ,TFORCE(j,1) ! Time-varying force at selected nodes
ALLSEL
SOLVE
j = j+1
*ENDDO
ALLSEL
SAVE
/EOF
Ich hoffe Ihr könnt mir helfen.
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP