aa (sys2.r-r0)/sys2.r "" bb (r1/(r1-r0))^2 "" Ccl_rrrr (lambda+2*mu)*aa "Elastic Moduli in cylindrical coordinates" Ccl_tttt (lambda+2*mu)/aa "" Ccl_rrtt lambda "" Ccl_ttrr lambda "" Ccl_rtrt mu*aa "" Ccl_rttr mu "" Ccl_trrt mu "" Ccl_trtr mu/aa "" Ccl_1111 Ccl_rrrr*(cos(sys2.phi))^4+(Ccl_rrtt+Ccl_rtrt+Ccl_rttr+Ccl_trrt+Ccl_trtr+Ccl_ttrr)*(sin(sys2.phi))^2*(cos(sys2.phi))^2+Ccl_tttt*(sin(sys2.phi))^4 "Elastic Moduli in Cartesian coordinates" Ccl_1112 (cos(sys2.phi))^3*sin(sys2.phi)*(Ccl_rrrr-Ccl_rrtt-Ccl_rtrt-Ccl_trrt)+(sin(sys2.phi))^3*cos(sys2.phi)*(Ccl_rttr+Ccl_trtr+Ccl_ttrr-Ccl_tttt) "" Ccl_1121 (cos(sys2.phi))^3*sin(sys2.phi)*(Ccl_rrrr-Ccl_rrtt-Ccl_rttr-Ccl_trtr)+(sin(sys2.phi))^3*cos(sys2.phi)*(Ccl_rtrt+Ccl_ttrr+Ccl_trrt-Ccl_tttt) "" Ccl_1122 Ccl_rrtt*(cos(sys2.phi))^4+(Ccl_rrrr-Ccl_rtrt-Ccl_rttr-Ccl_trrt-Ccl_trtr+Ccl_tttt)*(sin(sys2.phi))^2*(cos(sys2.phi))^2+Ccl_ttrr*(sin(sys2.phi))^4 "" Ccl_1211 (cos(sys2.phi))^3*sin(sys2.phi)*(Ccl_rrrr-Ccl_rtrt-Ccl_rttr-Ccl_ttrr)+(sin(sys2.phi))^3*cos(sys2.phi)*(Ccl_rrtt+Ccl_trrt+Ccl_trtr-Ccl_tttt) "" Ccl_1212 Ccl_rtrt*(cos(sys2.phi))^4+(Ccl_rrrr+Ccl_tttt-Ccl_rrtt-Ccl_rttr-Ccl_trrt-Ccl_ttrr)*(sin(sys2.phi))^2*(cos(sys2.phi))^2+Ccl_trtr*(sin(sys2.phi))^4 "" Ccl_1221 Ccl_rttr*(cos(sys2.phi))^4+(Ccl_rrrr+Ccl_tttt-Ccl_rrtt-Ccl_rtrt-Ccl_trtr-Ccl_ttrr)*(sin(sys2.phi))^2*(cos(sys2.phi))^2+Ccl_trrt*(sin(sys2.phi))^4 "" Ccl_1222 (cos(sys2.phi))^3*sin(sys2.phi)*(-Ccl_tttt+Ccl_rrtt+Ccl_rtrt+Ccl_rttr)+(sin(sys2.phi))^3*cos(sys2.phi)*(Ccl_rrrr-Ccl_trtr-Ccl_ttrr-Ccl_trrt) "" Ccl_2111 (cos(sys2.phi))^3*sin(sys2.phi)*(Ccl_rrrr-Ccl_trtr-Ccl_ttrr-Ccl_trrt)+(sin(sys2.phi))^3*cos(sys2.phi)*(-Ccl_tttt+Ccl_rrtt+Ccl_rtrt+Ccl_rttr) "" Ccl_2112 Ccl_trrt*(cos(sys2.phi))^4+(Ccl_rrrr+Ccl_tttt-Ccl_rrtt-Ccl_rtrt-Ccl_trtr-Ccl_ttrr)*(sin(sys2.phi))^2*(cos(sys2.phi))^2+Ccl_rttr*(sin(sys2.phi))^4 "" Ccl_2121 Ccl_trtr*(cos(sys2.phi))^4+(Ccl_rrrr+Ccl_tttt-Ccl_rrtt-Ccl_rttr-Ccl_trrt-Ccl_ttrr)*(sin(sys2.phi))^2*(cos(sys2.phi))^2+Ccl_rtrt*(sin(sys2.phi))^4 "" Ccl_2122 (cos(sys2.phi))^3*sin(sys2.phi)*(Ccl_rrtt+Ccl_trrt+Ccl_trtr-Ccl_tttt)+(sin(sys2.phi))^3*cos(sys2.phi)*(Ccl_rrrr-Ccl_rtrt-Ccl_rttr-Ccl_ttrr) "" Ccl_2211 Ccl_ttrr*(cos(sys2.phi))^4+(Ccl_rrrr-Ccl_rtrt-Ccl_rttr-Ccl_trrt-Ccl_trtr+Ccl_tttt)*(sin(sys2.phi))^2*(cos(sys2.phi))^2+Ccl_rrtt*(sin(sys2.phi))^4 "" Ccl_2212 (cos(sys2.phi))^3*sin(sys2.phi)*(Ccl_rtrt+Ccl_ttrr+Ccl_trrt-Ccl_tttt)+(sin(sys2.phi))^3*cos(sys2.phi)*(Ccl_rrrr-Ccl_rrtt-Ccl_rttr-Ccl_trtr) "" Ccl_2221 (cos(sys2.phi))^3*sin(sys2.phi)*(Ccl_rttr+Ccl_trtr+Ccl_ttrr-Ccl_tttt)+(sin(sys2.phi))^3*cos(sys2.phi)*(Ccl_rrrr-Ccl_rrtt-Ccl_rtrt-Ccl_trrt) "" Ccl_2222 Ccl_tttt*(cos(sys2.phi))^4+(Ccl_rrtt+Ccl_rtrt+Ccl_rttr+Ccl_trrt+Ccl_trtr+Ccl_ttrr)*(sin(sys2.phi))^2*(cos(sys2.phi))^2+Ccl_rrrr*(sin(sys2.phi))^4 "" rhocl aa*bb*rho "Density"