Hi,
I am working on heat transfer model, trying to evaluate its transient behaviour (within ABAQUS). Model is kind of complex and so far I've been using constant volume losses in term of heat fluxes. Since these thermal losses are highly dependent on temperature of the previous step, I am trying to write DFLUX subroutine using SOL as variable.
Everything would have been allright if I would have to define only one volume heat flux. But in my model, catch is that I have different components within which the heat is generated and for each different relation to SOL exist.
So I have to define multiple FLUX(1) for volume (JLTYP.EQ.1) and somehow define for which load or which set within abaqus it applies.
So far, I have managed to separate body and surface heat fluxes by using IF LOOP:
IF (JLTYP.EQ.0) THEN
FLUX(1)=0.2*SOL... (code as subject to SOL)
ELSE IF (JLTYP.EQ.1) THEN
FLUX(1)=0.02*SOL (code as subject to SOL)
END IF
Probably, Inside second part of JLTYP.EQ.1 loop the another loop separating different volume fluxes should be integrated. I would appreciate any help or idea how to solve this issue.
Thanks in advance,
Thanks for help. However, I have found other, more sufficient way:
If you use mesh global numbering control and if you set under model attributes "Do not use parts and assemblies in input files", you can explicitly call specific node in dflux by using if loop with NOEL variable. If you are interested in more detailed explanation, you are welcome to ask. :)
BR
- Badges
- Science topic
- Similar topics
- Engineering
- Thermal Engineering
- Thermofluid
- Heat Transfer