单元在整体坐标系下的刚度矩阵

import java.util.Scanner;

public class jbcs{
	static int N;//单元数
	static Double EI;//抗弯刚度
	static Double EA;//轴向刚度
	static Double α;//转角
	static Double l;//杆件长度
	static Double cos;
	static Double sin;
	static Double[][] T=new Double[6][6];//坐标转换转换
	static Double[][] k_ii;
	static int n=100;//结点位移数目
	Double[][]  KK = new Double[n][n];//单元贡献矩阵
	Double[][] K_z =new Double[n][n];//整体刚度矩阵
	Double[][] K_z1=new Double[n][n];
	Double[][] K= new Double[6][6];//整体坐标系下的单元刚度矩阵
	static int[] k= {0,0,0,0,0,0};//定位向量矩阵
	public  void kk(Scanner s) {
		jbcs jb=new jbcs();
		Double[][] k_ii= {
  {EA/l*Math.pow(cos, 2)+12*EI/Math.pow(l, 3)*Math.pow(sin, 2),(EA/l-12*EI/Math.pow(l, 3))*sin*cos,-6*EI/Math.pow(l, 2)*sin},
				{(EA/l-12*EI/Math.pow(l, 3))*sin*cos,EA/l*Math.pow(sin, 2)+12*EI/Math.pow(l, 3)*Math.pow(cos, 2),6*EI/Math.pow(l, 2)*cos},
				{-6*EI/Math.pow(l, 2)*sin,6*EI/Math.pow(l, 2)*cos,4*EI/l}};