Dear FEAP administrator,
I installed feap8.4's fe^2 , and I can run the Iblk3 example successfully.
Now I'm trying to write my own user element, so I'm wondering how to use FE^2
for my own subroutine.
I tried the following inp file, but not success:
**************************************************
noprint
mate,1
user,1
cc la la 0.0 0.0 0.0
la cc la 0.0 0.0 0.0
la la cc 0.0 0.0 0.0
0.0 0.0 0.0 mu 0.0 0.0
0.0 0.0 0.0 0.0 mu 0.0
0.0 0.0 0.0 0.0 0.0 mu
mate,2
user,2
cc la la 0.0 0.0 0.0
la cc la 0.0 0.0 0.0
la la cc 0.0 0.0 0.0
0.0 0.0 0.0 mu 0.0 0.0
0.0 0.0 0.0 0.0 mu 0.0
0.0 0.0 0.0 0.0 0.0 mu
finite
rve,Iinpt
******************************************
for micro's inp, I write the following line:
**************************************
64 27 1 3 3 8
parameter
mi = 3
!noprint
mate,2
...
***************************************
One example I'm trying to do is: for macro scale, I have governing equations for displacement U(linear elastic problem), and
concentration C0(diffusion problem); while for the micro scale, I have a different concentration C1(micro's diffusion problem).
Macro's U and C0 will influce the flux of the micro scale's particle surface.
My idea is, for each gauss point, compute macro scale's U and C0, and then give it as constant to micro diffusion problem, where the micro scale's diffusion equation need U and C0 to compute the flux for its boundary condition. Then we solve C1 and return it to macro scale.
Is it possible to do it in feap? Or any other good ideas?
Thank you!
The attachment is my input file.