using DISP subroutine

[sguaiz]

Hallo Everybody,

I would like to know how to use the subroutine abaqus:

SUBROUTINE DISP(U,KSTEP,KINC,TIME,NODE,NOEL,JDOF,COORDS)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION U(3),TIME(2),COORDS(3)
C

user coding to define U

RETURN
END
I try to put some values like U(1)=0.003 U(2)=0 U(3)=0 just to make it run but it gave some problems in compiling without say me why. I would like to put a function of displacement on certain nodes. Until now I just make a fortran subroutine that can write me the input file and put such boundary condition.

thanks

 

[xerf]

Look in the .msg file.

 

[sguaiz]

Dear xerf,

this is the msg file

1

ABAQUS VERSION 6.5-1 DATE 12-Feb-2007 TIME 23:55:09
For use by Legends Never Die under license from ABAQUS, Inc.

displacement
passo iniziale
STEP 1 INCREMENT 1 STEP TIME 0.00


S T E P 1 S T A T I C A N A L Y S I S


passo iniziale

AUTOMATIC TIME CONTROL WITH -
A SUGGESTED INITIAL TIME INCREMENT OF 1.00
AND A TOTAL TIME PERIOD OF 1.00
THE MINIMUM TIME INCREMENT ALLOWED IS 1.000E-05
THE MAXIMUM TIME INCREMENT ALLOWED IS 1.00

LINEAR EQUATION SOLVER TYPE DIRECT SPARSE

CONVERGENCE TOLERANCE PARAMETERS FOR FORCE
CRITERION FOR RESIDUAL FORCE FOR A NONLINEAR PROBLEM 5.000E-03
CRITERION FOR DISP. CORRECTION IN A NONLINEAR PROBLEM 1.000E-02
INITIAL VALUE OF TIME AVERAGE FORCE 1.000E-02
AVERAGE FORCE IS TIME AVERAGE FORCE
ALTERNATE CRIT. FOR RESIDUAL FORCE FOR A NONLINEAR PROBLEM 2.000E-02
CRITERION FOR ZERO FORCE RELATIVE TO TIME AVRG. FORCE 1.000E-05
CRITERION FOR RESIDUAL FORCE WHEN THERE IS ZERO FLUX 1.000E-05
CRITERION FOR DISP. CORRECTION WHEN THERE IS ZERO FLUX 1.000E-03
CRITERION FOR RESIDUAL FORCE FOR A LINEAR INCREMENT 1.000E-08
FIELD CONVERSION RATIO 1.00
CRITERION FOR ZERO FORCE REL. TO TIME AVRG. MAX. FORCE 1.000E-05
CRITERION FOR ZERO DISP. RELATIVE TO CHARACTERISTIC LENGTH 1.000E-08

CONVERGENCE TOLERANCE PARAMETERS FOR MOMENT
CRITERION FOR RESIDUAL MOMENT FOR A NONLINEAR PROBLEM 5.000E-03
CRITERION FOR ROTATION CORRECTION IN A NONLINEAR PROBLEM 1.000E-02
INITIAL VALUE OF TIME AVERAGE MOMENT 1.000E-02
AVERAGE MOMENT IS TIME AVERAGE MOMENT
ALTERNATE CRIT. FOR RESIDUAL MOMENT FOR A NONLINEAR PROBLEM 2.000E-02
CRITERION FOR ZERO MOMENT RELATIVE TO TIME AVRG. MOMENT 1.000E-05
CRITERION FOR RESIDUAL MOMENT WHEN THERE IS ZERO FLUX 1.000E-05
CRITERION FOR ROTATION CORRECTION WHEN THERE IS ZERO FLUX 1.000E-03
CRITERION FOR RESIDUAL MOMENT FOR A LINEAR INCREMENT 1.000E-08
FIELD CONVERSION RATIO 1.00
CRITERION FOR ZERO MOMENT REL. TO TIME AVRG. MAX. MOMENT 1.000E-05

VOLUMETRIC STRAIN COMPATIBILITY TOLERANCE FOR HYBRID SOLIDS 1.000E-05
AXIAL STRAIN COMPATIBILITY TOLERANCE FOR HYBRID BEAMS 1.000E-05
TRANS. SHEAR STRAIN COMPATIBILITY TOLERANCE FOR HYBRID BEAMS 1.000E-05
SOFT CONTACT CONSTRAINT COMPATIBILITY TOLERANCE FOR P>P0 5.000E-03
SOFT CONTACT CONSTRAINT COMPATIBILITY TOLERANCE FOR P=0.0 0.100
DISPLACEMENT COMPATIBILITY TOLERANCE FOR DCOUP ELEMENTS 1.000E-05
ROTATION COMPATIBILITY TOLERANCE FOR DCOUP ELEMENTS 1.000E-05

TIME INCREMENTATION CONTROL PARAMETERS:
FIRST EQUILIBRIUM ITERATION FOR CONSECUTIVE DIVERGENCE CHECK 4
EQUILIBRIUM ITERATION AT WHICH LOG. CONVERGENCE RATE CHECK BEGINS 8
EQUILIBRIUM ITERATION AFTER WHICH ALTERNATE RESIDUAL IS USED 9
MAXIMUM EQUILIBRIUM ITERATIONS ALLOWED 16
EQUILIBRIUM ITERATION COUNT FOR CUT-BACK IN NEXT INCREMENT 10
MAXIMUM EQUILIB. ITERS IN TWO INCREMENTS FOR TIME INCREMENT INCREASE 4
MAXIMUM ITERATIONS FOR SEVERE DISCONTINUITIES 12
MAXIMUM CUT-BACKS ALLOWED IN AN INCREMENT 5
MAXIMUM DISCON. ITERS IN TWO INCREMENTS FOR TIME INCREMENT INCREASE 6
CUT-BACK FACTOR AFTER DIVERGENCE 0.2500
CUT-BACK FACTOR FOR TOO SLOW CONVERGENCE 0.5000
CUT-BACK FACTOR AFTER TOO MANY EQUILIBRIUM ITERATIONS 0.7500
CUT-BACK FACTOR AFTER TOO MANY SEVERE DISCONTINUITY ITERATIONS 0.2500
CUT-BACK FACTOR AFTER PROBLEMS IN ELEMENT ASSEMBLY 0.2500
INCREASE FACTOR AFTER TWO INCREMENTS THAT CONVERGE QUICKLY 1.500
MAX. TIME INCREMENT INCREASE FACTOR ALLOWED 1.500
MAX. TIME INCREMENT INCREASE FACTOR ALLOWED (DYNAMICS) 1.250
MAX. TIME INCREMENT INCREASE FACTOR ALLOWED (DIFFUSION) 2.000
MINIMUM TIME INCREMENT RATIO FOR EXTRAPOLATION TO OCCUR 0.1000
MAX. RATIO OF TIME INCREMENT TO STABILITY LIMIT 1.000
FRACTION OF STABILITY LIMIT FOR NEW TIME INCREMENT 0.9500
AUTOMATIC TOLERANCES FOR OVERCLOSURE AND SEPARATION
PRESSURE ARE SUPPRESSED
GLOBAL STABILIZATION CONTROL IS NOT USED
FRICTION IS INCLUDED IN INCREMENT THAT THE CONTACT POINT CLOSES

PRINT OF INCREMENT NUMBER, TIME, ETC., EVERY 1 INCREMENTS

THE MAXIMUM NUMBER OF INCREMENTS IN THIS STEP IS 100

LINEAR EXTRAPOLATION WILL BE USED

CHARACTERISTIC ELEMENT LENGTH 2.50

DETAILED OUTPUT OF DIAGNOSTICS TO DATABASE REQUESTED

PRINT OF INCREMENT NUMBER, TIME, ETC., TO THE MESSAGE FILE EVERY 1 INCREMENTS

EQUATIONS ARE BEING REORDERED TO MINIMIZE WAVEFRONT

COLLECTING MODEL CONSTRAINT INFORMATION FOR OVERCONSTRAINT CHECKS

COLLECTING STEP CONSTRAINT INFORMATION FOR OVERCONSTRAINT CHECKS


INCREMENT 1 STARTS. ATTEMPT NUMBER 1, TIME INCREMENT 1.00


***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE
PART-1-1.647 D.O.F. 3 RATIO = 2.38309E+013.

EQUILIBRIUM ITERATION 1

AVERAGE FORCE 0.00 TIME AVG. FORCE 1.000E-02
ALL FORCE RESIDUALS ARE ZERO
ALL DISP. INCREMENTS ARE ZERO
ALL CORRECTION TO DISP. ARE ZERO
THERE IS ZERO FORCE EVERYWHERE

AVERAGE MOMENT 0.00 TIME AVG. MOMENT 2.500E-02
(FROM FORCE FIELD)
ALL MOMENT RESIDUALS ARE ZERO
ALL ROTATION INCREMENTS ARE ZERO
ALL CORRECTION TO ROTATION ARE ZERO
THERE IS ZERO MOMENT EVERYWHERE

ITERATION SUMMARY FOR THE INCREMENT: 1 TOTAL ITERATIONS, OF WHICH
0 ARE SEVERE DISCONTINUITY ITERATIONS AND 1 ARE EQUILIBRIUM ITERATIONS.

TIME INCREMENT COMPLETED 1.00 , FRACTION OF STEP COMPLETED 1.00
STEP TIME COMPLETED 1.00 , TOTAL TIME COMPLETED 1.00


THE ANALYSIS HAS BEEN COMPLETED



ANALYSIS SUMMARY:
TOTAL OF 1 INCREMENTS
0 CUTBACKS IN AUTOMATIC INCREMENTATION
1 ITERATIONS INCLUDING CONTACT ITERATIONS IF PRESENT
1 PASSES THROUGH THE EQUATION SOLVER OF WHICH
1 INVOLVE MATRIX DECOMPOSITION, INCLUDING
0 DECOMPOSITION(S) OF THE MASS MATRIX
1 REORDERING OF EQUATIONS TO MINIMIZE WAVEFRONT
0 ADDITIONAL RESIDUAL EVALUATIONS FOR LINE SEARCHES
0 ADDITIONAL OPERATOR EVALUATIONS FOR LINE SEARCHES
0 WARNING MESSAGES DURING USER INPUT PROCESSING
1 WARNING MESSAGES DURING ANALYSIS
1 ANALYSIS WARNINGS ARE NUMERICAL PROBLEM MESSAGES
0 ANALYSIS WARNINGS ARE NEGATIVE EIGENVALUE MESSAGES
0 ERROR MESSAGES



JOB TIME SUMMARY
USER TIME (SEC) = 1.0000
SYSTEM TIME (SEC) = 0.20000
TOTAL CPU TIME (SEC) = 1.2000
WALLCLOCK TIME (SEC) = 3






maybe I put all the displacement in all nodes, in case how can I put the displacement only in certain nodes?

thanks again

 

[johnhors]

Dear All,

Sguaiz is using a cracked version of Abaqus !

 

[sguaiz]

why are you saying that?

 

[johnhors]

Legends Never Die - is a well known team of software crackers

 

[sguaiz]

what's happen now?

I order the license but it will come in april, but I have to show results with abaqus before that.

 

[bassmanjax]

Squaiz

Is this a student project?

Regards

Martin

 

[sguaiz]

Dear Martin,

yes it is a student project but I will have on april the license as I say. the funding come from the industry so they ask me to do some work with abaqus and until april I cannot stay and doing nothing, otherwise I will have to looking for another job.
what should I do? leave the forum?

Regards,
Sguaiz

 

[bassmanjax]

The forum rules clearly state that student posting is not permitted;

http://www.eng-tips.com/market.cfm

Not only that, but you are clearly using an illegal copy of ABAQUS - that won't endear you to other forum members! I suggest that you get your legal copy of ABAQUS sorted quickly....