MAPLE can do it better

It is nice. Now I can effectively use maple to formulate my coefficinet matrix for FEM.
It took me less than few hours to modify a 2D code that I have written for elasticity for Piezoelectric coupling. I used maple and just add another degree of freedom to FEM code and it was done. The Next step will be writing the code for polar coordinate. This will be tricky since obviously in polar coordinate the coordinate variables enter the K matrix. I wish myself luck.
I was thinking to write the code using UEL for ABAQUS too. I guess it will be a very idea since I will learn the UEL too. But, I am planning to finish what I am doing now ASAP and write a report for it.

The Polar Coordinate

I am trying to model the behaviur of a telescopic actuator today. I found out that I can do it with a 2D model and only 3 degrees of freedom. It will be a 2D mesh with rectangular lagrangian elemtns. But, the problem right now is the weak form that comes from equilibrium equation. I had trouble working with polar coordinate. I am going to learn how to work with radial coordinate that comes into almost every equilibrium equation.
The Constitutive model wont be hard and I am sure that I can use the same thing as Cartesian for this. The main problem will be the equilibrium equation.